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THE MECHANIZATION AND AUTOMATION OF PRODUCTION- 
THE GENERAL LINE OF DEVELOPMENT OF INDUSTRY 


Results of the June Plenum of the Central Committee 


of the Communist Party of the Soviet Union 


On July 13-16, 1960, a Plenum of the Central 
Committee of the Communist Party of the Soviet Union 
was held, at which great attention was devoted to prob- 
lems in the automation of industry. 

The Plenum summed up the work of industry 
and transportation done in the year and a half since the 
Twenty-first Congress of the Communist Party of the 
Soviet Union, The basic conclusion of the Plenum was 
that the grand program for the creation of the material 
and technical basis of Communism adopted by the 
Twentv-first Congress of the Communist Party of the So- 
viet Union is being successfully accomplished, 

Work upon the introduction of progressive tech- 
nological processes, on the mechanization and automa- 
tion of industry, the general line of development of in- 
dustry, is expanding. 

Technological progress in all branches of the 
national economy is inseparably connected with success 
in Soviet sciences, In the area of automatic control our 
science has substantial attainments and, in a number of 
directions, is outstripping the science of more developed 
capitalistic countries, As a result of the theoretical and 
experimental work of Soviet scientists, there have been 
successful planning and building of many complex sys- 
tems of automation and remote control, the automatic 
control systems of cosmic rockets, for example, which 
assure precise escape of artificial satellites of the earth 
and sun into assigned trajectories, or automatic control 
systems for atomic reactors, metallurgical and other 
complex aggregates, 

All these tasks can be solved only on the basis 
of development of the theory of automatic control and 
the creation of automatic systems of high quality. 

Soviet scientists have put forward a number of 
new ideas which already have found application in the 
techniques of automatic control, Results have been ob- 
tained which lay a basis for the theory of discrete auto- 
matic systems and represent a fundamental contribution 
to the theory of automatic control, including relay, im- 
pulse and digital systems of automation,which are exten- 
sively applied in engineering. A number of the princi- 
ples revealed, to which automatic systems are subject, 
and also the methods developed for the investigation and 
calculation of these systems, are being extensively util- 


ized by scientific research and experimental design or- 
ganizations in creating systems of automatic control. 
Important results have been obtained in the area of non- 
linear systems of control, the theory of combined and 
self-contained systems, and in the theory of invariabili- 
ty. A beginning has been made in the development of 
automation in the designing of technical equipment. A 
firm theoretical foundation has been created for the 
building of optimal and self-adjusting systems of auto- 
matic control which establish the most rapid and econ- 
omical system of flow of a technological process, New 
principles of research have been worked out regarding 
semiconductors, magnets and pneumatic elements for 
the building of reliable automatic systems. 

The First Congress of the International Federa- 
tion on Automatic Control has been held with great suc- 
cess, More than a thousand scientists and specialists par- 
ticipated in it, including more than seven hundred dele- 
gates from foreign countries. 


The reports presented to the Congress, and the 
discussions which revolved about them, showed that So- 
viet scientists are working in the most up-to-date direc- 
tions in automation, Many very eminent foreign scien- 
tists unanimously noted the high level of Soviet science 
in the area of automation. 


It should be mentioned, however, that, side by 
side with great success of our scientists in the field of 
automation, the mechanization and automation of pro- 
duction are still being accomplished slowly. The Plenum 
of the Central Committee of the Communist Party of the 
Soviet Union emphasized that solving the fundamental 
problem of the current Seven-Year Plan —to gain the 
maximum amount of time in the peaceful competition 
of socialism with capitalism— is possible only through a 
universal raising of the tempo of technical progress, and 
on one basis — increase in the productivity of labor. 


The July Plenum of the Central Committee of 
the Communist Party of the Soviet Union bound the Par- 
ty, the Soviet and economic organizations to work out 
and put into effect measures directed toward the most 
rapid introduction into production of the attainments of 
science and technique, toward the creation and utiliza- 
tion of the output of the most improved machines, 
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equipment anu apparatus for complex mechanization 
and automation of productive processes, 

Decisive importance in accelerating the tempo 
of technical progress lies in advancing the growth of ma- 
chine-tool construction, the expansion of production and 
raising the technical level of machines, apparatus, in- 
struments and means of mechanization and automation. 
In the light of these tasks, in spite of the success obtain- 
ed, it is impossible to be satisfied with the condition of 
scientific research in the area of automation, It is ne- 
cessary to expand by several times the front of scientific 
research, both as regards problems whose solution will 
permit obtaining in the very near future a great econom- 
ic effect, and as regards long-range problems of the 
growth of automation many years from now. 

It is necessary to expand much more extensive- 
ly the fundamental theoretical and research work, with- 
out which a valuable solution of the problems of auto- 
mation is unthinkable, 

It is necessary to put into effect more rapidly 
the practice, recommended by the Plenum of the Cen- 
tral Committee of the Communist Party of the Soviet 
Union, of cooperative creative participation of scholars, 
engineer-technologists, designers, builders, and factory 
workers in the creation of progressive technological pro- 
cesses and new technique in all stages of their develop- 
ment from laboratory investigation to introduction into 
production, 

In recent years there has been a substantial 
change in the character of the problems facing science 
in the area of automation, Present-day complex systems 
of automatic control can no longer be built on the basis 
of mere engineering intuition or of experience. Without 
guiding principles basea on profound theoretical assump- 
tions, it is impossible to determine the most advantage- 
ous structures for control systems. Therefore, it is neces- 
sary to accelerate the elimination of gaps in the theory 
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of automatic control by concentrating our efforts on the 
solution of the main problems, such as problems in the 
creation of new control systems by complex aggregates 
with an abundance of elements and connections between 
them, 


Side by side with the further development of 
essential divisions in the theory of automatic control, 
special attention should be devoted to the creation of a 
theory of self-adjusting systems and systems of optimi- 
ization and a theory of modeling and the algorithmiza- 
tion of technological processes, 

It is necessary to develop the method and ap- 
paratus for automatic investigation of technological and 
industrial processes as the objects of automation, 

It is necessary to have considerable strengthen- 
ing of research work on raising the reliability of control 
apparatus, In the area of problems in the building of 
technical means of automation, this task is a decisive 
one, since, if there will be no automatic equipment of 
high reliability, the tasks of building completely auto- 
mated productive complexes can not be solved in prac- 
tice, 

The resolution of the Plenum of the Central 
Committee of the Communist Party of the Soviet Union, 
"On the course of fulfillment of decisions of the Twenty- 
first Congress of the Communist Party of the Soviet Un- 
ion regarding development of industry and transport, 
the introduction ot the latest attainments of science and 
technique into production,” has an exceptionally impor- 
tant significance for the creation of the material and 
technical basis of communism in our country. 

Soviet specialists im automation are applying 
every effort, all their knowledge, intelligence, and 
energy to bring to life the historic decision of our Party 
regarding further acceleration of the building of a com- 
munist society, 








AUTOMATIC OPERATING MODE OPTIMIZATION 
OF A CLASS OF SYSTEMS WITH RESPECT 
TO STATISTICAL PERFORMANCE CRITERIA* 


S. A. Doganovskii 


Moscow 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 8, pp. 1105-1114, 


August, 1960 
Original article submitted March 22, 1960 


This article is concerned with problems in automatic optimization of the operating mode of a system sub- 
ject to random disturbances. The design of a computer whose operation is based on assigned performance criteria 
is substantiated. Results obtained in investigating the optimization process are given. 


1, Formulation of the Problem 
One of the most important problems in com- 
plex automation consists in securing such a control of 
industrial devices that they operate under the most ad- 
vantageous (optimum) conditions. This problem is now 
being solved, due to the widespread introduction of vari- 
ous modern computer devices in control systems and the 
provision of self-adjusting automatic control systems. 
As is known [1], automatic optimization sys- 
tems, which select the system optimum operating mode, 
represent an important class of self-adjusting systems. 
The optimum operating mode of a system is 
determined in correspondence with the performance cri- 
terion, which is assigned beforehand on the basis of tech- 
nological requirements imposed on the quantity to be 
optimized in concrete problems of industrial process 
(device) control, The formulation of the required and 
the most suitable performance criterion is, in each case, 
a separate, and, quite often,a very complicated problem. 
In a number of industrial systems to be con- 
trolled, where random disturbances play an important 
tole, it is advisable to accept a certain average value 
for the deviation of the quantity x to be controlled as 
the performance criterion ¢, for instance, 
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where a; are certain weighting factors of the system, 
and Fj are disturbances acting on the system. Hereaf- 
ter, the optimization of the system operating mode will 
designate the minimization of the ¢ quantity, which is 
determined on the basis of an accepted performance cri- 
terion, 

Such statistical processing of data on the pro- 
cess (system to be controlled) state must be performed 
by the computer, which we shall call the unit for form- 
ing the quantity to be minimizedt (UFQM) [4]. The 
UFQM output signal ¢ determines the variable (or a 
number of such variables ¢j) to be minimized in the 
process of optimizing the system operating mode, Then 
the 9; quantities, together with the values H; of the re- 
strictions for the entire system which are to be taken in- 
to account [2], are fed to the inputs of an automatic op- 
timizer, which searches for the necessary system para- 
meters Z; for which the minimization of the UFQM out- 
put signal is secured, 


* The basic principles used in this article were dis- 
cussed on June 10, 1959 at the IAT AN SSSR (Institute of 
Automation and Remote Control, Academy of Sciences, 
USSR) seminar on the theory of optimum, self-adjusting, 
and interdependent automatic control systems. 

tT In the particular case of correction with respect to 
the minimum of the mutual correlation function, the 
forming of ¢ is performed by a correlator-type comput- 
er. 


779 











Fig. 


The control device of such an automatic con- 
trol system, which was developed at IAT AN SSSR dur- 
ing the years 1957-1959, was designated as the comput- 
er for optimizing the system operating mode (COSOM) 
(Fig. 1a). The automatic control system, the block dia- 
gram of which is shown in Fig. 1b, has the following 
characteristics. In the first place, the formation of the 
control signal ¢; is here effected after preliminary sta- 
tistical processing of the system state data, which is se- 
cured by UFQM, Therefore, the effect of random fac- 
tors (disturbances) acting on the optimizer A, which 
is the basic component in the control of the system, is 
considerably reduced, In the second place, the quantity 
to be minimized is not the controlled variable x itself, 
but the quantity ¢, which is formed on the basis of the 
performance criterion, which determines the optimum 
operating mode of the system, Therefore, the system 
characteristic with respect to ¢ is always at the opti- 
mum, even in the case where the characteristics of the 
system to be controlled are linear. Papers [2, 3, 4] are 
devoted to problems in the design of such a class of au- 
tomatic control systems. 

The main purpose of the present paper is the 
investigation of the process of the system's approach to 
its optimum operating mode under the action of random 
disturbances. In the general case, the process of the sys- 
tem motion as represented in Fig. 1b is an unsteady-state 
process even in the case of steady-state random distur- 
bances F;. , 

The consideration will be restricted here solely 
to the case where the system operating mode is 
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searched for with respect to one variable of the system 
to be controlled, without taking into account possible 
restrictions of the control system variables with respect 
to the integral of the square of the controlled variable 
x for a steady-state random input disturbance F. The 
case of the asymptotically steady-state (quasi-steady- 
state) process of the system motion is considered. The 
obtained results are then rendered more accurate by es- 
timation of errors to be found, which is connected with 
the limitation of the UFQM averaging time interval in 
actual self-adjusting systems. 


2. Equations of the System and the 
Investigation Method 

For an illustration of the method used in inves- 
tigating the system whose diagram is shown in Fig. 1b, 
consider an actual system O°, which is a system for com- 
pensating disturbances F, acting at its input (Fig. 2). 
The structural diagram of such a system is typical for 
control systems which respond to disturbances [5, 6]. 
Here, disturbance F is fed directly to one of the inputs 
of the system O to be controlled and, moreover, to its 
second input through compensating device ky. The ad- 
justment of the compensating device secures compensa- 
tion of disturbance F at the output of the system to be 
controlled, i.e., in such a system, conditions under 
which the controlled variable x is independent (invari- 
ant) of disturbance F can be secured. If the character- 
istics of the system to be controlled change and, conse- 
quently, if the conditions for the controlled variable in- 
variance are disturbed, COSOM must act, in the general 














case, on the parameter of the linear, as well as the non- 
linear section of the compensating device, thereby cor- 
recting its operation. This correction of the compensat- 
ing device parameters is performed in accordance with 
the beforehand-assigned performance criterion on the 
basis of statistical processing of the system state data. 
The tuning of the compensating device is completed af- 
ter the conditions for the controlled variable invariance 
are satisfied, i.e., after the system optimum operating 
mode has been secured, Thus, in the system shown in 
Fig. 2, the conditions for the controlled variable invari- 
ance are automatically secured by COSOM, which leads 
to the most advantageous (optimum) operating mode of 
the entire system. 
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Let us consider the case where the system to be 
controlled is an inertialess linear system whose channels 
have transfer coefficients ky and kp, respectively, and 
the compensating device is an amplifier with a variable 
amplification factor Ky. 

For the performance criterion of the system op- 
timum operating mode, we shall take the well-known 
criterion of minimum root-mean-square error (RMSE). 
In practice, it is more suitable to take the minimum of 
the controlled variable mean square, i.e., 


{ 


9 =lim 7 \ [x(t)}Pde. (1) 


ou 


Let us consider the search for the system opti- 
mum operating mode without taking into account the 
restrictions imposed on the system variables for the class 
of steady-state random input disturbances F with the 
spectral density Sp (w). The system control device, CO- 
SOM, operates in the following manner [4]. Let there be 
a certain initial value of the amplification factor ky(0) 
of the compensating device for which the invariance 
conditions in the systes: are not satisfied, i.e., kg 0) # 
k;/kp, Consequently, there will be random deviations of 
the controlled variable x(t) at the system output. At the 
beginning, COSOM forms the ¢, quantity in UFQM ac- 


cording to the assigned performance criterion (1) and 
memorizes %;, Then COSOM assigns a trial step by 
changing the amplification factor kg(0) of the compen- 
sating device by Akg; in this, the ¢= 9, +4¢ quantity is 
formed again, and the quantity 4¢ is determined, which 
is fed to the input of the final (integrating) optimizer 
device, 2nd which changes the amplification factor ky of 
the compensating device, 

For the shaping of ¢ signals, UFQM is provided 
with a generating unit, which is used as a quadrator, and 
an averaging unit, which consists of a pulse filter. The 
latter represents an integrator with the cancellation of 
readings after integration time T. The necessary divi- 
sion of the integration result by T is secured by a suit- 
able choice of the integrator scale factor. 

From a practical point of view, ¢ cannot be 
determined over an infinite integration time interval. 
Therefore, the determination quantity ¢* [7], which is 
given by the following expression, is formed in UFQM: 


T 
g° = + [x(t)}* dt. (2) 


0 
Assuming that the formed signals ¢ and ¢* are 
close to each other, we arrive at the system asymptotic 
steady-state operating mode. Then the ¢ quantity can 
be expressed in terms of the spectral density S-(w) of 
the input disturbance by the relation 


oo 


= — | Sr(w)| K (ja)|*do, (3) 
0 


where | K(jw)| is the modulus of the transfer coefficient 
of the system to be controlled, 


For the case of an inertialess linear system, expres- 
sion (3) will be written in the following form: 


co 


p= 2 Se (w) (hs — hgka)* den = D (ky — ake)? (4 
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where D designates the input disturbance dispersion 


"0 


D = —\Sp(w) do. tm 
0 


Since ¢ is the quantity to be minimized in the 
system, Eq. (4) determines the system characteristic as 
the object of optimum (extremal) control (Fig. 3), The 
amplification factor ky of the compensating device is 
the argument of this characteristic, and the transfer co- 
efficients ky and kz of the system to be controlled and 
the dispersion D of the input disturbance are the para- 
meters of the characteristic. A change in the system 
parameters will cause a shift of the characteristic extre- 
mum with respect to the ky axis and will change the 
slope of its branches. The minimum of the ¢ quantity 
will be secured for a given value of the compensating 
device amplification factor ky for which the following 
condition for the controlled variable invariance is satis- 
fied: 


1 
he; 


ky= (6) 


For ky > 0 and kp> 0, kg cannot be determined 
by values in the left half-plane in Fig. 3, which is 
marked by the shaded boundary of the region where 
changes of the compensating device amplification fac- 
tor are allowable, 

Figure 3 also shows the magnitude of the trial 
step Ak», which is assigned by the automatic optimizer 
in searching for the necessary value of the amplification 
factor kg which will satisfy Eq. (6). Thus, the forma- 
tion of the quantity ¢ to be minimized on the basis of a 
given performance criterion makes it possible to reduce 
the linear system O characteristic to the extremal value 
(Fig. 3). 

Figure 4 illustrates the structure of a system 
subject to extremal control with the characteristic of 
Fig. 3. 

Here, Ky, Kz, Ks, and Ky are contacts of the op- 
timizer control device A, MD is the optimizer memory 
device, and ID is the optimizer final (integrating) de- 
vice, A detailed description of such an optimizer system 
is given in [4]. 

Let us find the equations of motion for the sys- 
tem and determine the number of search steps n for which 
the system will be within the Ak)/2 zone in the neighbor- 
hood of the extremum point of the characteristic of the 
system to be controlled. The increment in the output 
voltage of the automatic optimizer final device will be 


determined by the expression 
Tos 
k k 
b= \ Agdt=—ae- Tos Agy (1 
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where Ag; is the function increment at the ith point, 
Tgs is the optimizer operating step time, kgq is the scale 
factor of the optimizer final device, RC is the time con- 
stant of the optimizer final device, and z; is the output 
voltage of the optimizer final device. 
By substituting in Eq. (4) 

ky om koko = Z, (8) 
we obtain the following expression for the characteristic 
of the system to be controlled: 


o = Dz’. 


Consequently, for small values of Az, 


Ao = 22,DAz. (10) 


By substituting in Eq. (7) the value of Ag; 
from Eq. (10), we obtain the equation for the output 
voltage increment in the optimizer final device: 


k : ‘ 
82, alee aS 7 os 2Dz,Az = 62; = BDz, (11) 
where 
k 
B= 2-N-T og As. (12) 


As it is assumed that the disturbance F which 
acts on the system belongs to the class of steady-state 
random disturbances, the quantity D is constant and, 
consequently, 9 is also a constant quantity. 

Starting with the initial value of kp(0), the val- 
ues of 6z; will be determined by expressions of the fol- 
lowing form: 

82; = 62; = 6(zi_, — 5z;_,) = 


= 0[z5(1 — 6) *—6z, (1—8)'4 - 
= 0z)(1 — 8)’, 


(13) 



































and the value of z, will be determined by the expres- 
sion 


Zn = 2 (1 — 0)". (14) 


By solving simultaneously Eqs. (8) and (14), we ob- 
tain the equation 


ky — kgkg[n] = [ky — kako (0) (1 — 8)” = 


‘ (15) 


= [Gi — ko (0) | (1 — 0)" he. 
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By solving Eq. (15) with respect to k(n], we 
find the system equation of motion as a function of the 
number of steps in searching for the required value of 
the compensating device amplification factor: 


ko{n] = <211 —(1 —8)"] + ko (0) 11 — 81". (6) 


From a consideration of the system equation of 
motion in the adjustment of the required amplification 
factor value (16), it follows that the following inequality 
has to be satisfied if the scanning process is to be conver- 
gents 


| 1-6| <1, (17) 


iJBans 


§< 2. 


(18) 


Figure 5 shows the diagrams of the system mo- 
tion during the adjustment of the required value of the 
compensating device amplification factor kp as functions 
of the argument n for jumplike changes in the k,/k, pa- 
rameters of the system to be controlled. 

Until now, our point of departure was the as- 
sumption that the dispersion D was formally determined 
by averaging the square values of input disturbances F of 
the system to be controlled over an infinite observation 
interval (T-@). However, the integration time interval 
is always finite in practice, and it is convenient to make 
it as small as possible, because, in the first place, the 
actual integrator circuits are imperfect and their errors 
increase with time, and, in the second place, a reduction 
of the scanning time is thereby achieved, Therefore, for 
steady-state random input disturbances, the magnitude 
of dispersion (D.), produced over the finite interval T, 
will also be a random quantity, which will fluctuate 
about its average value D, which is determined for an 
infinite integration time interval. 

The degree of approximation of D by the D+ 
function can be characterized by the magnitude of the 
mean square deviation of these functions: 


op = M [(Dy — DY}, (19) 


where M is the sign of mathematical expectation, ort 


T 
op = Fe \(T — t) 2A? (2) dr, 


0 


(20) 


where 
R(t) = M [F (t) F (¢ 4- +). 


Thus, if we know the input disturbance autocorrela- 
tion function R(T), we can determine o, . 


Since the input disturbance spectral density 
S; (w) is connected with the autocorrelation function by 
Fourier transformation, we obtain the following expres- 
sion: 


0 


0 


Let us consider a case often encountered in prectice, 


where the spectral density of the input disturbance F is 
of the following form: 


So 


Sp (a) = i407? . 
0 


(22) 


Then the equation for determining the root-mean- 
square error ©, will be written thus; 








T rs) 
== V \r-9 \ = dw] de. 
0 


The integral in square brackets in expression (23) 
can be determined as; 


ww t 
\ ilies a ae 
) 4 tw? OT, : (24) 





«0 y 





Fig. 5 





+ A normal probability distribution of the F(t) function 
is assumed. 
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By substituting (24) in (23), after suitable transfor- 
mations, we obtain the equation 





a 
where 
ES anerss, a is —2T/T, 4 5 3 
* (4) = 73 V ae) pe (26) 


Expressions (25) and (26) indicate that the root- 
mean-square error tends to zero for T*> ©. By assigning 
the maximum allowable value 9, ,, So. and To, we 


find the required value of the integration time T by us- 
ing the inequality 


Sp = 39 (Fe) <p me (27) 


The graph of the #(T/T,) function is shown in Fig.6. 


Let us now determine the number of search 
steps for which the system falls within the Ak,/2 zone in 
the neighborhood of the extremum point of the system 
characteristic. 

From Eqs. (13) and (14), we find 


zo(1— On = 5. (28) 














The number of search steps for which the system falls 


within the Ak,/2 zone is determined from Eqs. (28) and 
(12): 
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Fig. 8 


Figure 7 shows the graphs plotted with respect 
to Eq. (29), which indicate the effect of quantities B, 
C, and D on the number of steps required in searching 
for the system optimum operating mode, The quantity 
C shows how many trial steps Ak,/2 are contained in the 
interval between the initial and the required (for secur- 
ing the extremum) values of the compensating device 
amplification factor ko; it is determined by the assigned 
value of Ako, which is the initial deviation of the com- 
pensating device amplification factor from the required 
value, and it does not depend on D + fluctuations. 
Since, as was mentioned above (Fig. 3), changes in D., 
cause changes in the slope of the system characteristic 
branches, the number of search steps for which the sys- 
tem motion will settle in the neighborhood of the extre- 
mum point will also fluctuate about its average value 
M[n]. 

The mathematical expectation of n, which is, 
as can be seen from Eq, (29), a nonlinear function of 
D>, will be determined by the expression 


p(Dy) 


‘ 0° 


—oo 


adDr. (31) 


Here, p(D,,) is the probability density of the ran- 
dom quantity Di. 

The above investigation is theoretically valid 
for the case of an asymptotically steady-state process of 
the system motion toward its optimum operating mode, 
Nevertheless, the investigation method considered is ap- 
plicable to a number of practical problems. 

As an example, Fig. 8 shows a typical oscillo- 
gram of thickness fluctuations of a metal strip at the in- 
put of a continuous-operation cold-rolling mill [8], By 
using Eqs. (19) - (27), a preliminary statistical processing 
of data in this oscillogram will make it possible to find 
the necessary averaging time T of the UFQM integrating 
device, to determine the COSOM parameters with re- 
spect to Eqs. (29)-(31), and to determine the necessary 
number of search steps in designing a self-adjusting sys- 
tem for the automatic control of the strip thickness in 
continuous-operation cold-rolling mills [4]. 


SUMMARY 


1, Automatic »ptimization systems, including the 
computer for optimizing the system operating mode (CO- 
SOM), represent a new and promising direction in the de- 
velopment of modern high-quality automatic control sys- 


tenis. They can be used in various branches of technolo- 
gy, for instance, for the optimum control of welding con- 
ditions in electric tube-welding mills, for the control of 
the strip thickness in continuous-operation cold-rolling 
mills, in designing self-adjusting autopilots, etc. 

2. The performed investigation of the self-adjust- 
ment process in systems of this class makes it possible, 
in a number of practical cases, to determine the neces~- 
sary number of steps in the search for the optimum oper- 
ating mode and to determine the parameters of individual 
COSOM units, 

The author extends his thanks to A, A, Fel’dbaum 
for his valuable advice in the completion of this work, 
and to Yu. P. Leonov for his discussion of the obtained 
results, 
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The design of the two simplest types of extremal systems with inertial linear units is considered. 


In the system of the first type, the scanning process is completed in the minimum time for assigned res- 
trictions of the variable of the system to be controlled. In the system of the second type, the scanning process 
takes place for a minimum value of the dynamic error integral. The diagram of the optimum control equipment 
and the results of experimental investigations of the systems are provided. 


1. Introduction 


Systems which are designed for the automatic de- 
termination of the input in any industrial device for 
which its output, which depends on the input value, as- 
sumes either the minimum or the maximum value are 
presently used to an ever-increasing extent. Such sys- 
tems are called extremal systems, 

In [1], it was shown that extremal systems be- 
long to a separate class of automatic control systems, 
being conceptually different from most other such 





systems, At the same time, extremal and(other) automatic 


control systems have many characteristics in common, 
Therefore, in developing the theory of extremal systems, 
it is often possible to use certain principles and methods 
of the more developed theory of automatic control sys- 
tems, 

The present paper is concerned with the design of 
optimum extremal systems with respect to quickness of 
response and the integral value of the dynamic error. 

It is assumed that the scanning process for finding 
the extremum points occurs in an inertial system with a 
limited value of one of its variables, For automatic 
control systems with optimum quickness of response, this 
problem and its solution are known [2, 3]. 

The basic results obtained in these papers can also 
be used in designing extremal systems with the optimum 
quickness of response. 

In designing optimum extremal systems, a number 
of restrictions are imposed on the characteristics of the 
system where the scanning is to be effected. It is as- 
sumed that no random disturbances affect this system 
during the scanning process and that the system output 
quantity y is related to its input quantity x by a parabol- 
ic dependence of the following form: 


y=a(r— b)+c. (1) 


Quantities c and b can assume different values dur- 
ing the scanning process, The a coefficient is assumed 
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to be constant, although, for certain more complicated 
control systems, this restriction may not be introduced, 
It is assumed that b can change only stepwise in the sys- 
tem to be controlled, 

As a rule, the characteristics of actual industrial 
systems differ from those used in this paper, and, there- 
fore in certain cases, for instance ,whenrandom distur- 
bances of large magnitude are present in the system to 
be controlled, the extremal system may be inoperative, 
The characteristics of the system to be controlled are 
usually approximations of the above-indicated charac- 
teristics, and, therefore, the extremal systems under con- 
sideration will be only approximations of the optimum 
systems, However, even in these cases, although the 
control system described below is not actually an opti- 
mum system, it will still be possible to reduce the scan- 
ning time for finding the extremum, as well as the inte- 
gral dynamic error of the extremal system. 

We shall assume that the system to be controlled, 
where the scanning takes place, consists of two units; the 
linear inertial and the nonlinear inertialess units, The 
linear unit consists of two integrating elements, which 
are connected in series, The linear unit input is restrict~- 
ed with respect to its modulus. The nonlinear unit of 
the system is described by Eq. (1). We shall consider the 
case where it is possible to measure the rate of change 
in the linear unit output and, as a result of this measure- 
ment, obtain, for instance, a voltage proportional to this 
rate, 


2, Extremal System with Optimum 





Quickness of Response 


The structural diagram of the extremal system is 
shown in Fig. 1. The system consists of the scanning ob- 
ject O and the control equipment CE, We shall consider 
the case where the system input quantity x, is limited 
with respect to its modulus: | x, | <E. Let the coordin- 
ate origin be chosen in such a manner that b=0 corres- 
ponds to the extremum of y at a certain instant of time, 














which is taken for zero. Let us assume that, for t= 0, 
the value of x will deviate from the extremum point by 
b. It is required that the system find the extremum 
point in the minimum time. 

In this case, the deviation from the extremum point 
can be determined by the equation 


1 dy 
6 =~s.az° (2) 


It is known[2,3] that, in an automatic control 
system which controls a system that includes two inte- 
grating elements and has an input quantity whose modu- 
lus is restricted, the deviation is eliminated in the min- 
imum time for the following control conditions; 








+ E for m |x| -+ =O > 0, 
Ig = (3) 
—E for 2|2,|+ =ar < 0. 


a 

Here, E designates the maximum allowable value 
of X20 

This control condition can be used in designing an 
extremal system with optimum quickness of response. In 
an extremal system, it is first of all necessary to reduce 
most effectively the value of x; when the system moves 
in the wrong direction. This operation must be effect- 
ed at definite times while the system moves in the cor- 
rect direction, inorder to prevent the transition of the 
extremum point due to the presence of integrating ele- 
ments in the system, In order to determine the switch- 
ing time, according to (3), it is necessary to know the 
value of x, and the deviation of x from the extremum 
point, which is given by (2). 

Let us express b in the following form: 


pa tay _ 1 dy de _ 1 dy 


"ra "“thaa: * 


2a dt 


Hence, by taking into account (3), the control con- 
ditions for the extremal system can be written thus; 




















( 4 E,when |2}| $2 4-52 >0 for 2,<0 
and [at] Get + GE <0 for x, >0, 
r= 
| _. xe,when | 23) 0 YY >0 for 1 >0 
a d 
[ and atl gar tear <0 for 1, <0 
(5) 
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Fig. 1. 


The electric circuit of the equipment securing opti- 
mum control of system O is given in Fig. 2. The equip- 
ment includes a polarized relay P2:, a two-coil polar- 
ized relayP,,a de amplifier A, [4], with diodes D, and 
D, for the shaping of |x,| , and a differentiator, which is 
based on an amplifier Az. For obtaining the |x{(aa,/Ea,)) 
term, a carborundum resistor T with a nearly third- 
power voltampere characteristic is used [5]. The char- 
acteristic of the nonlinear resistor T can be brought clos- 
er to a third-power characteristic by means of resistor 
Rg. If the output y of the system to be controlled chances 
slowly, difficulties are encountered in obtaining the 
dy/dt derivative, The output signal y usually contains 
high-frequency disturbances, which are intensified by 
the differentiator and which suppress the main signal pro- 
portional to the derivative. Therefore, instead of the de- 
rivative, it is advisable to determine the increments 


Ay = Yi — Yo (6) 
for equal time intervals 


At = t; — ty. (7) 


In determining Ay, it is necessary to calculate the 
small difference between two large voltages and also to 
memorize the voltage over the 4t period. This problem 
was solved by means of a capacitor circuit [6]. 


3. Extremal System with the Minimum 
Integral Dynamic Error 


The following quantity will be designated as the in- 
tegral dynamic scanning error: 








T 
[= \(y — %) dt. (8) 


Here, T is the scanning time, and y, is the extremal 
value of y. 
According to (1), 


Ye = ¢. (9) 















































Fig. 2. 
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By substituting (1) and (9) in (8), we obtain 
T 
l = lace -- bp dt. (10) 
0 


For the system considered above (Fig. 1), we shall 
first find the control conditions for which the error will 
be at the minimum: 


I= Imin - 


In automatic control systems, the I,,j, quantity cor- 
responds to the minimum value of the integral of the 
quadratic dynamic error, 

We shall use L. S, Pontryagin’s maximum principle 
[7]. We shall substitute the functions and write the sys- 
tem equation in terms of symbols used in [7]; 


I =z a(x—b)? = ai, 
V aa,x, = Ly, ad V az, =U. 


As a result, we obtain 


——— > ey 


T 

a ye ate” =e ld 

iy = \ aidt = \ 7d, ae in 
0 0 

2, = 2, = f, T, = u(t) = A 


By using the L, S, Pontryagin maximum principle, 
we shall write the equations for the ¥)=¥%9(t), ¥1=¥,(t), 
and ~,=¥,(t) functions: 

. 2 of 
%=—V be (i=0,1, 2), 


Ox; 
a=0 


bo one 0, dy re 2210, be ae oe Pr. 





(12) 


The function H=H(Xp, Xp Xe, Vo, V1, Vo, u) is of the 
following form: 


H = dof? + dif" + def? = dori + dite + dau. (13) 


If u=u(t) represents optimum control, then, ac- 
cording to the maximum principle, the H function must 
attain its maximum with respect to u, when u changes 
within the | u(t) |<1 and 0<t=T intervals and ¥ = 
= const=—a*S0. Hence, for optimum control, it follows 
that 


u=sign pp. (14) 
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Thus, the totality of optimum trajectories must 
satisfy the following system of equations: 


I = 2g, Tg=—Signds, $, = 2a°X, by = — hy. 
(15) 

Thus, in order to secure I,)jp, as well as the mini- 
mum scanning time, the system input whose modulus is 
restricted must assume either the x,=+E or x,=<E val- 
ue, 

Let us consider the motion system in two intervals 
corresponding to x,=+E and x,=-E, The velocity at 
the beginning of the first interval and at the end of the 
second interval is assumed to be equal to zero, In this 
case, the motion times for the first (ty) and the second 
(tg) intervals are equal to each other: 


t= th. (16) 


The dynamic error is 
ty 
I ae \ (aye — b): dt 4 
0 
ty 


Et? (17) 
+a \ (“3 ‘1 —b “PF a,a,Et,t — see) ae = 








2 2 
0 


= Fe ajayaE*t} — 2a,a,abt? + 2b7at,. 


Let us determine t, for which I=I,,jn. By finding 
the partial derivative of I with respect to ty from Eq. 
(17), 


= aia; E* tt — 6a,a,Ebt? + 2b7 = 0 (18) 


we obtain 


gk b 
ty = 1.04 a3a,8 . (19) 


By substituting (19) in (17), we obtain 


1... — 2:772ab* Vo 
min Vaat (20) 


It should be noted that, according to [2], 


° b 
. 2 V wz (21) 


for a similar system for optimum control conditions with 
respect to quickness of response. 

In the last case, for two intervals, the system oper- 
ation will coincide with the extremum point. 








ith 





From a comparison between (19) and (21), it is ob- 
vious that, for control conditions securing a minimum 
value of the dynamic error integral, the retarding time 
is shorter than for optimum control conditions with re- 
spect to time. Consequently, the system must first pass 
the extremum point, then change direction of motion, 
and again pass through the extremum point from the oth- 
er side, etc, 

We shall demonstrate that the system motion during 
the intervals following the two considered above hardly 
affects the value of the dynamic error integral. 

From (19) and (21), we find 


t = 1.04 t,° ‘. (22) 


The system transition beyond the extremum point 


Ab = (1.04t,)* a,a,/ — tia,a,E = 0,08a,a, Et?. (23) 


By using (20), we find the ratio of Impin= min (4b) 
to Gnin = tmin(b): 


~ 


min — (0.002. 
Jenin (24) 


The Tnin quantity can be neglected. Fig. 3 shows 
the dependence of GI on Bty, where 


a,a.k 
any ae (25) 
a,0,8 (26) 
pay 


It is obvious from the graph that, for equal & and 8, 
the dynamic error integral I, for a system with opti- 
mum quickness of response [which corresponds to t,8=1 
according to (21)] differs little from its minimum value, 
which corresponds to ty8= 1.04 according to (20), 


a] 





Fig.3 


By substituting (21) in (17), we find 


Joqg= 9.778 a. (27) 
According to (20), 
fain = 0-772. (28) 


Let us denote the ratio of I, to I oar by 6. From 
(27) and (28), we obtain 


I 
8% = ~paip == 99.23%. (29) 
oar 


Thus, the considered extremal systems, the one with 
the optimum quickness of response and the one with the 
minimum value of the dynamic error integral, are close 
to each other and they have practically the same control 
equipment, 


4. Experimental Investigation Results 

The above optimum control scheme for an extremal 
system was experimentally investigated by me~.s of an 
electronic simulating device. The quadratic character- 
istic of the nonlinear part of the model simulating the 
system to be controlled wasreproduced by means of a car- 
borundum resistor [5]. The oscillogram in Fig, 4 shows 
two scanning processes, In the first section, for t < t,, 
the scanning was effected for the open carborundum re- 


sistor circuit [the term [ |2e| eo*] was not included in 





control conditions (5)J. The scanning system did not 
find the extremum pint, and the amplitude of the sys- 
tem oscillations increased, 


For t=t,, the | 2° | 7 term was introduced in the 
1 


control program, and the system attained the extremum 
point in the minimum time and performed small oscilla- 
tions about this point. 

We also investigated the optimum scanning process 
in an extremal system whose inertial part consisted of an 
inertial element with the time constant T and the ampli- 
fication factor a, and an integrating element with the 
coefficient a2, which were connected in series, The op- 
timum conditions for controlling a regulating system with 
such an inertial part have been considered in [8]. 


jE —— 











Fig.4 


789 











aS hn eS nn 


For the experimental system, these conditions can 
be written in the following form by taking into account 
(4): 


(+-£, when 2aa,7 | zi + a,Ez, In 

















mec | and 2aa,T | a, + a,Ez, In 
ra ! — E, when 2aa,T | a, + a,Ezx, In 
| and 2aa,7'||2} + a,£2, In( 
a » EPL 
t t 
PAA 2 
- Fig. 5 b 


The oscillogram of the scanning process in an ex- 
tremal system where nonlinear velocity feedback was 
not present in the control equipment is shown in Fig. 5a. 
The oscillogram of the scanning process in the same sys- 
tem with the optimum control conditions (30) is shown 
in Fig. 5b, 

As in the above-considered case, the system was 
able to find the extremum point in the minimum time 
(for a given restriction of x2), and it performed small 
oscillations about this point. 

K. B. Norkin and V. A. Yakovlev collaborated in 
the experimental investigation of this system. 
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The problem of optimization (in the specified meaning of the term) is considered in tracking a trajectory 
of a body in accordance with the rule of proportional search [1], The control function related to the position of 
the direction control is given on an interval of the real axis. The method of parametric representation of the 
control function [2, 3] is used in solving the problem; this enables one to use the methods of classical calculus 
of variations. The optimum control is found and its sufficient conditions are analyzed. 


1. Statement of the Variational 





Problem in the General Case 





The equations of motion of a point A* attacking 
goal Bt in accordance with the rule of proportional 
search, are as follows [1}: 


oD, = v, sin (p,-— W,) — vsin(g — ®,), 
- =F(t) d, (1,1) 
p = ¥, COS (e-— ,) — v cos(¢ — @,), 


where p is the distance from the point A to the goal B; 

¢ is the climbing angle of the trajectory of point A; ¢, 

is the trajectory angle of the goal B (¢, = const)— putting 
%> 0 does not reduce the generality of the case; ©. is 
the sighting angle (Fig. 1); v and y are the magnitudes 
of the velocities of the point A and the goal B,respective- 
ly (they are assumed constant, and y.< v); &t) is the con- 
trol function connected with the autopilot and with the 
position of the direction control. 

When p< 0, i.e., when % sin (Q@— @-)— v sin (¢-@)= 
=0, it is the case of the so-called “tracking triangle,” 
i.e., a relationship between ¢ and ® such that it ensures 
the interception of the goal. At the initial moment of 
time f9g@)= 9() and o@ = &.(1) do not usually 
satisfy the above relation, and thus arises the optimum 
(in the specified meaning) control problem of getting 
point A on the tracking triangle. 


We shall regard the following as the optimum tra- 
jectory, If at the initial moment of time the tracking 
triangle occurs, then ®.=0, i.e., the tracking line does 


not change its direction in space and = =const 
t 


({ ®, dt =0): Otherwise, we have &. #0 and the 

¢(1) 

tracking line turns the angle | 2) —@2) | (this is a mono- 
tonic change in the angle ®., as can be seen from the 


first equation (1.1)]. The less |@2-—#™ | becomes, the 
nearer is the transition path of point A to the path &.=0, 
i.e., to the tracking triangle whose occurrence ensures 
(within the limits of the proposed problem) a mathemat- 
ically precise interception. It is proposed? that the 
equation bringing point A from its initial state to the 
tracking triangle be called an optimum one if L= 
(2) 
= | \ @,dt| attains its minimum value; min L = 
(1) 
= min e290) , this min jo2-¢| being equiva- 


lent to min &” (& >0) or max 9 (&,<0), 

The trajectory corresponding to a controlling opti- 
mum (in the above meaning) will be called an optimum 
path. 
Note. a? = or ), where (is the moment of 
time either of the point A entering the tracking trian- 
gle, or of intercepting the goal B, 

We introduce the notation ® =x,, ¥ = Xz, pip X= Xz, 
p» =pce const, A parametric representation of &t) 
is introduced as done by Miele [2, 3}. We represent &t) 
as a smooth function of a new controlling parameter p 
satisfying the following conditions; 


¢<=0 when —o< Pp pi, 
Ex= 4 when PaSp<+, 
— = E(p)when Pi SPS Pz 


(a single-valued transformation), 


* Subsequently, we shall understand by “point A® the 
center of gravity of the controlling body A. The equa- 
tions of motion of the body about its center of gravity 
(about the point A) are not considered, 

Tt The goal B is understood to be a body. 

+ The criterion of optimum control was introduced by 
B, V. Shirokorad. 
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Fig. 1. 


Due to this representation, the role of the control 
function is taken over by p. Owing to the way the func- 
tion is introduced, one is able to apply the methods of 
classical calculus of variations to find the optimum so- 
lution. 

It should be noted that the Miele method does not 
pretend to be unique. In the last few years Soviet and 
foreign research workers [4-8] have submitied quite ef- 
fective methods of finding optimum processes in systems 
with restrictions, The Miele method, however, is espe- 
cially suitable, as it permits reduction of the problem 
with restrictions to Meyer's problem, and thus allows ap- 
plication of the full apparatus of the classical calculus 
of variations. 

The simultaneous equations (1.1) are now written 
in the form 


: i : : 
4= zp) [ve sin (pe — 2%) —v sin (%—%)], 


E(p) 


t, = 
xp? 


[ve sin (~_e — 2%) —vsin(%,—2,)}, (1.2) 


r= Sur [v_ Cos (Pe — X71) — V COS (X_ — 2)], 
where € = E[p(t)]. 
The boundary gonaitions are: 
At t=t'".xj", X2°, Xs are given; 
At t=t" we have 


v_ Sin (pe — 2\”) — vsin (2 — 2) = 


a” <0 


0, (1.3) 


a 


is not known). 
Meyer's variational problem [9] can be stated in 
the following way. 

Within the class of functions x(t), x{t), Xs(t), p(t), 
satisfying the equations of motion (1,2), and the boun- 
dary conditions (1,3), a system of functions is to be de- 
termined such that it minimizes |x; )- ad 

The solution is sought within the class of sectional- 
ly smooth functions x,, Xz, and x3; the function p(t) may 
have discontinuities of the first kind. We shall call the 
path as given by x, Xg, Xs and p admissible when Xj, 

Xg, X3 and p satisfy all the conditions of the Meyer's 
problem. A control program corresponding to an ad- 
missible path shall be called admissible. We shall con- 
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sider one particular case prior to solving the general 
case. 


2. A Particular Case 





We assume that | q%—x,| and |x,—x,| are small 
(for example, not exceeding 10°); this is often encoun- 
tered in practice. Now the simultaneous equations (1,2) 
and the boundary conditions (1.3) can be expressed in 
the following way: 


Zy = az, —br,+¢, Ie = (ar, — br, + c)&(p), 


vn U o--F, 














Is = Z3 = 4 ———_—-I, (2,1) 
op)’ pit) 
Asis, v Ye Pe 
ae SO 8a 2, Ce Sb 
arr Pant zap) a, zap) “ 
or 
v—v, a v 
fe )_(v—v,)t’ a p) —(v—v,)t’ 
Ve (2.2) 
“7 M_ (w—v)t 
The boundary conditions are; 
At t=t% Xy, xy are given; 
At t=t2 


~ le) (2) - 
az, — bry’ +c= 
1 Ig c=(@ (2,3) 


(a =azp™, b = bry, ¢ = crap). 


We introduce the notation ® = ax,—bx_+c. The 
Meyer's variational problem, as applied to this particu- 
lar case of the simultaneous equations (2.1), can be de- 
scribed in the following way. In the class of functions 
X,(t), Xt), p(t) satisfying the equations of motion (2.2) 
together with the boundary conditions 5 -3) 4, tem of 
functions is to be found minimizing [x -x¢ |. 

In order to solve this, we introduce Lagrange's func- 
tion: 


F =}, lax, — bx, + ¢ — 24) + 
+ hg [E(p) (ax, — bx, + c) — 23). (2.4) 
Euler's equations take the form 
ah, + a (p)%, + 1i=0, 
ah, + bE(p) dy — 2, = 0, (2.5) 
Xe (ax, — bay + 0) % = 0. 


From the last Euler equation of the (2.5) system 
it can be seen that the optimum path (when it does ex- 
ist) may consist of three branches; a) dt/dp = (A, #0), 
i.e.,€ = 0 and — = 1 (two branches); b) A, = (dE /dp # 0), 


i.e., the controlling function — is not constant(ax,— bx,+c #0 


as the tracking triangle does not occur). 








), 
é 0), 
+c #0 





We shall deal first with the case of 4, = 0(dEAp#9. 


It follows from the equations (2,5) that when A,= 0 then 
also A4y=0, But both the Lagrange multipliers A, and A, 
cannot vanish simultaneously [9] at any point of the in- 
terval A), t ) as it would eventually be inconsistent 
with the boundary conditions (2.3) being independent. 
This means that a curve with a variable €,and such that 
A,= 0, does not exist. Therefore, a possible optimum 
trajectory consists of two branches: €=0, & 1dEAMp=9. 
We shall now show that the branch € = 0 does not 
satisfy the boundary condition at the upper end. As as- 


sumed, the initial point is in the region of ® #0 (the 
tracking triangle not increasing). 

It should be @ “= 0 for an admissible path at 
r=t”, Consequently, if, say,@ )> 0 then © must, in the 


course of time, decrease (not necessarily monotonically). 


How this function changes in time during motion along 
a trajectory of the simultaneous equations (2,1) can be 
found by taking the derivative of @ with respect to 
time; due to the equations of motion (2,1), the latter is 
obtained as 


q® _ = _ Fp, — (be (ar—ba + _ a—Be @ 


xp) zap 





When £=0 we have ? _ 2 @ As %oa=v-,), 


= 2 
x3p)> 0, the sign of d®/dt(—=0) agrees with the sign of 
®, Therefore, if ® ’>0 in the case of motion along the 
path with the control €=0, the function ® must in- 
crease monotonically, but if ® Yeo it decreases mono- 
tonically. This shows that the branch €=0 cannot satis- 
fy the boundary condition at the upper end, 


mg « v 
When £=1 we have d® = s—S@ ies ae ee 
dt aap) 


zap) 


Consequently, d®/dt(€=1) is opposite in sign to ®, i.e., 


the branch € =1 could be a solution of the boundary prob- 


lem, 

We shall show that the branch €=1 is actually an 
optimum one, i.e., gives the min | xi"—xy"| . In the 
system of equations (2.1) we shall take x, as the inde- 
pendent variable (we shall show later that this is per- 





\, 





: 
. 
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Fig. 2. 


missible). Then the simultaneous equations (2,1) reduce 
to a single equation: 


dz 4 
fe Fay? OSMa)<4 (2.6) 


We shall determine the admissible displacements 
in the (x2, x4) plane as those which are not consistent 
with (2.1) or (2.6) (10). Let I represent point corres- 
ponding to the initial moment of time t (Fig. 2). We 
shall now determine where the point gras be displaced 
during the subsequent time interval t" + Ot, 

It follows from (2.6) that 1 S(dx,/dx,)S+@ . The rep- 
resentative point can therefore at the time t de be 
found only in the region bounded by the straight lines 
(dx/dx,) =+ and (dx,/dx,)=1, i.e., in the zones A’ or A 
(Fig. 2). B is clearly the prohibited zone, 

We draw the straight line ax} bx) +2 =0 in the 
plane (X2, x4), on which the representative point at the 
time moment Pig. 3) should lie. The gradient of the 
E > 1, as v>v%. There are 


v-—v, 





straight line is 3 = 
a 


three possible positions of the initial point I with respect 
to the straight line @**=0; 

1) the initial point is below the straight line #”=0, 
therefore, o%< 0. 

2) the initial point is above the straight line 6X0, 
therefore, ®**>0, 

3) the initial point is on the straight line #0, 
i.€., ©) also and thus, the tracking triangle condi- 
tions take place straightaway, with no transition process, 
This case of course, need not be considered, 

We now return to the first two cases, 

1. The point ly is below the straight line #=0, 
i.e., in the region of ®<0, It follows from equations 
(2.1): 


—,, &®. (2,7) 











/ Fig. 3, 
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As xsp"20, and £20, then below the line #”)=0, 
X, and xX, are monotonically decreasing functions of 
time, This makes it possible to take x, instead of time 
t as an independent variable, and also turns the zone A‘; 
into a prohibited one, as displacements consistent with 
the equations of motion (2.1) and (2.6) are found in the 
zone A; (Fig. 3). 

This confirms the existence of an optimum trajec- 
tory in the case under consideration. It will consist of a 
straight line 1 (€=1) until the latter intersects the 
straight line @*”=0 at the point M;. The intersection 
occurs because the gradient of the boundary line @@)9 
exceeds unity. The optimum character of this trajec- 
tory is obvious (Fig. 3). From the point I, (Fig. 3) we 
draw a straight line parallel to the boundary line @*’=0, 
Its gradient is vv— Vc). It was assumed that 0< y < v; 
hence 1<vAv—Ve)< +, If the control €= const takes 
place (the path consists of one branch) it is obvious that 
the admissible controls must satisfy the conditions 

v—v,. 


<t<i, — = const. (2.8) 





v 


Thus, the control € =0 (if this one occurs only) 
proves not admissible. The corresponding path (the 
branch €=0) must, therefore, be discarded, This is in 
agreement with the result obtained analytically, In ad- 
dition, €=0 signifies that the body A is not controlled, 
Disregarding the exceptional case when it is directed 
toward the goal already at the initial moment of time 


(the tracking triangle has occurred), we can draw the 
conclusion that an uncontrolled body A cannot intercept 
the goal B. 


2. The initial point I, is above the straight line 
oo, i.e., in the region of @>0, The equations (2-1) 
indicate that in this region x, and x2 increase monoton- 
ically with time; thus, the zone A} becomes prohibited, 
and the displacements consistent with the equations (2.6) 
and (2,1) are in the zone Aj. The existence of the op- 
timum path is confirmed, This path runs along the 
straight line I,M, (€=1) until it meets the straight line 
©®)-0 at the point Mz, The intersection can take place 
as the gradient of the straight line @“’=0 exceeds unity. 
This path is obviously an optimum one (Fig. 3). 

From the point I, (Fig. 3), we draw a line parallel 
to the bounding straight line eo, Its gradient is 
vA v— v,)and 


l< as O<v,<»v. 





< +0, 


A control €= const is obviously admissible provided 
it satisfies the condition (v—v,)/v < € =1(€= const). 
Therefore, the control €=0 is not admissible. The cor- 
responding path (the branch €=0) is not admissible ei- 
ther and should be discarded. Thisis in agreement with 
the result obtained analytically. 
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Note, The equations (2 1) are valid in the region 
|%—x,| 10°, | x»—x,] =10°, i.e., within and on the 
perimeter of the parallelogram A’B"E'D' shown in Fig. 4, 
The center of the parallelogram lies on the bounding 
straight line ©) (x = and x= G satisfy the equa- 
tion @*’=0), It is obvious that the solution of the varia- 
tional problem is possible when ths initial points lie 
either in the region ABD (when ®@**>0) or in DBE (when 
<0) (ABIIA'B’, DEND'E’, ZBAD=n/4), 

It would be interesting to find the time of entry 
on the tracking triangle under the optimum control €=1. 
When €=1, the equations (2.1) can be integrated. We 
write the optimum curve equations in the form: 


Ve 





v—U_ 


1 = Mi (—y) + Ny, 


Ve 
D—_ 


p 
m= M(t) + Nz, (2.9) 
M, = My = +{(ve —v) 2 + val — .v:}, 
Ny = [0 (2 — 2) + evel, 
c 
Nz = = [(v atin Ve) (x) sian a{?)) + Pec]. (2.10) 
As the optimum path satisfies the boundary condi- 


tions, one is able to find t 2 the time of entry on the 
straight line #0, We have 











ov 
— 
Pe * pi?) = o =~ 
(—5) Ms; (ay) + aN, — 6N,+c=0 
" § 
ss 
SAY 9/8 
4 
4 
gid iy € f° 
oF 
go 
0 %. ~ 
Fig. 4. 








or, taking into account the expressions for 2%, 6 and 7, 
Ve 

(2) v—Ve 

— vel (ay) + (v —v,) Ny —vN, + ore = 0. 


(2,11) 


It is easy to show that (v-\y,)Ny—vN2+@V, =0 for 
any xj’ and x from the region where the equations (2.1) 
are valid. Then it follows from (2.11) that 


U 





(2), °%e 


m() = 


We shall distinguish two cases; 

a) M,=0, p%#0, But My=—(2%%,); therefore, the 
vanishing of M, denotes that already at the initial mo- 
ment of time t*”, the tracking triangle occurs,which is 
contrary to our assumption. 

b) M, #0, p= 0. Then t®)is equal to the time 
‘> of intercepting goal R; the latter time can be obtained 


from the equation p = oo) —(v—v,)t for p=0 as 


ty = PY - .); therefore, 1%). tp =P MW): 


3. General Case 





We now proceed to the equations (1.2) and the gen- 
eral variational problem. There are the same stages in 
solving it as in the case considered previously, The La- 
grange function F becomes 


F=i,{— + =a [v, sin (p, — 2,;)— v sin(z,— )I\+ 


Ey, sin (p, — 2,)— vsin(z,— z,)I}+ 
sp 
+ ds {— 23 + —ar le COS (%¢ — %)— v Cos(1,— z)}} 


+{-% + 


(3.1) 
The Euier equations are 


3 [v cos (2, — 23) — vcos (Pc — Z;)] + 


4. ath) [v cos (%_—2}) — U,.COS (%e—2;)] + 


+ dg [ve Sin (p, — 21) — v sin (%_ — 2,)] + 4yo™ = 0, 
E (3,2) 
v 
— Ay {008 (%2 — Xi) —1, =) cos (%_ — 2) + 


+ Agu sin (1, — 2) + Ap = 0, 


4 [v sin (x, — 2.) — v, Sin (P_ — 2)] — 
3 


— 1,22 |. sin (ge — %1) — 


— vsin (1, —2,)] + Ap) = 0, 


d : . 
4s & {v, sin (pe — 21) — v sin (zy —2,)] = 0. 


The last of the Euler equations (3.2) shows that 


the theoretically possible optimum trajectory may con- 
sist of three branches: 


a) (d&Ap)=0, (Ag #0), E=0, E=1. 


b) A,=0, d&/dp #0, i.e., € is not constant 
Lv, sin (%—x4)—v sin (x2—x)) # 0, as the tracking tri- 
angle does not occur in accordance with our assump- 
tions ]. The branch with variable € satisfies the follow - 
ing simultaneous equations obtained from (3. 2) by put- 
ting A2= 0. 


a [v cos (x7, — 2) — v, cos (p_ — 2,)] + 


+ hg [ve Sim (Pe — 21) — 
— vsin(z,— %)} + Ap = 0, 
—h, ~ cos (12 — 2) + Agv sin (7, — z,) = 0, 

73 
"i eo. em 
= [v sin (x7, — 23) — v, sin (~, — 24)] 4- Ago™ = O. 
3 

The function F of (3,1) obviously does not depend 

on time; thus, Euler"s equations (3.2) have the following 
first integral: 


RiZy + Age + Ags = C = const. 


To determine the constant C of the first integral, it 
is necessary to consider the transversality conditions [9] 
which easily provide that C =0 (the function 1) 
to be minimized does not depend on t””, i.e., the varia- 
tional problem is time-independent). Making use of 
equations of motion (1.2) we write the first integral in 
the form: 


ms [v_ sin (~- — x) — v sin (xz, — z,)} + 


(3.4) 
+ ef (p) [ve Sin (pe — 71) 


— vsin (1,— 2))+ 
+ hg [v_ COS (Ge — 1) — V COS (7, —z,)] = 0. 
When A,= 0 we have 


M [ve sin (p. — 2) — vsin (x, — 2,)] + 


+- hg lv, cos (p, — 21) — veos(z,— 2,)} = 0. (3,5) 
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From the second equation in (3,3) and from (3.5), 
the relation between x, and x, is found which must take 
place when € is not constant: 


t, COS (@,, — 21) —v COS (rg — 24) 
v, sin (9, — 2) — v sin (z3 — 7) * 
(3.6) 


Here vy, sin (%.—x,)—Vv sin (x2—x,) #0, as the track- 
ing triangle does not occur; in addition, Ay and Ag cannot 
vanish simultaneously on the curve with variable €, be- 
cause when A,=0, it follows from the equations (3.3) 
that also Ag=0, But the Lagrange multipliers cannot 
vanish simultaneously [9] (on the branch with variable € 
we have \,=0) at any point of the interval (t"’, t'’). 

From the equations (3.6) we obtain the relationship: 





bg (%—-%4) = — 


cos (Pc _— 1s) = -. > 4 (3,7) 
c 


The inconsistent relation (3.7) was obtained as a 
result of the assumption that there is a curve with vari- 
able —€(A,=0). There are no branches, therefore, when 
the controlling function is variable, The theoretically 
possible optimum control is € =0 and € =1, i.e., the 
same as in the particular case above. 

Further, it is necessary to clarify as to whether the 
controls € = 0 and € =1 are admissible, i.e., whether 
the corresponding paths satisty the boundary conditions. 
It is obvious that both paths (€ = 0 and & = 1) pass through 
the initial point, We introduce the notation: 


rv, sin (p, — %) —- vsin(7, — 27,) =G, 


(3.7') 
Gg at G (t), g™ = G (1) 


(according to our assumption oz 0). An admissible 
trajectory should have 


G™—0, z%<0, (3.8) 


at t=1%, and if, say, o> 0, then the function G will 
finally decrease along the path (but not necessarily 
monotonically), To study how the function G varies 
with time, its derivative with respect to time is con- 
sidered, and in view of the equations of motion (1.2) 
we obtain 


dG d ‘ . 
a = a e8in (Ge — 2) — vsin(z, — 2,)|= 


= —- Ue COS (Pc— 7) 2, = 
~~ ¥ C08 (%_—2,) (2 — .1;). 
But x= 0 when € = 0; thus 


(Fane == — (24290), = — ( Gr; his _ (3.9) 
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and Xy= Xz when € =1; hence 


dG 2 
(ar). -. [11% COS (P_ — Xl, “ay 
(3.9°) 
Bie cos (@, — 21) | ; 
7m b=} 


zp) 
(G)z—, ((G):.-,) denotes that the function G is 


taken on the path with control € = 0 (€ =1). 

We consider the condition (3.9) when — =0. It in- 
dicates that the path with the control € = 0 can satisfy 
the first boundary condition in (3.8) (G® = 0) in the re- 
gion x,> 0 only. Therefore, the second boundary con- 
dition in (3.8) (x 0) cannot be satisfied. The branch 
— = 0 thus proves inadmissible and must be discarded, 
The corresponding physical meaning is the trajectory of 
an uncontrolled body A, It is obvious that if the body A 
is not directed at t= t®) toward the meeting point (the 
tracking triangle does not occur) then the further motion 
along the path with € = 0 cannot take it to the goal B. 

The condition (3.9°) when € =1 shows that the path 
with the control € =1 can intersect G*’=0 at t=t only 
in the region where cos (@.—x,) > 0. The significance 
of this condition will be explained later, But it can al- 
ready be said now that it is not a very severe one and 
that it can be satisfied, 

We shall show that the path € =1 is actually an op- 
timum one, We consider all the admissible displace - 
ments in the plane (x2, x,) [i.e, the displacements which 
are not inconsistent with the equation of motion (1.2)]. 
This was certainly possible in the simpler cases of the 
simultaneous equations (2,1) or (2-8). 

In the general case, the possibility of such a con- 
sideration is determined by the following facts, First, 
by studying the Euler equations (3.2) it was found that 
such a curve does not exist when the control function € 
is not constant (d&/dp #0). In an unfavorable case when 
A2,= 0 one could obtain a relation between x4, X, and Xz 
which would be fulfilled on the curve with a variable 
control function €. Then it would be necessary to con- 
sider all the admissible displacements in the space of 
the variables x;, Xz, X3, and this would make it rather 
involved, 

Secondly. the first boundary condition at the upper 
end Gj" =0 is independent of xg". The function xs de- 
pends only on x, and x, and the other boundary condi- 
tion xs <0 can, therefore, be investigated in the plane 
of (Xg, xy) only. In the initial condition, a vital role is 
only played by the initial values of the functions x, and 
X2 (this is indicated by the physical aspect of the prob- 
lem). All this enables one to consider admissible dis- 
placements in the plane (x2, x) which are given by the 
first two equations in the system (1.2), the equation 
dxy/dx, = /E (x2) on condition that 0 S€ (x,)=1. It is 
reasonable to refer our discussion to the square of side 2, 
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its center coinciding with the origin( Fig. 5); x, and x, be~ 
ing angles, one can consider x, + 2km as coinciding with 
Xy, and X2+ 2mm with x, (k and m are integers), The 
boundary branch G“’ = 0 (3,7") has three branches in the 
square: GY’ =0, GY) =0, GR) =0 (Fig. 5). The branches 
are in fact arcs of some periodic curve with period 
equal to 242m along the appropriate axis. The straight 
line xy= x, which is cut by Gf?) = 0 at the points(¢,., ¢) and 
(-1/ 2-4, — / 2— ) serves as the axis of G\”) = 0. 
The width of the belt (along axis Ox,)containing 

c?) =0 is 2arc sin(v,./v).The other two branches cP) =0 
and GS" =0 “meander® around the straight lines Ow Xot+ 7 
and Xy=X,g—-m, We denote the region between G7” = 0 


and c? =0 by Ay, and the region between c? = 0 and 


co? =0 by Ay, the region above c= 0 by A,, and the 
region below cy =0 by 4. 

The function G is positive in the regions 4y, and 
4,, and negative in 443 and 4), In addition, the curves 
X3=0 are also shown in Fig, 5, i.e., the curves 
Ve COS (%—X4)—V COs (X_2—X,)=0. We introduce the 
notation v, cos (G%—x4)—V cos(X,—x) =), 4 The curve 
g =0 has four branches in the square of side 27; 


81=9 g2=0, gs=0, gy=0. 


Notation; 4y, is the region between g,=0 and 
Z2= 0, dy the region between g,=0 and g,=0, dy3 the 
region between g,=0 and gs=0, 53 the region above 
g3= 0, 5g the region below gg=0. We notice that %;> 0 
in the regions 5y3 and 5, and &3< 0 in the regions 49, 
55, and 54. 


The tangents to gy and g» are drawn at the point 
X1=% (their gradients are unity), The strip between 
the two tangents parallel to the straight lines x; =x, has 
the width equal to 2 arc cos (v,/v) (if measured along 
Ox,). When v, <( 42/2) v (usually there is v, = 4v) we 
have in the I part arc sin (v./v)< are cos(v,/v); there- 
fore, the strip embracing the branch € =1 lies entirely 
within the region 542(%3< 0). An analogous approach 
will show, on the contrary, that the branches G2“ =0 
and cP) =0 are er in the region of X3> 0. Thus 
the paths which at t=t ) are located on those branches 
are not admissible as obviously xf?) is positive, i.e., the 
second boundary condition (3.8) (x3 < 0) is not satis- 
fied, Thus the solutions of the boundary problem can 
ony be represented by the paths crossing the curve 
G}?=0 at some t=t* , the entire one lying in the re- 
gion 542(%3< 0). 

We denote the strip containing the curve c? ) 0 
by 4%) , the strip containing the curve fee =0 by 


42) and the strip containing the curve cf ) by 4.2) . 
2 3 


We now determine the admissible displacements in 
the plane (x2, x3) [10]. The admissible displacements 
are located in the regions 44, and A; —both within and 
on the boundary of the angle formed by the rays x; =1 
and x, = © (the angle opening upward ‘as in these re- 
gions on the function G> 0); in the regions 44, and 4,, 
inside and on the boundary of the angle formed by the 
rays xj=1 and x; = (the angle now opening downward 
as in these regions the function G< 0), 
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Further, we shall consider various positions of the 
initial point I in the square with side 2m (Fig. 5). The 
variational problem admits no solution for all the points 
I not in 4,0. as no path starcing from such an initial 


point can be admissible: none satisfies all the boundary 
conditions (3,8). Now let 1€4 cf All the admissible 


trajectories starting from I must stay in 4,0) , other- 


wise they would cease to be admissible. First it is ob- 
vious that the trajectory € =0 is not admissible for any 
be I, as it stops neither on the curve = 0 nor on 
= 0 in the region of %3> 0, i.e., the second boundary 
pe of (3.8) is not satisfied; secondly, if 1¢4 ys 


the branch € =1 is admissible and gives min (max) xy ) 

When the point I approaches the boundary of the region 

4 > the region containing the admissible paths shrinks 
1 


and in the limit (I is on the boundary of 4 ) it degen- 
1 


erates into the line € =1. 

Notes, a) Straight lines x,= %, +(1/2) and 
X1=—@,~(m/2) are drawn in Fig. 5, In the region be- 
tween these straight lines cos (¥.—x,)> 0. It is demon- 
strated in Fig. 5 that in the region of cos (%.—x,)< 0 the 
path € =1 cannot cross Gj’ =0. Such an intersection is 
only possible in the region of cos (@%.—x,)> 0. 


b) The tangent to the curve ,G c? = 0 at the point 
(%» %) has a gradient equal to x} = “ twiv- v,)] and also 
xj =0, i.e,, the point (%, %) is a point of inflection. 
Coming ee to the particular case in section 2 we no- 
tice that ®*”=0 has the gradient [v/(v— vj)} and that it 
passes through the point (@, 4%). The equations (2,1) 
are thus valid exactly in neighborhood of the point 
(%» @) where the curve y=0 has a point of inflec- 
tion. 

c) The general variational problem, as stated in 
section 1, can also be formulated (with the same opti- 
mum requirement) for a yeas pount) condition at 
the upper end, namely: p? )=0 when t= e The opti- 
mum path will now consist of two branches; the branch 
where TNS control € = 1 over, and the branch 
x, =x} ) = const, x a= xs =const, where xfFR) and 
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xfTR ) satisfy the trackin Coane forme i, ety: 


conditions v, sin(%, —X4 )y- v sin (xft 
&3< 0. The optimam character of this wen ‘a “x 
vious as the control € =1 gives min (max) xy”, and x, 
does not vary when the tracking triangle condition is 
satisfied, In addition, both branches can also be obtain- 
ed from the Euler equations (3.2). 

In conclusion the author wishes to express his thanks 
to B, V. Shirokorad and A, M, Letov for formulating the 
problem and constant advice given during the work. 
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ON THE STABILITY OF THE PERIODIC MOTION 
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The stability of periodic motions of piecewise linear automatic control systems is investigated on the basis 


of the resulting equation L(D)z = -K(D) f(z) + F(t). 


1, Considered in [1] is a dynamic system described 
by the equations 


n 


Zi= Dy Oy Zj+ Mf (21) + Fit) (= 4,2,..., 2), 
j=1 
(1) 


where aj, A; are given numbers ( some of them are actual- 
ly zero), F(t) is a periodic function of time having the 
total period T (in particular a constant number, which 
may be even zero), f(z)is a given piecewise linear 
function. The same reference gives a method for con- 
structing the linear approximation equations permitting 
solution of the problem of the stability of a periodic sol- 
ution %(t) of the given system. 

In the determination of the periodic solutions of the 
system (1) it is customary to start from the resulting 
equation [2,3,4] 


L(D) 2z,= — K (D)f (%) + F (), (2) 


where F(t) is a periodic function of time having the pe- 
tiod T, and 


L(D) = a,D"+ a,D"~'+ ...+ an, 
K(D) = boD®+ b,D—14...+6, © 


are linear differential operators (D is the symbol for the 
generalized derivative with respect to time), Usually, 
the degree m of the operator polynomial K(D) is lower 
than the degree n of the operator polynomial L(D), and, 
consequently, some of the first coefficients of K(D) in 
(3), ie., bg, bi, ...+ bp—m—y are equal to zero. 
Equation (2) is obtained immediately from system 
(1) after eliminating all the coordinates z; and their de- 
tivatives, except z,, Consequently, the coefficients ao, 
ayy. . +s Ans Dnems Dp—msy+---+ bp Of the polynomi- 
als (3) depend only upon ajj and A;. It is even easier to 


write equation (2) if the transfer function of the linear 
part of the control system is known: 


K 
W (p) = ree 


The polynomials (3) are here the denominator and the 
numerator of the transfer function of the linear part, re- 
spectively, 

The periodic solutions of equation (3) are obtained 
in the form %(t) or F(Z,(t)}. In order to solve the prob- 
lem of the stability of the periodic solutions Z, , f(Z), 
according to[1], one must use the reverse transformation, 
namely, express all the coordinates Z through Z,, and 
then introduce the linear approximation equation. 

Here we suggest another method for constructing 
the linear approximation and a direct method for solving 
the problem of the stability of the periodic solution Zp 
starting from the resulting equation (2). 

2. Therefore, let us consider an automatic control 
system described by the equation 


L(D)z=—K(D)f(2)+F(). 4) 


Let Z be a known periodic solution of equation (4), 
which can be ascribed either to an external periodic 
force F(t), or to the dynamic properties of the nonlinear 
system itself for F(t) =const. 

Let us denote by 7%, 7, .. » 7) the times of pass- 
ing through the links of the piecewise linear character- 
istic F(Z) for a period of motion T, where 7 5) ts 


j= 


let us also denote by t; the instants. when the coordinates 
%, JC) and their derivatives may be discontinuous with- 
in a period, i.e., t are determined by the formulas 


j 
=u G=1,2,...,0, 
i=] 











Se 





where {, = T. The subscripts of T and t are always in 
correspondence with the order in which they are met, 
both in one and in the other direction during the period 
of motion, and not with the geometrical sequence of the 
joints, 

If x(t) is a sufficiently small deviation of the per- 
turbed motion z(t) of the system under study from the 
known (unperturbed) motion Z(t), i.e., 


x(t) = 2(t)—z(t), 


this deviation is determined by the equation 
L(D) x(t) = — K (D) {x(t)[D.f(2)|,_+}, (5) 


where [D,/(z)],+ is the generalized derivative of 
the function f(z) with respect to its argument z for z= 2. 


If one takes 
} (2) = 4j2z + Yi» 


one finds 
[D./ (2)), + = 


‘et l 
=> > {ete —w —4_.)— 
siaaachar (6) 
Af: 
— 1(¢—wT —1;)) + aa 8(t—vT —1)}, 
Z P 
i 
where 1(t) is the unitary function of Heaviside, 5(t) is 
the Dirac momentum function, and Af. denotes a dis- 
continuity of the characteristic [ f(z) 13-3 at the time 
t= vT+t; » ABes 





Aj; = f|2(vT + t3+0)]— 


—/[z(T +t3—0)] = 
we be (7) 
=f [2 (t; + 0) —/ [2 (45—0)] = 


= (hits —hj)z(t) + Yin — te 


Formulas (6) and (7) are effective only when Z(t) 
is a continuous function of time. As is well known, this 
condition is fulfilled form =n~1. 

As to the derivative Z(t), it has discontinuities of 


the first order in the points of discontinuity of [ f(z)],-%. 


Consequently, for these points one must replace tj in 
formula (6) either by t;—0 or by t;+0. In the following 
we shall always yee t, +0, i.e., we shall consider the 
right-hand value z(t;+0). 

The product x(t{D,F(z)],.% will be written in 
the form 


co «(Cot 
=(0)1Dsf (2)),z = >, 5) [Ast (¢) + 


v=0 j=1 
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A 
+ et zr + t; +0) 8(t— wT — t)| : 
z(t; + 0) 
(8) 
where 
yj (t) = z(t) 


Zyj (t) =Q for 


for v7 + t).,§<t< vw +4;, 

Wi<t<vwP +t. (9) 
and vl + t;<t<vwW+T. 

If (8) is introduced into (5), we get 


o l 
L(D)2(t)=—K(D){3} 5} [Ayes (¢) + 


v=xQ j=] (10) 
Af; 
+2—— 2 (oP + 5+ 0)8(¢—v—1) | 
z(t; + 0) 


Let us show that equation (10) is really a linear ap- 
proximation in the sense of Lyapunov, and that, there- 
fore, it can serve for deciding whether the periodic sol- 
ution Z(t) is stable or not. 

It was shown in [5,6] that the linear approximation 
for solving the problem of the stability of a periodic sol- 
ution of the system (1) has the form 


Du = > Qijtj;+ [Dif (2)) 7% (i = 1, 2, .- +2 2) 
j=1 (11) 


As has been done on system (1), where, by elimin- 
ating all the z;"s except z,, and by replacing z, by z, 
we got the resulting equation (2), system (11) can give 
the linear approximation (10) corresponding to the re- 
sulting equation (2). 

In fact, if Xg, X3, ..., xX, are eliminated from (11) 
and x, is denoted by x, and if the expansion (8) is taken 
into account, the actual linear approximation of (11) 
takes just the form (10), After all, equation (10) is noth- 
ing but another way of writing the linear approxima- 
tion (11). Consequently, equation (10) solves the prdéb- 
lem of the stability of periodic solution Z(t) of the re- 
sulting equation (2), similarly to a well-known theorem 
by Lyapunov: if the zero solution x(t)=0 of the linear 
approximation (10) is asymptotically stable, then also 
the periodic solution Z(t) is asymptotically stable. 

3. On the basis of (9), for each period of motion, 
ie., for vy T< t=vT+T, (10) gives the equations 


[L (p) + kjK (p)] x(t) = 0, 

d (12) 
P= (f=1,2, ...0, 
which describe the perturbate process in the given in- 
terval of time, except for the discrete points vT+t.. In 
these points the deviation x(t) and its derivatives lose 
the discontinuities ax) » which are determined by the 
system of equations (see the Appendix): 


0= —Af; -—*— 2 0T +t; +0), 


z(t; + 0) 











(@o + kybo) Axyj = — Akjbox (vT + t; + 0) - 


-Af;x 2(vT + t;+0), 
z(t; +0) 


(do + Ajbo) Ax; + (a, + hjby) Ax; = (13) 
= — Akj[box’ (vT +t; +0) + b,2(vT + t; +0)} — Af; “ x(vT +t; + 0), 
Z 








(ao + kjbo) ax” ++ +++ (Gn—kjbn_) Az; = 
= —Ak;[box'™—» (vT+ tj-+ 0) +--+ by_yr (vl + t; + 0)) — 


—Af;z——_ 2(vT +-t, +0), 





where we have used the notation Akj = kj+,—kj. 
The first equation of (13) means that the coefficient by must be equal to zero, as has already been men- 


tioned, The other n equations ¢ of (23) are a system of linear equations which can be solved with respect to the 
jumps Axy j, Axyj, a6 , ag j ). 


Let us write the corresponding solution in matrix form 



































AX; = A;X (vT f t; -t- 0), (14) 
where 
Ax; x(t) 
Az,; x’ (t) 
AX; = : » X(t) =] ° ’ 
Ax(y—0) ai™—1) (t) 
0 0 0 0 0 0 
aii + Bir 0 0 0 0 0 
A;=|| ®t+Ba a OO 0 0 Oo] , (15) 
ais + Bis Cy, 86 0 0 0 
Ajn—1 + B;, n—1 @j,n—2 %j,n—g Zina. - - Bin 0 
The elements of Cine aT of the matrix (15) are x(t) is reduced to a storing of the solutions of (12) taking 


into account the jumps (14) in the interval of a period. 
4. After integration of the linear approximation 
(12), (14) in a period of motion, i.e., for vT<tSvT+T, 


obtained from the recurrent formulas 











fa * Ak, s a4 + kb, _4 3 we get 
To” ee Pei yt 
ge a dg X(T + T) =UX (v7), (16) 
bass Af; — a, 4+ k by, where U is a square matrix. 
Bin aaa “— ot — >} a Bits Let us determine the matrix U. 
t (t5) =1 Let x;(t) be a solution of the equation 
where [L (p) + kjK (p)) z(t) = 90 
Af. Let us write this solution in the form 
bg j 
a, = — Ahi Pa — 
“ Z(t) n 
x;(t) = Sy ei" Cy, (17) 
i=] 


Equations (12) and (14), taken together, are equiva- 
lent to equation (10), and represent another form of the 
linear approximation when x(t) is varied. On the basis 
of these equations, the determination of the variation of L (dj) + kjK (Ay) = 0. 


where Aj; are the roots of the characteristic equation 
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ee eee 
mn 











(case of simple roots) and the C;;'s are integration constants. 


From (17) we can build the matrix 








XOT +t, +) =X, = M; HC, (18) 
where 
ent eriat erin Ci 
M;(t) = hues hne viet Mine hal Cj= es . 
Mrte*at arterial Rot dint Cin 








The integration constants are determined by the 
equation 


C; = Mj*(0)X (wT + 3.40). (ag) 


On the other hand, if we store the solution (18) for 
two successive intervals j-1 and j, we get 


or, taking (14) into account 


X (wD + t3-1 + 0) = (E — Aj)? Mj (tj) Ci, 


(20) 
where E denotes the unit matrix. 
If (20) is introduced in (19), we get the recurrent 
formula for the integration constants 


Cj = Mj" (0) BAMja(tja)Csa, ©) 
where Bis is a nonsingular triangular matrix 


Bi = E aad Aj_,. 


Taking into account the fact that, according to 
(19), 


C, = M;*(0) X(T + 0), 


formula (21) can be written in the form 


C; = U;X (vT +0), (22) 
where ie 
U; = Mj" (0) Il j—nM jn (tj-n) Mj—n (0) 
n=1 
(23) 


(f=22, 20.404 Me 














The product of the matrix under the sign II is made 
to the right with n increasing. Here, My. + (0) can be 


replaced by Mj; 4(T). 
For j=2+1 and t =0,formulas (18) and (22) give 


X (vP + T + 0) = Miz (0) Ci4: = 


= Mi+; (0) U14.X (vT + 0), 


whence, according to (16), we get 


U = Mi4; (0) 014: (0) 


or, taking (23) into account, 


i 
os ll Bien M jn (Tj—n-+1) Mints (0). (24) 


n=l 


Formula (24) gives the matrix U of transformation (16). 

5. It is well known that the zero solution of equa- 
tion (10) is asymptotically stable, i.e., the variation 
x(t) and all its derivatives down to the (n—1)th deriva- 
tive tend to zero when the time increases beyond any 
limit if the characteristic numbers of the matrix U of 
transformation (16)( i.e., the roots of the characteristic 
equation of the linear approximation 


det (U — pE) = 0, (25) 


where E = the unit matrix,p= a scalar parameter) all lie 
within the unit circle in the plane of the roots, 

In this case, the periodic motion Z(t) is asymptot- 
ically stable. 

Reference [6] showed that in the autonomous dis- 
continuous case the periodic motion is stable also when 
one of the roots of the characteristic equation (25) is 
equal to unity, provided all the other roots are less than 
one in modulus. But if only one of the roots of the 
characteristic equation (25) has a modulus greater than 
unity, the corresponding periodic motion is unstable 
(this is a theorem analogous to that of Andronov and 
Vitt). 








APPENDIX 
Let us substitute in equation (10) 


co 


D%z(t) = p*z(t) + > >) [Az, 5°» (t —vT — ty) fone 4xig-Y) S(t 


v=0 j=1 


—vI —t,)], 


co lt co il 
p|>) > kx (| = hyp (t) + >) Dy (4 [eye () 82—Y (@ — we — 8) +--+ 


v=0 j=1 v=0 j=] 


+A [ke (QI 50¢—wF— th, (p=), 


Ax] = 


where 
A [kx (t)) = k,, 2 (vT + t+ 0) — ko) (VT +t;— 0) = 
= ko") (vT + t;+0)—k; jz (vT + t;,+ 0) — 

and = 2") (vT + t; +0) Aky + kjAzy? 


ml 


+ 0) 


Bata 


To 


a (vT + t;+0)8(t —wr—1)] = 


A 
=>} hath ete + 5 +0) 6@ (¢ vr —t)), 


vary jung 2 (8; + 0) 


Then let us equate the coefficients of p and gf), As a result of this operation, we get equation (12) and, 
correspondingly, for each of the discontinuity points t = vT + tj, we obtain the equations 


0=—Afjz 


agAz,y = — by [Akye (VT + t, +0) + ky] — Af, ——"X 


by 


F(t; +0)’ 


z(vT +t, +0), 
z(t; +0) 


aAz,;+ aAz,, = — by [Ak,a’ (vT + - + 0) + k;Ax,;] —b; [Ake (vT + t,+-0) + kyAz,.)— 


aAzxi?—Y +...4 4 


b,_, [Ak x (vT + t; + 0) + k,Az,,] — Af; ~ 


By combining the terms of these equations according to the discontinuities ax? 
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The reduced equivalent amplification factor of a nonlinear element in the presence of noise is found; it 
is suggested that it be used in studies of periodic modes in self- oscillating automatic control systems when high- 


frequency noise effects them. 


Let us examine an automatic control system com- 
posed of a linear and a nonlinear part, and described 
by equations 


a'x aly 
he te hie cy 
Y = fu (2), (2) 
Z(t) =X (t) + Xn), 9 


where (1) and (2) are the equations for the linear part of 
the system and the nonlinear inertialess part, respective- 
ly. 

A vibratory signal, or noise in the form of a period- 
ic disturbance X70, whose frequency spectrum lies be- 
yond the pass band of the linear part of the system, w., 
acts on the input of the nonlinear element; according to 
equation (3), the output function X(t) also acts on this 
element, Later, we shall examine the case when X(t) 
is a random function. 

We will be interested in the periodic oscillations of 
X(t), whose frequencies lie between w, > w > (w,/2) (Wy 


is the cutoff frequency of the linear part of the system). 
In this case, X(t) can be represented as a harmonic func- 
tion of the type 


X (t) = Acos ot, (4) 


where A and w are the amplitude and frequency, respec- 
tively, of self-oscillation X, The oscillations of X can 
be determined by the harmonic-balance method, and 
the presence of noise is reflected only in the form of an 
equation for the conditional equivalent amplification 
factor of the nonlinear element. This equivalent ampli- 
fication factor will depend on the type of nonlinear ele- 
ment, the amplitude and frequency of the unknown peri- 
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odic oscillations, the shape of the disturbance, and also 
on the phase shift between the disturbance and the un- 
known periodic oscillations. When the frequency of 

X7(t) is high, the effect of the phase shift between the 


functions A cos wt and X;,(t) and their frequency shift 


is not very great, and one can introduce the concept of 
an average (reduced) equivalent amplification factor. 

The function Z(t) may be either periodic or al- 
most periodic, but in either case, the first harmonic am- 
plitude of X(t) is nearly A at high frequency. 


We will assume that the frequency spectrum of dis- 
turbance XW (0 is large enough so that A can be taken 


as the first harmonic amplitude of the function Z(t). 
Let us approximate the function X(t) by a step 


function Xo (Fig. 1a). Then expression (3) for Z(t) 


takes the form 


Z(t) =Z(t)= Acoswt+ Xy(t). (5) 


The function X during period T can be repre- 


sented as the sum of n rectangular pulses, with width 
At and height (i=1, 2, 3, ...., Mm), equal to the or- 
dinates of the steps of function X,,(t) in the time inter- 


val At(i-1)<t<idt, The function eG) can be repre- 


sented as a series 


Xn (t) = 2 Xn, (0), (6) 


i=] 


for an infinite time interval; the series consists of im- 
pulses 








r- 





Xu. (t “ for (¢ —1) At + $t >t Didt + Br, Just as for the function § X7,(0), an infinite num- 
uy I]; for (i — 1) At + Br ct < idt + Br ber of values for 6 =1, 2, 3,. 
(7) each value of i. 
As can be seen from expression (9) and Fig. 1c, the 

Here i is the ordinal number of the impulses of (7), impulses of 6X;(t), in contrast to those of Xy,(0, are 
and of the | pulse of Xq (0, as shown in Fig. 1a and 1b; not rectangular. 
6 =1, 2,3, ..., © is the period number of X qn): In In case the ratio T/r = y is an integral number, 
the period of the functions 5X;() is equal to T, and 
each of the functions 5X,(t) contains y pulses during 
period T. If y is not an integral number, the period of 
the functions will equal cT, where c is the least inte- 
ger for which the product cy is also integral. 

The impulses of 5X;(t), like those of6 Xp, (0), 
(8) will follow one another at intervals Tt , and will be con- 


— current when indices i and 6 are identical (Figs. 1b, 
a lc, and 1d), In this case, 


+» @ is assumed for 


the function Z(t), A cos wt can also be represented as a 
series which consists of n periodic impulses 


6X,() (i= 1,2,3,....m) 


Acoswt = >) 8X; (t), 


i=1 


8X. = 
i (t) (9) 8Z; (t) = 8X, (t) + 8Xn,(t) (i = 12,3)... 0) (10) 
O for (i — 1) At + Bt St Didt + Br, 


= 
Acoswt for (i—1)At+Pr<t<idt+fp. 4 impulse functions defined by expressions; 








32 =| 0 for (i — 1) Ot + Bt >t Didt + Br, 
tw | Acoswt +Xq for (¢ — 1) At + Bt << t < iAt + Br, ee 
and the sum of all functions 6Z,(t) (i= 1, 2, 3,. n) _ series, at the output and input of the nonlinear element. 


gives Z(t), Le., In this case, the equivalent amplification factor is 


defined by the expression 


n 
it os T T 
= Ss wt t)= 8Z,(t). (12) 
Z(t) = Acosot + Xn(t) a (t) In = Fa\ Yecosutdt + j 7q\Yesinwtdt, (15) 
0 0 
Let us write expression (10) in the form 
where 
8Z; (t) = (A cos wt + TH) (¢—ayat+pe<tciat+p:, (13) 
where II; is a constant equal to the height of the ith Yu = fu(Z,) = (16) 
impulse of Xy(0 during the time being considered, 
which depends on i and 6 (Fig. 1). cs. 
Let us shift X(0 in expression (12) relative to it- = fu >} (A cos wt + Tlh)g+e—nat+erct<+n t+ 
i=) 


self by an amount equal to k At, where k=0,1,2,... 
(n- 1). 
Then expression (12) takes the form 


Since any phase shift of X_(0 with respect to A cos wt 


Z,(t) = Acoswt + Xn (t —kAt) = is equally probable, and has little effect on the 
(14) equivalent amplification factor if X n® is high fre- 
oe quency, we will define the mean po. Nan amplifica- 
™ Py (A cos wt + Ti) 42—rat+pe<tcei-phy attire. tion factor as the arithmetic mean 
We will define the reduced equivalent amplifica- ‘ n—1 17 
tion factor of the nonlinear element as the ratio of the hy = 4 > I, ay 


first harmonics of the functions, expressed as Fourier 





oe 


se es - 














for time shifts of Xy(0) by kAt as k changes from 


0 ton-1, Putting expressions (15) and (16) in formula 
(17) in place oft, we get 


n—1 


hy=2,\'t 2 fs 5! A cos ut + 


0 L fond 


+ Th) -pn—y atgpect<(i+n at+pe | C08 wt dt +. 


z 
(18) 


T n 
+igq\$ Dh > (A cos wt + 


o k=O f= 


ro 





+ Th) t+n—1at4pectcpey Alpe sin wt dt. 


after transformations. 
Since series (14) is formed by adjacent pulses of 


(13), the relation 


Ty= he |S A (cos wt +- 


i=] 


+ Thi) 4-4—s) at4pr<t <n) At+6e = (19) 
n 
= >) fal(Acos wt + Ty)apn—ayarypect<ui+mat+se)- 
i=] 
is justified. 
Analogously, one can write 
n—1 


>) fal(A cos wt + TH; )¢44-4—1) at4-62<t< (+n) 41481] = 
k=0 


k=0 


= fx (A cos wt + Il). (20) 


On the basis of (19) and (20), expression (18) can be 
put in the form 


i -13 (4) fu (A cos wt + Tl,) cos wt dt+ 
” tas 
me (21) 
+f qa\ fu(A coset + Th) sin ot de 
0 
or 
eS | A, 
av a et ( Thi), (22) 
where 


1 
- 
» (A cos wt +- Tiesa-oaisnrcicesmaies er 





I (A, 1,) = 
T T 
iS i3°f 
= 74 \/n(Zn)) cos wt dt + Fa) Ju (Zn,) 8in wt dt, 
0 0 (23) 
Zn, = Acos wt + II. (24) 


In this case, I(A, T1,) is the equivalent amplifica- 


tion factor of the nonlinear element, calculated under 
the condition that the function Z= 27,9 at the output 


of the nonlinear element is defined by expression (24). 
Let us assume that X(0 is approximated in such a 


manner that it has different heights m during period T 

and that their amplitudes are equal to 
ll.(q=1,2,3,..., m). 

The term dt, will denote the sum of time intervals At 


long, which occur in one period r , during which the 
function X69 assumes the value Tq. 


Then, saying that n= r/At, expression (22) will be 
written as 


by= DATA), 8 
q=1 


where Tg can equai zero for one of the values 


q=1,2,3, soe, M, ae 
If f,,(Z) has even symmetry, and the term X_,(t) 


in Z(t) has such a form that for every value of II; > 0 
during a period there is a Il; < 0 of equal absolute mag- 
nitude, one can easily see that expression (25) assumes 
the form 


= a 
= >) 2— I (A, II,). (26) 


q=l 


Here the value of q corresponds to the sections 
where Tg = 0 or 1, =0 for the function Ky(0s 


o =(m/2) for even m and o=(m+1)/2 for odd m. 

for which Tg is zero (m odd) is taken n when 
4ty is opm Pendant during period tT. The period 
= the function X 19 was broken into n equal parts in 


determining the mean equivalent amplification factor; 
time shifts of X(t) by k At were considered, where 


k=0,1,2,3, ...,(m—1)(n values in all). 


Let us divide At into \ equal parts, so the period 
T will be divided into nA equal parts. We will define 
the reduced mean equivalent amplification factor I, , 
for such time shifts of the function 
Xyplt-k4t—(e/A) At), 








z 





where a series of intermediate values e=0, 1, 2,..., 

X-—1 are assumed at each value k=0, 1, 2,...,(n-—1). 
Since expression (25) does not depend on the phase 

shift of XA , one can easily see that the expression 


for the unknown mean equivalent amplification factor 
will be the same as (25). 

As \— @, the expression for = is defined as the 
mathematical expectation, since any phase shift of 
X(0 is equally probable, Therefore, the equivalent 


amplification factors developed above may be consid- 
ered not only as the arithmetic mean of all possible 


discrete initial phase shifts of XA (0, but as the mathe- 


mnatical expectation also. 
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In order to get the expression for the mean equiva- 
lent amplification factor for a continuous function, At 
must be made to approach zero as a limit in formulas 
(25) and (26). The expression (4t,/r) appears in (25) 
and (26); it depends on the value of ni, , and is «a dis- 
crete function. 

As At-+-0,n->, Tl, 1, and the ratio At,/t 
tends toward a continuous function: 

(At,/r) -W (I) dll. 
In this case, expression (25) assumes the form 





}-od T 
Ly =75 \ W' (II) att | fu (A cos ot + Il) cos wt dt ++ 
—0o 0 


(27) 
+00 T, 
tidy \ W (II) amt | fn (A cos wt + Il) sin wt dt. 
—oo 0 
It follows from the formulas obtained that the re- 
duced equivalent amplification factor does not depend 
on the frequency of the noise, but is determined by its 
form, Hence, W(II) can be treated as the probability 
density that the function Xy(0 assumes the value II 
during a period, Therefore, if a random stationary pro- 
cess, whose frequency spectrum lies beyond the the pass 
band of the system, is considered to be the noise, then 
the formula obtained would hold, but the probability 
density W(II) in this case must be preliminarily com- 
puted from the given random process, 
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Fig. 3, 
Elementary terms of the type I(A, Ti, )= NIo(A, Th) 


appear in formulas (25) and (26); expressions for these, 
which pertain to certain types of nonlinear elements, 
are presented in the Table. 

It is convenient to characterize the quantity I_ by 
the ratio 6= I1_/a, 9 

When systems are graphically analyzed, one can 
use a graph of the function .* NI af (A), whose 
form depends on X7(0. 


The graphs can easily be obtained for any type of 
step disturbance Xo or random function, whose pro- 


bability density is approximated by discrete functions, 
by means of a family of curves of Ig(A/a), plotted for 
elementary terms, where 8 = 8;, B,, B3,. . . = Const. 
One of these families, plotted for the case when the 
nonlinear element has the characteristics of a linear 
amplifier over a limited range,is shown in Fig. 2. 
When X71 (9 is shaped as in Fig. 3, where Tr is the 


period, II the height, At the step width of Xn loay™ 


these families are coincident with the families of |, 
at B= B;, Bo, By, ... = Ccomst. 

We will get one curve for the dependence of the 
mean equivalent amplification factor of the nonlinear 
element for a given disturbance. 

The analysis of stable self-oscillation in a disturb- 
ed system is done by means of a graph of the disturb- 
ance and a graph of the linear part of the system in the 
same manner as for an analysis of system stability with- 
out perturbation. 
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ON THE CONTROL OF CERTAIN INDUSTRIAL PROCESSES 
BY MEANS OF MINIMUM CHANGES OF THE CONTROLLING INPUTS 
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The problem of developing automatic control systems which ensure the required transient response with 
minimum change of the controlling inputs, necessary only for compensating for disturbances (without measuring 


the latter), is discussed. 


It is shown that this problem can be solved by using automatic control systems which have a variable 


structure, 


A nonlinear industrial electrical regulator with a variable structure, which has been developed, is described 
as an example of one of the many practical solutions of the stated problem. 


In the control of certain continuous industrial pro- 
cesses with interrelated parameters, the best technical 
and economic features are achieved when the control- 
ling inputs need only be changed, in accordance with 
the static characteristics of the systems, enough to com- 
pensate for disturbances, We will call these changes of 
the controlling inputs minimum. 

In the use of existing industrial regulators, it is, in 
principle, possible to have transients for which the 
changes of the controlling inputs will not exceed those 
necessary to compensate for disturbances. However, 
these modes are realized only for very small amplifica- 
tion factors in the control circuit, and this causes a 
large dynamic error and increases the control time. 

The transient response of systems having delay is espe- 
cially bad when the above-mentioned limitations on 
the magnitude of the controlling inputs are imposed. 

Hence, a contradiction arises between technologi - 
cal requirements and the dynamics of control in the use 
of existing industrial regulators for a number of systems; 
for technological reasons, it is best to minimize the 
given value of the controlled parameter by means of 
minimum changes of the controlling inputs, but the er- 
ror will be maximum when the change of controlling 

input is minimum, according to the conditions of the 
dynamics of control, 

: It is for this reason that it is necessary to devise au- 
tomatic control systems which ensure the required qual- 
ity of control for minimum changes of the controlling 
inputs, necessary only for compensating for disturbances 
in the static modes (without measuring the latter). 

Let us clarify the necessity for stating such a prob- 
lem by a series of technological examples. 

When the hot mass does not burn completely in fir- 
ing in the boiling layer, it is best to control the tem- 
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perature of the boiling layer by changing the consump-~- 
tion of air (the draft), both from a transmission coeffi- 
cient standpoint, and from a standpoint of the dynamic 
characteristics of the process. 

However, a change in the draft also strongly affects 
the hydrodynamics of the boiling layer. This effect 
is especially strong in the treatment of polydispersed 
materials obtained in industry; the range of change in 
the draft is extremely small for these materials. 

Stable “boiling” cannot be achieved if the draft 
velocity is low. Stable boiling occurs when the ef- 
flux of small fractions reaches 30-50% of the load. 

If the draftis further increased, the efflux of small 
particles from the voiling layer increases sharply, and 
at a certain time its quality begins to get worse. 

The air draft at which the quality of the small par- 
ticles which are carried off begins to get worse corre- 
sponds to the ultimate (maximum) productive capacity 
of the apparatus (furnace); the most intensive (stable) 
boiling is also ensured at this draft. 

Any change in the draft will either lead to a reduc- 
tion in the yield of the furnace and intensity of boiling, 
or will cause the quality of the small particles being 
borne off to worsen. In view of this, it is clear how im- 
portant it is to ensure that disturbances are compensated 
for by minimum changes in the draft, in the operation 
of an automatic system which controls the temperature 
of the boiling layer. 

In controlling the boiling- layer temperature when 
the hot mass burns completely (firing of zinc concen- 
trates, pyrites, etc.), the changes of the controlling 
flow (load) which are necessary for a stable tempera- 
ture in static modes does not exceed +5—8%, In this 
case, the fluctuation in the concentration of SO, in the 
burned gases occuring in the sulfuric acid industry 








ar 





should not exceed + 0.3-0.6%, which is completely per- 
missible. 

The most important thing in the reliable use of 
straight-flow evaporative elements, which are used to 
control the heat output from the boiling layer, is to 
have a mode in which the changes in the aqueous lead 
of the element are minimum, since it is known that 
straight-flow evaporative elements do not work well 
when the load changes sharply [1,2]. 

In controlling the temperature of flame-heated, re- 
flective, Martin, rotary, and other furnaces, it is best to 
minimize the changes in fuel and air consumption [3]. 
It is known that the greater the change in the consump- 
tion of fuel and air, the more the nydraulic mode of 
the furnace deviates from optimum. This causes the 
thermodynamic performance to become worse, and 
complicates the operation of the control system because 
of a displacement of the flame root with respect to the 
pulse~selection points. 






































discharge density of classifiers with minimum fluctua- 
tions in the pulp speed in the classifiers, controlling 
the superheat of steam with minimum variations in the 
consumption of the injected water, etc. 

The problem stated was successfully solved by the 
development of a new regulator which was based on 
the theory of nonlinear automatic control systems with 
variable structure [4—7], 

Let us examine a very simple automatic control 
system as an example illustrating the synthesis of a 
nonlinear control law which ensures the necessary 
change in the structure of the regulator; this system is 
described by a nonlinear differential equation of the 
second order, The block diagram is shown in Fig. la. 

Without the nonlinear device which changes the 
structure of the closed loop, the control system is de- 
scribed by the following linear differential equations; 

equation of the controlled object: 

































One could give analogous examples from a num- Tt + pe = — +f (t); (1) 
ber of othertechnological processes — controlling the 
“4 f(t) 
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%, (9.9) 
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equation of the master element : 


t— s(t) =4%, (2) 
equation of the astatic regulator: 


veer 
T Mt = - F- (3) 


Here x is the relative change in the coordinates 
being controlled, ¢ is the relative change in the error 
signal, w is the relative movement of the controlling 
mechanism, f(t) is the disturbing input, g(t) is the mas- 
ter input, T,/p is the time constant of the controlled 
object, Ip is the static transmission coefficient of the 
object, T, is the integration constant of the astatic 
regulator, and 1/6 is the static-transmission coefficient 
of the sensing device. 

Let us synthesize the nonlinear control law of the 
linear automatic control system being examined, with 
changes of the controlling inputs minimized, We will 
use the method of multileafed phase surfaces to derive 
this nonlinear control law [8]. Then the problem re- 
duces to one of constructing a phase portrait of the non- 
linear system, which reflects the succession of structure 
changes necessary to meet the requirements formulated 
above. 

Eliminating the variables x and » from equations 
(1), (2), and (3), and setting f(t)=const, and g(t) = const, 
we get a second order linear differential equation 





9 + 2op + wey = 0, (4) 
where 

‘ 2 1 

ob = 7» ane ee 


Three quite easily realizable system structures are 
possible for a given object structure, when other influ- 
ences on the control law (derivative, integral, etc.) 
are absent, 

1, The object is in a negative integrating feed- 
back loop. In this case, the system is described by 
equation (4). 

2, The object is in a positive integrating feedback 
loop. In this case, the system is described by the equa- 
tion 


9+ 2p — wp =0. (5) 
3. The system is opened at the error signal chan- 


nel, In this case, the equation for the control system 
has the form 


9 + 2bp = 0. (6) 
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Eliminating time from equations (4), (5), and (6), 
we get an image of the control process on phase planes 
? and ¢, 

It can be shown that only the two system structures 
described by equations (4) and (6) need be used to real- 
ize the control process with the restrictions imposed 
earlier. In this case, the relation between coefficients 
2b and wi in equation (4) is chosen so that the process 
will be oscillatory. 

The phase plane for equation (4) will be filled by 
a family of spirals (Fig. 1c), The equation for the spi- 
rals will take the form: 


(G+ bp)" + wtp = C oxp| > (a+ arctg etie)|, 


i? 


(7) 


where w, = 4w} =b*, and C is an arbitrary constant. 
The phase plane for equation (6) will be filled by 
a family of straight lines (Fig. 1b), whose equation is- 


9 = —2bp + A, (8) 


where A is an arbitrary constant. 
Let us plot a phase plane, so that if(To+ Kj) 9)9=0, 


where T and K; are defined by the ratio 2b =(K;/T,; ), it 


will correspond to a phase plane of a system with a 
structure described by equation (6) (leaf I), and if 
(T,9+K 9)9> 0, it will correspond to a phase plane of 


a system with a structure described by equation (4) (leaf 
II). In this case, the new nonlinear control system, 
whose phase plane is shown in Fig. 1d, meets all the re- 
quirements stated above. The process goes on without 
overcontrol with a small error signal ¢. The controlling 
mechanism is deflected only by an amount necessary to 
compensate for the disturbing input. This follows from 
the transient curves (Fig. le) occurring in the finished 
variable-structure control system. The control processes 
have been found for T, = 7000 sec, T,=50sec, and 
p =1, 

Hence, the equation of a nonlinear regulator, which 
ensures the necessary structure change during the tran- 
sient, will have the form: 


7 | ° 4 
Tep=Fe—Lu@oeze (9) 


where 


O for (%i9+Ki—e>0 
f. (9,9) =i for if 
Saute (1 for = (T}9 + Kj) p<0 


Let us compare the properties of this nonlinear sys- 
tem with those of a linear control system based on the 
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IR - 130 isodromic industrial regulator and consisting of 
the static and astatic parts with transfer functionsK,, and 
1/T,p, respectively. 

By means of simple transformations, the IR-130 reg- 
ulator can be made into an integral regulator with de- 
rivative action (Fig. 2a). The derivative coefficient T 








f(t) 









































is related to the constants of the static and astatic parts 
by the relation T= T, Key. 

The comparison is made for systems which operate 
as stabilizers; consequently, the maximum dynamic er- 
ror is the main index of comparison, 

The parameters of these automatic control systems 
are chosen so that they will not have overshoot, and so 
that the movement of the controlling mechanism will 
be minimized. Since the coefficients p and 6 are set 
equal to unity, the coefficient in the linear variation 
must not exceed unity, on the basis of the restrictions 
imposed on the movement of the controlling mechanism. 

Since the transient is not oscillatory, the integra- 
tion constant of the regulator T, is related to the time 
constant of the object T, by the equation: 


4 
t= 34 K_7 Te (10) 


The maximum error for the linear system is found 
from the formula 




















‘ 2(p + Kg) 7,8 
= 4 sip 
Pmax %o / ToT bexp 7 (47, —(p + Kgy)*7,) 7,8 
(are SYVUT,—otk, TAT? 2 
(11) 


and for the nonlinear system from the formula 





2pT 6 














ae pT’, ~ 
= f T,7 8 ~ t : — ae ‘ 
¥ max Po VY _ ly (47, — p*7’,5) Tb (are 8 V (47 ,— p*7',) Td 2 )| (12) 


where % =f(t)/T, . 


The formulas were derived by means of point trans- 
formation, by reflecting the +¢ half-line into the +¢ 
half-line. 

Graphs of the dynamic error for an object without 
(1), and with (2),differentiation are shown in Fig. 2b; 
the derivative action was chosen so as to result in the 
least values of ¢ for minimum movement of the con- 
trolling mechanisii. Curves for an object with a non- 
linear regulator are shown also in Fig. 2b for T, =50 sec. 
(3), and T,=1sec (4). 

A comparisun of the results obtained shows that in a 
system with nonlinear correction, the error can be made 
adequately small as the rapidity of regulator response is 
increased, in contrast to a linear system, where the er- 
ror magnitude is constant for different values of T, - 
Such a response cannot be achieved in a linear system, 
since T, is related to T, by formula (10). Moreover, the 
error in a linear system decreases as the time constant 
of the object increases. 





It is difficult to obtain “pure” (undistorted) differ- 
entiation for the formulation of a nonlinear control law 
as applied to slowly occurring processes, In making a 
linear actuating mechanism, special linearizing ele- 
ments are built into the regulator; in particular, thermo- 
elements are used in the IR- 130 regulator. Therefore, 
it would appear wise to build such apparatus by using a 
nonlinear regulator with a variable structure, in which 


it would be possible to eliminate the derivative and lin- 
earizing elements of the actuating mechanism. The 
circuit of a regulator which has been developed is shown 


in Fig. 3, and anover-all view of the experimental mod- 
el in Fig. 4, The regulator consists of logic circuits, in- 
corporating a constant velocity integral link, which en- 
sures, in conjunction with relays Pj, BR, Py, &, %, and R, 


the necessary structure change of the system during the 
transient. The regulator includes a constant-velocity 
integrating actuating mechanism. 
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Fig. 3. Nonlinear variable structure regulator; 

(a) block diagram (T, time constant of the IM 2/120 actuator, T— time constant of 

the D - 32 capacitor motor), K;, K,—constant acceleration factors, v — feedback factor 

of the correcting circuit; (b) schematic diagram: R,- R,=470 ohms (type SP - 2), 

Rg - Rg = 470 ohms (VS- 05), Rs = 390 ohms (VS), Rg = 1.5 kohm (SP-2), Ry = 150 ohms 
(VS-1), Rg = 120 ohms ( VS-1), Rg= 1.2 kohm (wire-wound), Ryg-Ryy = 60 ohms (wire- wound), 
Rp Bro Rr3~ Ryp4 = 4000 ohms(PEV-12), Rrs5—Rrg= 1800 ohms (PZ-15), transformers Tr,. 
Trg, Trg, and relays P, and P, are standard units from the IR-130 regulator; relays Ps, P,, 
Ps, and Pg are type RPT-100; CB 1 and CB 2 are control buttons; US-1, and US-2 are 
universal switch gangs, LS 1 and LS 2 are limit switches. 
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Fig. 4. Over-all view of the nonlinear regulator, A,~adjusting 
knob for the T coefficient; A,—adjusting knob for the Ky and Ky 
factors; Py and P, are the error channel and correcting loop channel 
relays; Ps, Pg, Ps, and Pg are type RPT - 10 relays. 


The experimental data obtained on the transients 
showed that the regulator ensures a high control quality 
of stable objects without delay, and of objects having a 
relatively small delay, with minimum changes in the 
controlling inputs, which are practically equal to the 
amounts necessary to compensate for disturbances. 

If significant delay is present in the object, the 
structure of the logic networks, and also other system 
characteristics, are chosen in accordance with the dy- 
namic properties of the object. 

The regulator successfully passed industrial tests 
and is being used experimentally, The test results per- 
mit one to predict that this type of regulator will find 
use in the automation of continuous industrial processes 
in which control should be accomplished with minimum 
changes in the controlling inputs, because of technologi- 
cal reasons, and also in cases where high static and dy- 
namic control accuracy must be ensured without intro- 
ducing derivatives of the parameter being controlled in- 
to the control law. 

A theoretical analysis of a nonlinear regulator with 
variable structure ,which was actually built,is an inde- 
pendent problem and will be considered separately. 

The electrical type of nonlinear variable-structure 
regulator, which was described above, is one of many 
possible constructions based on the principle being con- 
sidered, namely, a method for construction systems 
which ensure a high control quality with minimum 
changes in the controlling inputs, This principle may 
also be applied to pneumatic, hydraulic, as well as other 
electrical apparatus. 


M, S. Lodyseva and O, P, Zelentsov took part in 
this work, 
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A method is proposed for calculating a program of discrete calculating devices for linear servomechanisms, 
based on the use of statistical results on possible external influences. The method ensures that the systems ob- 
tained will satisfy the "roughness" condition in the presence of unstable elements in the control circuit. 


Special importance has recently been attached to 
methods of calculation of automatic control systems and 
elements, based on the theory of stationary processes. A 
significant contribution to the development of these meth- 
ods has been made by the Soviet scientists A, Ya. Khin- 
chin [1] and A, N, Kolmogorov [2]. There is a good de- 
scription of the theory in a publication by A. M. Yaglom [3]. 


The application of the theory of stationary processes 
in technology was initiated by N. Wiener. Wiener's in- 
vestigations in [4] played an important part in the de- - 
velopment of practical methods of designing selective 
and smoothing filters, which are important com- 
ponent parts of many automatic systems, Statistical ineth- 
ods in the designing of various technical devices have 
recently been developed in the works of V. V. Solodov- 
nikov [5], V. S, Pugachev [6], Ya. Z. Tsypkin [7], L. T. 
Kuzin [8], and others, 

Calculating devices of the discrete type (pulse fil- 
ters and digital machines) have recently been used in 
various technical systems, In the present work we will 
consider the problem of the synthesis of servomechanisms 
containing such devices, A peculiarity of this problem 
lies in the fact that only part of the system is synthesized 
—the correcting element supplied by the discrete calcu- 
lating device, while the object of the control is given. 

It is necessary, in the synthesis of the correcting element 
to the object, to consider the possibility of destroying the 
“roughness” of the system by a divergence between the 
parameters of the correction link and the controlled ob- 
ject, which leads to the accumulation of errors and in- 
stability. 

For systems with a stably controlled object, the prob- 
lem of synthesizing a discrete calculating device ,and 
taking into consideration the “roughness” condition, was 
solved by S, Chang [9]. Here we will solve the problem 
for the control of an arbitrary object, stable or unstable. 
In contrast with Chang's method, which involves labori- 
ous calculation, the proposed solution is based on apply- 
ing the method of equations involving polynomials de- 
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veloped in [10]. The solution of these equations is noi a | 
mathematically difficult operation, and can be compiete- 
ly carried out by hand, The problem we consider reduces 
to a system of two equations of this type. A method of 
solving these equations involving polynomials is given in 
the Appendix. 


1.Some Results from the Theory of 





Discrete Stationary Processes 





The results that we give here are available in one 
form or another in [3, 4, 7, 8, 9]. The various authors 
have, however, used different operational methods. The 
method used here also differs somewhat from those used 
previously, 

We will use the delay operator z as the basic opera- 
tor, where this operator is related to the differentiation 
operator p and the delay interval T by the equation 
z=exp(—pT). The use of the physically realizable op- 
erator exp(—pT) instead of exp(pT) simplifies the de- 
scription and design of the dynamic systems. 

A stationary discrete process will be a sequence of 
random variables X={x,} (i=— @,..., +), possessing 
the following properties: 

1, the mean value of the random variable x; does 
not depend on the subscript i, and for simplicity we will 
assume that this mean value is zero; 

2. the correlation coefficient of the random vari- 
ables x; and x; , , depends only on the difference k be- 
tween the subscripts; we ‘enote it by 

Wk) =XiX; , ,(kK=-™,... 


The quantity (0) “3 is called the dispersion of 
the stationary discrete process, 

The double infinite z-representation of the sequence 
of these coefficients will be called the correlation-repre- 
sentation of the process 


» +), 


+00 
¥(z)= > p(k) 2%, () 


k=—0oo 





\“y 





Since the correlation coefficients are even functions 
of k, the correlation-representation of the process pos- 
sesses the symmetry property relative to the unit circle 
described by the relation 


f(z) = F(z"). (2) 


A sufficient condition for the convergence of the 
series (1) is 


be 1) 


(k > M). (3) 
If this condition is satisfied, then the function ¥ (z) 
is analytic in the ring 


1—é<|z|/<1+6 (>0) 


and (1) is the Laurent expansion of the function in this 
ring. The correlation coefficients ~ (k) can be expressed 
in terms of ¥(z) as the coefficients of this Laurent 
expansion 





(H=a, § "OS 


| > {=1 


(4) 


We will limit ourselves to the consideration of pro- 
cesses with correlation-representations that are rational 
functions of z. 

The following method is proposed for obtaining the 
correlation-representation in terms of the sequence of 
correlation coefficients. We first of all sum the part of 
the series (1) with nonnegative powers 


WY (z) = >) o(k)z*. (5) 
k=0 


Condition (3) ensures that this series can be summed 
for all | z|<1+5. We obtain as a result a rational func- 
tion in the form 


K 
(2) = egy (6) 


where K(z) and U' (2) are polynomials in z, and the 
polynomial U (z) has no zeros inside the unit circle. 
The part of the series (1) with negative powers of k 


¥_(2)= Di g(k)z* (1) 


k=] 


can be obtained from the part with nonnegative powers 
from the relation 


¥_(z) = '¥, (2) — (0). 


It is also a rational function of z 


(z) 
¥_(s) = - (z)’ (8) 





where Q(z) and U"(z) are polynomials in z, and the 
polynomial U (z) has no zeros outside the unit circle. 

The function ¥(z) is obtained as the sum of (6) and 
(8) 


¥ (z) = ¥_(z) + Y,, (2). (9) 


This function can be reduced to the form 


_ 1%) 
F(2)= 7G" (10) 


where I(z) and U(z) are polynomials in z . 
The factorization of this function 


Y (z) = ¥* (2) ¥> (2) (11) 


is carried out by factorizing its numerator and denomin- 
ator: 

I (2) = I* (2) I-(2), (12) 

U (2) = U* (2) U~ (2). (13) 


We introduce the notation 





y | = 
Yj=7o. 4 @=74. ay 


For the addition of statistically independent pro- 
cesses, the correlation-representation of the sum will be 
equal to the sum of the correlation-representations of 
the two components, 

When a stationary discrete process passes through a 
stable system with transfer function H(z), the output is 
also a stationary process with correlation-representation 
[9]: 


VY, (zc) = A(z) A (2?) ¥; (2). (15) 


2. Statement of the Problem 





We will consider a servomechanism with a discrete 
calculating device having the transfer function 


— —W (2) G(s) 16 
1(@)= > Tweew ’ sr 
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where G(z)=[P(z)/R(z)]isthe transfer function of the 
given part of the system, W(z)=[C(z)/D(z)) is the line- 
ar program of the discrete calculating device, and P(z), 
R(z), C(z), and D(z) are polynomials in z. 

The input of the servomechanism is the sum of two 
processes ~the useful signal S(z) and the disturbance 
N(z) 


X (z) = S(z) + N (z). (17) 


The program of the discrete calculating device 
must be chosen so that the output of the system 


Y (z) = H(z) X (z) (18) 


reproduces as well as possible the given function of the 
useful signal 


V (z) = A(z) S(z). (19) 


The operator A(z) characterizes tracking or track- 
ing with lead 


A(z)=2z-* (=0,1,2,..--) (20) 


We will assume that the signal and the disturbance 
are stationary statistically independent processes with 
. Zero means and the given correlation coefficients 


bs (k) = SiSitky dn (k) = MWNi+k- (21) 


We construct the correlation-representations of these pro- 
cesses, and factorize them: 


If (2) IF (2) 
Ut (2) U5 (2) ’ 
IK (2) Ay (2) 
Ut@UR«) © 
We calculate the correlation-representation of the 
sum of these processes and factorize it: 


Y.(z) = ‘V,(z) + Py (2z)= 
It (2) IZ (@) 
Ut (s) UZ (@) UF (2) UN (2) © 


The polynomials U,(z) and Un, (2) are assumed to 
have no common factors 


¥,(z) = 





(22) 





Vn (z) = 


(23) 





(U,, Un) =1. (24) 


The error of the system is the difference between 
the actual and the desired signal at the output, The re- 
presentation of the error signal has the form 
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E=Y—V=H(S+N)—AS*= 


(25) 
=(H—A)S+ HN, 
If we take into account the statistical independence 
of the signal and the disturbance, and use the rule (15), 
we obtain the correlation-representation of the error 


Y.=(H—A)(HW—A)¥,+HAVy, (26) 
where for brevity we have used the notation H= H(z74y 
A= A(z". 

As a measure of the quality of the system, we take 
the dispersion of the error J= ¥,(0), which is obtained as 
the zeroth coefficient of the Laurent expansion of ¥, (Z) 
and which is equal to 


(27) 


In choosing the program for the discrete calculating 
device, it is necessary to ensure that the following con- 
ditions are satisfied: 

1. The synthesized system must be stable; this 
means that the function H(z) must have no poles in the 
region | z| = 1. 

2. The program of the discrete calculating device 
must be physically realizable; this means that the func- 
tion W(z) must not have a pole for z= 0. 

3, The “roughness” condition must be taken into 
account in the synthesis of the system, i.e., the condi- 
tion that the reduction of the zeros and poles of H(z) lo 
cated in the region |z|< 1 must be excluded, 


4. The synthesized system must be optimal, i.e., 
a minimum condition must be satisfied for the disper- 
sion of the error 


" 1 . dz 
J = ) a4. (2). 


|zj=1 


(28) 


In the following section, we give an algorithm for 
choosing W(z) so that it satisfies these conditions. 


3, Algorithm for Choosing a Program 
for a Discrete Calculating System 








The algorithm is as follows; 
1, Factorize the polynomials P(z) and R(z) 


P=FP, R=a RHR. (29) 

2. Form the system of polynomial equations 
UnP04+R-O® = It, (30) 
o,P0+4+U,M=1Uy. (31) 








We solve this system for the unknown polynomials 
6, ®, and Il, taking into consideration the condition of 
physical realizability ©(0)# 0 (see the Appendix to this 
paper and the Appendix to [10)). 

We then choose a solution in which the lowest pow- 
er of @ occurs, 

3. We find the program for the discrete calculating 
device from the formula 


Un Re 


, . 
W =. 


(32) 
We will show that the above algorithm yields a so- 


lution that satisfies all the stated conditions, If we sup- 
stitute (32) in (16) and use (30), we obtain 


UnP-o 
= r 





; (33) 


The expression (33), obtained without reduction of 
the-zeros or poles in | z|< 1, satisfies the “roughness” 
condition. Since t has no zeros for |z| =1, the system 
obtained is stable. 

The condition W(0)< © for physical realizability is 
satisfied by the method of solution of the polynomial 
equations. It remains to prove that the system is opti- 
mal. 

We find for the variation of the correlation-repre- 
sentation of the error (26) the expression 


oY. = (HV,— AYV,) 3A + (A¥,—AY,) 8H. 


We now calculate the integrals of each of the terms 
of this expression 


r= $ (#¥,—Av)3AaS , 
[zj=1 
y (AY, —AY,)sH© 


jzl=1 


r, = 


It should first of all be noted that the second inte- 
gral is reduced to the first by the substitution u=z~*, 
Thus, to prove that the system is optimal, it is sufficient 
to prove that [;=0. If we substitute in I’, the expressions 
in (20), (22), (23), and (33), we obtain 


21> P-9—1,Uy 


* dz 
P. - cw. 
: 2“U+U-UZ 2 


[z|=1 





Using (31), we can write this expression in the form 


r,=- $ ——*_3AF. 
tz 2°U, UN 


The value of this integral is equal to the sum of the 
residues of the integrand at its poles outside the unit cir- 
cle, This function, however, has no poles for |z|>4, 
and soly=0,. The validity of the algorithm is thus proved. 


In the case when the given part of the system is 
stable 


f=1, (34) 


the calculation of a program for a discrete correcting de- 
vice is somewhat simpler. In this case, condition (30) 
imposes no limitations on 9. Therefore, when 6 has 
been obtained from equation (31), it can be substituted 
in (30) and the corresponding ®@ can be calculated. The 
program W(z) is thus obtained more quickly in this case. 


SUMMARY 


The condition that the “roughness” of the system 
be maintained, when calculating devices are introduced 
into it, has a serious effect on the method of calcula- 
tion of the system. The theory shows that the introduc - 
tion of calculating elements in servomechanisms can, 
in all cases, be used to obtain optimal systems, satisfy - 
ing the “roughness” condition, but in the case of un- 
stable objects of control,the system will not be as good 
as in the case of stable objects of control. The control 
program for unstable objects also turns out to be more 
complicated. In order to obtain a high-quality servo- 
system, it is desirable first of all to make the object 
being controlled stable at the expense of supplementary 
internal coupling. In those cases where the introduction 
of this coupling is impossible, we can have recourse to 
stabilizing and optimizing the system by a discrete 
calculating device. 


APPENDIX 

The solution of the system of polynomial equations 
can be reduced to the solution of one polynomial equa- 
tion. Assume that there are two equations: 


A,X aa M,Y,= B,, (D 


A,X + MsYo= B, (Tl) 


in the unknown polynomials X, Yj, Yg. We will assume 
that the polynomials M, and M; have no common fac- 
tors 


(My, M,)=1- 
We will show that this system is equivalent to the 
“combined” system 


(A\M2+ AgMy) X + MyM2Y - (ByM,4+- ByM)). (1) 
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This equation is obtained by reducing the original 
equations to the samme modulus M= MyMz and then add- 
ing them, We denote Y by the sum of Y; ard Y,. 

It is clear that every solution of the original equa- 
tions will be a solution of the combined equation. 
We will show that the converse is also true: every solu- 
tion of the combined equation is a solution of the 
original equations. Inorderto do this, we write equation 
(III) in the form 


(A,X — By) My+ (AX — B2) My=—M, MY. 


Since the right-hand side of this equation is divisi- 
ble by My, and the second term on the left-hand side is 
also divisible by My, then the same is true of the first 
term on the left-hand side. But, since M, and Mg, are 


relatively prime, then (A;X~B) must be divisible by M,. 


We thus obtain (I). Equation (Il) is obtained in a simi- 
lar way. Thus the equations (I) and (II) follow from 
(III), which means that a solution of (III) will also be a 
solution of (I) and (Il). 
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Two digital devices, based on the principle of integrating voltage pulses, are considered. 


At the present time, in addition to digital comput- 
ers based completly on the application of devices using 
crystal semiconductors, there are also many computers 
using a combination of the latter with magnetic ele- 
ments with cores having rectangular hystersis loops (fer- 
rotransistor elements). 

One of the directions being taken in the develop- 
ment of ferrotransistor elements is toward the creation 
of junctions and elements based on the principle of inte- 
grating voltage pulses. The principle of integrating vol- 
tage pulses by using saturated magnetic cores has been 
known for a long time. It has been applied in designing 
magnetic amplifiers, pulse transformers, and other simi- 
lar devices, 

In digital technique, the first development of ele- 
ments based on this principle was relatively recent [1, 
2). 

We will describe briefly the basis of this method. 

If the input of an R-L electric circuit with a satur- 
ated magnetic core (Fig. 1a) consists of pulses of arbitrary 
form from a voltage pulse generator with low internal resis- 
tance, and if the energy of a single pulse is not suffi- 
cient to cause a complete magnetic pole reversal of the 
core, then the magnetic state of the core-material will 
vary by steps, and change from a state of negative mag- 
netic saturation to the opposite state of positive magnet- 
ic saturation, In this case the core accumulates or inte- 
grates the pulses arriving at its input coil (Fig. 1b). 

If it is assumed that the operating resistance R of a 
circuit containing saturated magnetic cores is so small 
that it can be neglected, and that all the voltage U is 
available for magnetic pole reversal of the core, then, 
starting from the equality 


dD 
U =s €== Ww, a 9 
we can obtain known expressions for the variation of the 
induction 4B due to one voltage pulse, and for the maxi- 
mum number, n, of pulses that can be accumulated in the 
integrating core: 


t 
on 
0 


2B,W; S 
Ws’ ” 


a Rua : 
0 


Here W, is the number of turns in the integrating 
coil, ® and B, are the current, and the residual induc- 
tion of the core-material, B,, is the maximum induction 
of the core-material,e the emf arriving at the coil, S 
the cross section of the core, and tb the duration of the 
pulse. 

It is evident, from the approximate expression for 
the maximum number of pulses that can be accumulated 
by the integrating core, that this number is given in 
terms of the amplitude and duration of the pulses, the 
number of turns in the integrating coil, the cross section 
of the core, and the residual induction of the core-ma-~- 
terial, 

A more accurate and complete expression for the 
number of pulses accumulated in the integrating core 
can be obtained from the expression for the energy com- 
municated to the core during its pole reversal. For a re- 
view of the essential principles and basic properties of 
elements constructed on these principles, however, this 
question need not be considered, 


2 je 


0 '-6, 

Fig. 1. (a) Electric circuit with saturated mag- 
netic core; (b) hystersis loop for the core-ma- 
terial, 


AB = n 


















Up to the present time, there have been investiga- 
tions of two types of elements based on integrating vol- 
tage pulses. The first type is based on the accumulation 
in theintegrating core of a definite, previously given 
number of pulses, and the subsequent delivery of a signal 
to the output when the integration has been completed. 
These elements are called pulse-accumulators or inte- 
grators. They can be widely applied as frequency-divid- 
ers, calculators, decipherers, and linear relays. The sec- 
ond type is characterized by being able to receive an 
input of not more than three pulses, and, depending on 
the number accumulated, yielding different output sig- 
nals, This type of element is used in the adders in arith- 
metic devices applied in building digital computers. 





1, Voltage Pulse Accum@lators 
or Integrators 





We show in Fig, 2 a circuit*diagram of an integrator, 
which explains the principle of its operation, A similar 
scheme, using material in the form of a ribbon, was 
first proposed in 1955 in[1]. We will describe briefly 
the principle of operation of the integrator. It consists 
of two basic parts: a pulse-former and an integrating el- 
ement, The action of the pulse-formers is to generate 
pulses of constant area (in the time-emf plane) for ad- 
missible variations in the parameters of the device and 
in the conditions in the surrounding medium. 


The purpose of the integrating element is to accum- 
ulate a definite number of pulses, to switch the circuits 
(from the integrating circuit to the regeneration circuit), 
return the core to its original state of magnetic satura- 
tion, and transmit an output signal with energy equal to 
the sum of the energy absorbed by the core during the 
accumulation period, 

In Fig. 3 and Fig. 4 we show the two most common 
integrating circuits, 

Certain modifications are possible in the integrators: 
the introduction of a bias in a return circuit to a battery 
or condenser, an automatic bias on the regeneration coil 
(Fig. 3; 1, 2, and 3 respectively), and finally, the latter 
possibility with minimal resistance in the integrating 
circuit (Fig. 4). 

The notation in Figs. 3 and 4 is as follows; W, is 
the generating coil of the pulse-former, W, the regener- 
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Fig. 3. General circuit diagram of an integrator. 
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Fig. 4. General circuit diagram of an inte- 
grator with minimum resistance in the in- 
tegrating circuit. 


ation coil of the pulse-former, W, the feedback coil, 
W, and W, are pulse-former output coils, W is the inte- 
grating coil of the integrating element, We the regen- 
eration coil of the integrating element, W, and Wp are 
the feedback windings for the integrating and regener- 
ating circuits respectively, W, is the output winding of 
the integrator, R,, the resistance in the circuit-base of 
the pulse-former, R;- the resistance in the regeneration 
circuit of the pulse-former, R, is the circuit-base of the 
integrating element, R, the resistance in the integrating 
circuit, R, the negative-feedback resistance in the 
pulse-former emitter-circuit triode. 

All the integrator circuits used ferrite cores of type 
FP-7A and FP-7C and the triodes P13, P14, P15, and 
P16. 

In practice, a large number of pulses can be accum- 
ulated in a single integrating core. The number depends 
on the parameters of the signal generated by the pulse- 
former, and on the properties and dimensions of the in- 
tegrating element, 

Under laboratory conditions, we have been able to 
accumulate up to 120 pulses in one element with a core 
of mark FP-7 ferrite 9mm in diameter. However, when 
the reliability of operation of the element is considered, 
and imperfections in the ferrite and the triodes are taken 
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Fig. 5, FP-7C, 5mm diameter 
W, = 20 


into account, the number of pulses that can be accumu- 
lated in practice is limited. 

Figure 5 shows an oscillogram of the pulses in the 
output coils of integrators accumulating 5 pulses, 

At the present time, standard integrators are de- 
signed that can accumulate 5 and 10 pulses, the cores 
of these integrators having a diameter of 5mm and be- 
ing made of FP-7C ferrite. Their nominal collector 
voltage is15 v, with a permitted variation of + 20-25% 
from the nominal value. 

Several independently operating accumulators can 
be connected in series. In this way a large capacity can 
be attained with a relatively small number of pulse- 
counters. Each stage of such a counter works from the 
output pulses of the previous stage. The pulse-formers 
of the integrators are transitional elements between the 
stages. Each stage can be regulated by varying the re- 
sistance in the base of the triode in its integrating cir- 
cuit. Several such counters were constructed and tested. 

As an example we mention a four-stage adjustable 
counter capable of counting 600, 400, and 200 pulses, 
The second stage of this counter was adjustable for the 
accumulation of 6, 4, or 2 voltage pulses, The first, 
third, and fourth stages were not adjustable, and could 
accumulate 10, 5, and 2 pulses, respectively. The 
whole counter operated stably in a range of collector 
voltages from 11 to 22v with a nominal value of 15v. 

We will now discuss the advantages and disadvan- 
tages of these integrators. 

A basic and indisputable advantage is their econ- 
omy in the number of elements needed in their con- 
struction —cores, triodes, and resistances, For example, 
a 65536 frequency divider can be easily constructed by 
using four integrators in series, where each of the stages 
divides the frequency by 16. For this, it is necessary to 
have only 8 cores made of material with a rectangular 
hysteresis loop, 12 triodes, and 12 resistances, To carry 
out the same operation using triggers, 16 double stages 
are needed, and 32 ferrite cores, 32 triodes, and 32 re- 
sistances would be required with this method, 

A second advantage of tne integrators is the absence 
of generators of powerful variable pulses. Each integrat- 
ing stage in a counter made up of integrators operates 
independently, using only the output pulses from the 
previous stage. The number of stages of such a counter 
can be made as large as we wish. 

One of the disadvantages of the integrators is the 
high rejection rate of triodes, Triodes must be rejected 


in connection with both their collector current and the 
value of their input and contact resistance in switching 
operations, A second disadvantage is the length of the 
delay between the input and output signals. When mul- 
ti-stage counters are used, this delay must always be 
taken into account, 


2. Adders 


We will consider the essential basis of a method of 
constructing adders, using a principle of integrating vol- 
tage pulses proposed in 1958 [3]. 

In Fig. 6 are shown the main elements of the adder 
—the integrating magnetic element and the selecting 
circuits, The object of the integrating magnetic ele- 
ment is to add the items arriving at the input. Depend- 
ing on the number of pulses arriving, the polarity of the 
core of the integrating element is reversed through a 
partial cycle to one of three possible states: a, b, or c 
(Fig. 7). The quantities to be added can arrive at the 
input of the adder simultaneously or displaced in time, 
but they must arrive at the integrating element before 
the arrival of the calculating signal. This signal pro- 
duces a magnetic pole reversal to the original state of 
magnetic saturation, and at the same time causes pulses 
accumulated in the integrating element to be added to- 
gether. This calculation can be carried out in two ways: 
by three consecutive voltage pulses, the action of each 
of which is equivalent to the action of one recorded 
pulse (Fig. 8a), or by the action on the core of a single 
powerful current pulse (Fig. 8b). 

Here the emf ¢ arriving at the integrating-element 
output-coil will be proportional to the magnitude A@ 
of the variation of the current in the core, 
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Fig. 8. The calculation of information recorded in an integrating core: (a) Three suc- 
cessive voltage pulses; (b) A powerful current pulse. 


Depending on the method of calculation used, ad- 
ders can be four-cycle or two-cycle. Both these types 
record the incoming signal in one cycle, but the calcu- 
lation is carried out in one cycle in one of them, and in 
three cycles in the other. 

One or other time-sequence of pulses (in the case 
of a four-cycle adder), or a single pulse of a definite 
amplitude (in the case of a two-cycle adder), arrive at 
the input of the selecting circuits, The task of these se- 
lecting circuits is to distribute the results of the summa- 
tion arriving from the integrating element into the ad- 
der-sum and adder-transfer channels, and then to record 
the results and transfer the corresponding values to the 
adder-sum and adder-transfer outputs. 

We will consider two types of adder, a four-cycle 
and a two-cycle, 

In Fig. 9 we show the block-diagram of one of the 
possible variants of four-cycle adders, In the integrating 
magnetic element there is a Y1 amplifier, and a strictly 
saturated magnetic core, The selecting circuits contain 
a Y2 amplifier with divided outputs and two cores for 
the intermediate storage of information. 

The adder operates in the following way. In the 
first cycle, with a system of pulses I, the information 
passing through the input amplifier is recorded in the in- 
tegrating element. Simultaneously, the calculation of 
the result of the previous addition is performed from the 
selection-circuit cores. 


In the second cycle, cophasally with the system of 
pulses II, the integrating element counts a single unit, 
and records it in the sum magnetic core. In the case 
when a single unit is recorded in the third and fourth op- 
eration cycles of the adder, there will be no signal at the 
output of the magnetic element. The unit recorded in 
the sum-core B appears in the output of the adder-sum 
when a system I of pulses arrives again. 

When two units are recorded in the integrating core, 
their counting proceeds in the second and third adder op- 
erational cycles, while the unit counted in the third cy- 
cle, after passing through the amplifier Y2, reverses the 
magnetic polarity of the transfer-core A and returns the 
sum-core B to its original zero state, In this case, with 
the arrival of a system I of pulses, a signal appears at 
the output of the adder-transfer. 


Finally, in the case of the additon of three units in 
the integrating magnetic element, the counting of the 
third unit and its transcription on the sum-core occur in 
the fourth operational cycle of the adder, In this case a 
system I of pulses will produce signals at both the sum- 
and the transfer-output. This type of four-cycle adder is 
the most economical from the point of view of the num- 
ber of elements needed for it. A disadvantage of all 
four-cycle adders is their slowness of operation. 

In Fig. 10, we show one possible block-diagram for 
a four-cycle adder, Its selecting circuit consists of three 
triodes and three cores, Counting in the integrating ele- 
ment proceeds in one cycle by a current pulse, A signal 
of completely definite amplitude will be obtained at 
the output, and this signal will depend on the sum of the 
pulses arriving at the adder-input (Fig. 8b). The select- 
ing circuits of the adder perform an amplitude separation 
of the pulse arriving at their inputs. This is done by the 
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Fig. 11. 


choice of the corresponding bias values for the amplifiers 
Yl, Y2, and Y3, The quantity counted by the integrat- 
ing element and selected by the amplifiers is recorded in 
the corresponding sum- or transfer-core, 

In Fig. 11 is shown the principle of the electric cir- 
cuit of a two-cycle adder. It consists of four FP-7C fer- 
rite cores 4mm in diameter, and four P14, P15, or P16 
triodes, The adder operates at a frequency of 50-100kc. 
The signal-noise ratio (for the noise, we have taken the 
signal at the adder-sum output for the addition of two 
units arriving at the input) was not less than 10. 

As an example, we show in Fig. 12 oscillograms of 
signals at the output of the adder-sum for the addition in 
the latter of one (a), two (b), and three (c) units, Tem- 
perature experiments were carried out both for separate 
adders, and for adders in series. Stability of operation 
was checked for a series of temperatures differing by 10° 
C. The collector voltage was also varied within limits 
varying by + 20% from its nominal value. 

The experiments gave very satisfactory results, and 
it was found that the adder operated stably in a range of 
temperatures from —40 to + 55° C, and for variations of 
the collector voltage from 20% below to 27% above the 
nominal value, 

The main advantage of adders constructed on the 
principle of integrating voltage pulses, as in the integra- 
tors considered earlier, is the small numbers of elements 
needed in their construction. Adders with three inputs 
and two outputs can be built with three or four cores and 





Fig. 12. 


the same number of triodes, They are simple to con- 
struct, reliable, and permit a sufficiently great variation 
from their normal operating regime. The frequency of 
operation of the adders, and consequently their rate of 
operation, can be significantly increased, 

In conclusion, we should say that magnetic elements 
built on the principle of integrating voltage pulses would 
appear to have good future prospects in applications both 
as integrators and as adders, They can be widely applied 
both in calculating devices and in various designs of au- 
tomatic control and remote control mechanisms, 
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New, highly reliable arrangements for selection and control of objects located at different remote places 
by a binary frequency code system are discussed in this paper. It is shown that it is advantageous to use two-sec- 
tion filters to increase the effectiveness of channel utilization when the frequency method of selection is employed. 


Until recently, the methods and technical means of 
remote control were developed mainly for applications to 
large and especially important objects in power systems 
and certain other branches of industry, said objects being 
concentrated each time at a single location. 

Now, after introduction of more perfect automatic 
control methods of production processes, the problem has 
arisen of controlling from a distance a great number of 
relatively small operational units carrying out certain 
functions at different places; as, for instance, to control 
objects distributed along pipelines, railroad tracks, a can- 
al, etc.; objects located at various distant points distrib- 
uted over the wide area of an oil refinery, a gas produc- 
tion plant, a large factory, an irrigation system, etc.; or 
objects distributed throughout the great underground 
space of a mine, 

The volume of work connected with development of 
remote control of such mu!tiple objects dispersed in area 
(if for nothing else than their great quantity) will in the 
future exceed by tenfold or better the work connected 
with remote control of objects “concentrated” (lumped) 
in a single location. 

At present, remote control has already been intro- 
duced at several thousand points of oil-well operations, 

‘and a great amount of development work is in process on 
remote control of multiple or dispersed objects in the 
petroleum and gas industries, in mines, factories, city 
public utilities, etc. 

When relatively large concentrated objects in elec- 
tric power stations or substations are remotely controlled, 
a control sending or dispatching point (DP) is connect- 
ed with each of the points to which the control commands 
are executed ("Execution Points,” (EP) by one or several 
communications lines, When disperse objects are re- 
motely controlled, each communication line is utilized 
for transmission of information concerning several “exe- 
cuting” objects, which are connected to the line at vari- 
ous points; the number of such points may be very large. 
This increases the effectiveness of the communication 
line's utilization, and is the main distinction of remote 


control systems for disperse objects in comparison with 
systems for objects concentrated in one place, 

It is, therefore, natural that rational methods of in- 
formation transmission between the dispatch point and 
disperse control objects may differ substantially from 
transmission methods between the DP and concentrated 
control objects. 

“Execution” points (where the controlled operation 
is carried out) usuaily work automatically without any 
attending personnel, and, therefore, much higher de- 
mands must be made of them with respect to dependa- 
bility. It is even advantageous to introduce a somewhat 
higher or additional complexity into the apparatus at the 
dispatching end, in exchange for a simplification, and 
increased reliability, of the remotely controlled arrange- 
ments at the execution points. Thus, it is advantageous 
to employ simple and highly reliable completely auto- 
matic arrangements for operation selection at the E P, 
while using, in many cases, a system of interrogation by 
choice or a cyclical scanning system, controlled to run 
at a predetermined sequence, for transmission of the re- 
mote measurement results and of alarm signals. This 
makes it possible, when remote measurement data and 
alarm signals are transmitted back from the EP to the 
DP, not to accompany them by the address of the EP, 
which considerably simplifies the coding arrangement at 
the EP, 

Interrogration of the objects is carried out either at 
the choice of the dispatcher, automatically, according 
to a predetermined program, or may be started up by a 
signal from the E P indicating a change in the regime of 
one of the objects. The latter signal is a common one 
to all execution points (indicating an accident, a break- 
down, or other trouble at any point). The period of in- 
terrogation or scanning cycle must be much smaller (at 
least by several times) than the time needed for the par- 
ticular technological process to reach equilibrium; the 
scanning cycle may be from a few seconds to ten min- 
utes, or even more. However, needless decrease of the 
cycle period can lead to a widening of the frequency 








band needed for transmission, to an unjustified accelera- 
tion of the remote control system's operation, and to a 
lowering of the stability with respect to interferences. 

In [3] it is shown that for disperse objects distributed 
over an area, it is advantageous to use communication 
lines of the cluster (or bunch) type. In laying out the 
communication lines, all execution points comprised in 
the remote control system are arranged in groups (clus- 
ters) according to certain characteristics or arbitrary 
rules, and each such group is connected with the dis- 
patcher point by a separate communication line, There 
is an optimum for the number of such groups if it is ne- 
cessary to have, simultaneously, a high degree of struc- 
tural reliability and relatively short lengths of commun- 
ication lines, 


In [4] it is shown that to increase the reliability of 
a remote control system, the execution points must be 
connected in parallel to the communication line, just 
as in power distribution and broadcasting over wire net- 
works the subscribers are connected in parallel. 


1. Selection of Objects When Using 
a Dual-Frequency Code 








In [5] several call systems for making connections 
with disperse objects,using a single frequency assigned to 
each particular object,were compared to each other. 
Some frequency selectors which operate by means of 
electrical L-Cresonant circuits without making or break- 
ing of contacts were offered in that paper. Such con- 
tactless frequency selectors are of considerable interest 
with respect to comparatively simple systems of remote 
control having a relatively small number of execution 
points (up to 12—24) in a bunch, with communication 
lines of 15~30 km in length, and with several objects 
of remotely controlled actions and of remote measure - 
ments on call at every EP. At the present time, there 
are already in operation several thousand execution 
points at oil wells equipped with such remote controls 
using single-frequency selector systems, Frequency se- 
lector systems with electrical L-Cnetworks are beginning 
to be employed in other fields of industry as well. 

When remote control is provided for pipelines, gas- 
distribution systems, irrigation systems, and similar en- 
terprises of the national economy, the necessity of con- 
necting many execution points to a common communi- 
cation line arises and, additionally, in a number of 
cases, of operating numerous objects at each execution 
point. In these applications, operation over considerab- 
ly greater distances is also required, 

Presently existing contactless devices for selection 
and switching-on of a controlled object, such as those of 
the tuned-frequency type, those of the exponential-trans- 
former type [6], and those of the magnetic type with 
cores having a rectangular hysteresis loop (RHL), are all 
about equal to each other with respect to reliability and 
complexity. Thus, each RHL cell [2], as well as each 


tuned-frequency relay (including the electric circuit as- 
sociated with it) [5], consists of a transformer with ferro-’ 
magnetic core, of condensers, and of diodes or transis- 
tors. However, in some circuits employed with the RHL 
cells, no condensers are used, For the same output pow- 
er, the over-all dimensions and other characteristics of 
the RHL and tuned-frequency systems are about equal. 
The substantial difference between the two systems lies 
in the fact that for the first RHL cell there are only two 
operational states in which it can possibly be, while the 
frequency relay can be tuned to a large number of fre- 
quencies and can be in two different states at each of its 
operational frequencies. 


The number of all feasible discrete-selection index- 
es for a frequency relay is n times larger than for an RHL 
cell, where n is the number of discrete working frequen- 
cies available in the transmission band when single-fre- 
quency selection is used, 

Contactless selectors using a dual-frequency code 
differ still more from the other devices; the number of 
possible states of such a selector (or decoder) to which it 
can be adjusted is n(n—1) times larger than for a single 
device using RHL, while with respect to complexity it is 
about equivalent to two single-frequency relays with 
their associated electric circuits, This makes it possible 
to carry out the selection of an equally large number of 
disperse objects by a substantially simpler contactless 
arrangement, 


It should also be taken into account that when the 
signal level is limited, stability of transmission against 
interference is much better with frequency diversifica- 
tion used for multiplexing than with time separation, 
and that in connections by wire lines and in many others 
the permissible signal level actually is limited, 

Selectors (decoders) based on tuned-frequency ele- 
ments are not the only devices which can secure the 


simplest and most dependable solutions of these prob- 
lems. 


For selection of disperse objects it is advantageous 
to employ arrangements using elements which are able 
to be in numerous diverse states (providing numerous 
characteristic indexes for selection), Frequency selectors 
have such properties, but exponential transformers with 
time separation of signals also have them [6]. 

When stability against interference is of importance, 
the methods used must also be compared with respect to 
this property. 

As a result of investigations carried out at the Insti- 
tute for Automation and Remote Control of the Academy 
of Sciences of the USSR, schemes of dual-frequency codes 
were suggested and developed in which damage to any 
one of the components cannot cause a wrong s lection 
or execution of a wrong command signal, but produces 
—in the Wérst case — only a precautionary refusal to ex- 
ecute any command, 








One such scheme is presented in Fig. 1. In this 
scheme, no local power supply is required. Two frequen- 
cies, f, and f,, are transmitted consecutively from the 
dispatch point. 


A voltage divider stepdown transformer decreases 
the effects of the line on two series resonance circuits 
LyC, and L2C2. 

The first circuit is tuned to frequency f,, the second 
to fp, When frequency f; is transmitted, capacitor C is 
charged by way of the diode preceding it up to a voltage 
u, at which the energy (cu’/2) is sufficient to operate re- 
lay P. The emitter-collector resistance of the transistor 
in the other circuit is then sufficiently high. 

When the second oscillation of frequency f, is 
transmitted, the output voltage of the second circuit is 
applied either directly or through a diode to the emitter- 
base terminals of the transistor, and opens it up for oper- 
ation of the relay: specifically, capacitor C is discharged 
through the transistor and the discharge current going 
through the relay coil causes the armature to make or 
break contact, The relay may be automatically locked 
in, when this is needed, The results proved definitely 
that, when the reverse sequence occurs (first frequency 
fz, then f;), or where damage occurs to any of the com- 
ponents, false operation of the relay does not take place. 

The number of objects N available for selection is 
determined by the code formula: 


N=n(n~1), 


where n is the number of discrete resonance frequencies 
into which the available frequency band can be subdi- 
vided for use as remote control channels without any sub- 
stantial interference occurring between adjoining fre- 
quency channels. 
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Fig. 1. Circuit diagram for control-ob- 
ject selection without local power sup- 
ply, using a dual-frequencies code. 
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Fig. 2, Object-selection circuit by codes of 
dual frequencies with amplifier. 


Thus, for instance, when the frequency band avail- 
able is 300 to 3000 cps, inductances, having cores of 
a silicon-iron-aluminum alloy (similar to the one known 
in the U.S.A. under the trade name of “Sendust) are 
used, n=10 to 20, 

This scheme can also be used when the control sig- 
nals are transmitted by radio. The above code formula 
determines the number of selections available, all oper- 
ating on a single carrier frequency; the receivers may be 
different. 

The time t needed to charge condenser C to a po- 
tential equal to 90% of the ultimate value is obtained 
when R, > Z, and Z.” Ry [1], by means of the formula; 


t => ont -CZy, 


where Zy, = Quy is the equivalent total resistance of the 


circuit, ny—~the transformation ratio of the autotransform- 
er, R a7 the resistance of the diode, and Rs~the emitter- 
collector resistance of the transistor in the nonconductive 
state, Parameters are chosen so that t is equal to 0.3 to 
1,5 sec, 

The discharge time constant T of the condenser C 
through the transistor while the relay is being operated is 


ta C(Rp + Rx), 


where R. is the resistance of the relay and R,,—the emit- 
ter-collector resistance of the transistor (in the high con- 
ductivity state), 
The energy needed to actuate the relay is 
72 
W > r- te, 
P 


where E. is the voltage needed for causing the relay to 
operate, t..—the time it takes the relay to operate. 

On the basis that the condition T=kyt where 
k,=1.5 to 2.0 is fulfilled, we have: 

kit 
C>a— > 
“ Ay t- RR, 

The potential E,. on the condenser at which the re- 
lay will operate (disregarding losses in the condenser) can 
be expressed as follows: 
yt 


Ry 


E. a E> 


The same resonant circuits LyC, and LC, can be 
utilized also for selection and control of a second object 
by means of transmitting the frequencies in inverse order 
(first f,, then f;), This decreases the average number of 
resonant circuits per control object to one; it is accom- 
plished by addition to the scheme represented in Fig. 1 
of_a.second transistor, one or more diodes, another con- 
denser C, and‘a second relay, Such a scheme was used 
in one of the control systems developed by the Institute 
for Automation and Remote Control of the Academy of 
Sciences of the USSR [IAT]. 
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In the scheme of Fig. 1, the energy for charging of 
condenser C is derived from the first tuned circuit, and 
in consequence, a higher audiovoltage must be transmit- 
tedover the line. To remove this shortcoming, an ampli- 
fier is added at the output from the first tuned circuit. 
One such scheme is represented in Fig, 2, 

In distinction from Fig. 1, here the output voltage of 
frequency f, from transistor II; is applied through trans- 
former tT, to a pair of halfwave rectifiers, the output 


from which is connected to condenser C, The oscillation 
of frequency f,, transmitted from the dispatch point, only 
opens the transistor, This permits keeping the amplitudes 
of the control frequencies transmitted over the line at a 
lower level. 

The distant operational control units at the execu- 
tion points may differ from the control units that serve 
for selection of the control objects,in that they are pro- 
vided with an additional resonant circuit of frequency fy 
and another transistor, connected in a manner similar to 
that of the L,C, circuit and transistor I». 

In that case the distant control command, “switch on? 
is transmitted by a combination of frequencies f, and fp, 
but the command, “switch off* —by consecutive sending 
of frequencies f, and fs. 

It is thus obvious that damage of any component of 
the diagrams shown in Fig. 1 and Fig, 2 or opening of any 
of the circuit connections can result only in a preventa- 
tive noncompliance (rejection) but not to a false selec- 
tion or execution of a false command, 

2, Two-Section Filters 





The L-C resonant circuit as a means for frequency 
selection is distinguished by the highest simplicity; how- 
ever, it has relatively low selectivity for off-frequencies 
far from resonance, i.e., when there is large frequency 
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where Q =(wL /R) is the figure of merit of the circuit. 

To decrease the interaction between the tuned relays 
and the resonant circuits,it is necessary to make the inter- 
space between adjoining resonant frequencies several 
times larger than the bandpass width of the circuit. The 
utilization of the available frequency range then be- 
comes very ineffective. 

Consequently, it is advantageous to use resonance 
circuits in those cases where the available frequency 
band is wide enough to accommodate the entire control 
system in question,with something left to spare, This us- 
ually is the case when one is concerned with relatively 
simple control systems, where the problem of economy 
in the transmission of the control signals is not of prime 
importance. 

High efficiency of transmission-channel utilization 
is achieved by the use of multisection (iterative) band- 
filters having high attentuation in the nonpass regions. 
And the higher the requirements for attenuation in the 
nonpass regions are, the more iterated sections the band- 
pass filter must have, 

In a number of cases it is advantageous to use the 
standard four- or six-section filters of voice-frequency 
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telegraphy, which are produced commercially as compon- 
ents of systems for multiplexing of telegraph channels. 

Increase of attenuation in the nonpassing regions is 
accompanied by decrease of efficiency and increase of 
attenuation in the passing range of the filter. Therefore, 
for many remote control systems,a reasonable compro- 
mise must be found between a single-section filter and a 
complex filter having a great many inductances and ca- 
pacitances, 

For the remote control of disperse objects,a good 
solution of this problem is provided by the two-section 
filters which are widely used in radio circuits; they make 
it possible to increase the number of frequency relays in 
a given frequency band several times with respect to the 
number of relays feasible when resonant circuits are used. 
When a relative attenuation between adjoining operation- 
al frequencies of the order of 2 nepers (which is equiva- 
lent to a factor of about 7) is permissible, such filters are 
in no way inferior to more complex filters. 

Two-section filters can be installed to replace re- 
sonantL-C circuits without substantial change in the op- 
erations of the system, and their design and adjustment 
are distinguished by simplicity [7]. 

The circuit diagram of a two-section filter adapted 
to replace a resonant circuit arrangement is shown in 
Fig. 3. When, in this diagram, Cy=C,=C, the coeffi- 
cient of coupling k is determined by the formula: 


Cc 
k=—e. 
where C. is the coupling capacitance. 

It is known that when the coefficient of coupling k 
is larger than(1//Q,0,)the amplitude versus frequency- 
response curve of the filter has two humps, and that the 
width of the pass band increases with k*0Q,Q., where 
Q, =(WL4/R ) and Q,=(wL,/R2) are the figures of merit 
of the first and second sections, 

In Fig. 4 two curves are presented: one showing the 
power P, in the second section, and the second the ef- 
ficiency , of the entire filter (exclusive of the associ- 


ated tuned relay); both curves are plotted versus the 
quantities Q,Q,k* as abscissas. A maximum of power is 
transferred to the second section when k*Q,Q, = 1. 

In Fig. 5 are presented generalized amplitude -fre- 
quency characteristics of such filters when the two sec- 
tions have components of equal value (L,=L2, Cy=C2, 
Qi=Q2=Q). In this figure the abscissa = Q[(f/fy)—(fp/f)) 
represents a quantity which is characteristic of the “off- 
tune” condition with respect to the resonance frequency 
fy, © is the ratio of the output voltages of the filter or of 
the circuit at resonance to the voltages in the “off-tune” 
condition, and 6 = kQ, 

In Fig. 5 there is also shown a generalized charac- 
teristic curve for the tuned-to-resonance circuit. It is 
evident from Fig. 5 that a decrease of 6 increases the 
selectivity of the filter when the off-tune frequency is 
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far away from the resonant frequency; it is therefore ad- 
vantageous to choose for the tuned relays 6< 1, 

Let us compare the properties of the resonant circuit 
scheme and the two-section filter scheme with respect to 
their use for selection and operation of control objects. 
We assume that technical requirements make it neces- 
sary to have at the adjoining frequency an attenuation 
which is at least 10 times that of the frequency in opera- 
tion. The horizontal line o = 10 on Fig. 5 corresponds to 
this condition, For o =10 the resonant circuit scheme 
requires an off-tune value of a=10, while the two-sec- 
tion filter secures the same 9 = 10 at &=3,6 and 8 =0.5; 
i.e., with the two-section filter about 1/3 of the inter- 
space between adjoining operational frequencies is need- 
ed, 

It should be noted that the attenuation of a two-sec- 
tion filter increases much more rapidly with the distance 
from the resonance frequency than that of the tuned cir- 
cuit, It the attenuation of the tuned circuit increases in 
proportion to the first power of the off-tune characteris- 
tic &, then the attenuation of the two-section filter in- 
creases nearlyin proportion to the second power of this 
quantity. 

When 8 is decreased below 0.5, the qualifications of 
the two-section filter for selection of control objects im~- 
prove only insignificantly, while the transmitted power 
and the efficiency decrease sharply. 

In Fig. 6 are shown three curves expressing the ratio 
(/'g) as a function of o for various values of B. Here 


Oy is the off-tune value of the tuned circuit, and Oe, the 


off-tune value of the filter. The curves of Fig. 6 show 
how many times greater the interspaces between adjoin- 
ing operational frequencies must be when tuned circuits 
are used than when filters are used, to obtain equal dif- 
ferences of attenuation 0 between adjoining operational 
frequencies, 

The load can be connected into the second section 
of the filter in the same way as into the output of a sin- 
gle-section tuned circuit, and the design computations 
are also the same. The same is true for the input connec- 
tions. 

To obtain higher sensitivity it is necessary to employ, 
in both the filter and tuned-to-resonance circuits, induc- 
tances with high figures of merit Q=(wWL/R). At present, 
in the audiofrequencies range, these requirements are 
met by inductances having a core of an aluminum sili- 
con-iron alloy ("Sendust") or a ferrite core (such as 
the one known under the trade name of “Ferrocart"), The 
former are especially advantageous when relatively high 
power has to be handled in the output from the tuned cir- 
cuit or filter (as, for instance, when a neutral relay is 
used), 

Inductances with ferrite cores have smaller over-all 
dimensions and permit simpler coil windings; therefore, 
they may have relatively wider applications, The neu- 








tral relay in the output from a tuned circuit or bandpass 
filter is usually connected by way of a transistor ampli- 
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The results of the study of the noise in an 0.4/6 kv oilfield powder network at 150 cps to 100 ke are given. 
Experimental data is analyzed statistically and approximate formulas for the functions of the pulse noise am- 


plitude distribution found. 


Side by side with the utilization of high-tension 
power lines of hundreds of kilovolts for communication 
and remote-control purposes ,greater attention has been 
paid in recent years to low-voltage power networks of 
0.4 - 35 kv. The interest they have aroused is due first, 
to the great extent and branching of these networks and, 
secondly,to their sufficiently high reliability. The use 
of power lines for communication and remote-control 
purposes provides in many cases considerable saving in 
equipment, owing to the elimination of separate commu- 
nication channels, 

The superposition of remote-control channels over 
power lines meets in practice with many difficulties, 
which are connected with the simultaneous transmission 
over the line of a commercial frequency voltage, and 
with a large attenuation as compared with normal trans- 
mission channels. The high attenuation of the power 
lines is due to nonhomegeneity and branchings,as well 
as to the large number of loads which shunt the line. 

It is only natural that under such conditions the adoption 
of these lines for high-quality communication chanaels 
becomes too expensive and complicated. Moreover,the 
difficulty of installing highly reliable remote -control 
systems over these lines is due to their high noise level, 
produced by switching and various electrical devices. 

Recently,several organizations studied the para - 
meters of 0.4-35 kv networks and found that it was pos- 
sible to install communication and remote-control chan- 
nels over such lines. Work in this sphere has been carried 
out in mining and agriculture. The main attention in 
these investigations was paid to the measurement of at- 
tenuation and input impedances of lines , and the least 
attention was devoted to noise, whose levels were esti - 
mated approximately. The knowledge of the nature and 
level of noise, however, is one of the conditions for a 
correct design of a communications system and,especial - 
ly, of a remote-control system. 

The constantly growing interest in using power net- 
works for remote-control channels and the absence of 


sufficient data on the noise in these lines, made it neces- 
sary to carry out a fuller investigation of the nature of 
noise in power networks, In the present work we give 
some of the results of such investigations carried out on 
a 0.4/6 kv oilfield power network. 


1. Noise in an Oilfield Power Network 





A 0.4/6 kv power network of one of the oilfields in 
the Groznyi area was chosen as the subject for investiga- 
tions. The construction of the network is sufficiently 
typical for a system with distributed objects. The 6 kv 
open wire (sometimes cable) lines connect a 35/6 kv 
step-down substation to distributing transformers. In ad- 
dition to the 6/35 and 6/0.4 kv transformers, powerful 
induction motors for compressors and boring installations 
can also be connected to the 6 kv line. Each transformer 
point has several radially distributed 380 v lines for sup- 
plying groups of 15-30 kw induction motors of the rocker 
platforms, Welding equipment is also fed from the low- 
voltage network. The 6 kv lines are usually several km 
long, and the 380 v lines,several hundred meters. 

Both the high-and low-voltage lines work in a star con- 
nection with an ungrounded neutral line. The noise in 
these lines was studied ona phase-to-phase basis. 

The above construction of a power network is also 
typical for other industrial establishments.such as mines, 
large factories, chemical combines, etc. 

The investigation of noise in power networks, like 
that of any other communication or remote-control 
channel, consists in the first place in determining the 
nature of the noise,and then measuring its frequency and 
amplitude characteristics, Both the qualitative and quai- 
titative noise characteristics must be studied over a wide 


" range of frequencies for several fixed bandwidths of the 


receiving devices,normally with a simultaneous measure- 
ment of the attenuation and the input impedance of the 

lines. 
The noise type, which later determines the methods 
of measuring the noise amplitude and frequency charac- 
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Fig. 1. Harmonic components of 50 cps. 


teristics, can be conveniently determined from an oscil- 
logram of the noise. The recording of the noise can, for 
instance, be made with a loop oscilloscope whose vibra - 
tors are connected to the output of a selective device. 
When operating with a multiloop oscilloscope,it becomes 
possible to record noise from several receivers with dif- 
ferent bandwidths and at different frequencies. This is 
important for comparing noise at various frequencies for 
a transient condition in the line. 

Noise was measured in the range of 150 cps to 100 
kc. A preliminary recording of the noise on an eight- 
loop oscilloscope type MPO-2 revealed several bands 
with noise of different types. In the lower frequency 
part of the range (up to 5-10 kc), the high noise level is 
due to the harmonics of the commercial frequency of 
50 cps. The level) of this noise can be measured in ef- 
fective values directly at the output of narrow-band de- 
vices. Fig. 1 shows a discrete noise spectrum in this 
range for a 6 kv (a) and a 380 v (b) network. The level 
of harmonics was measured by means of a selective 
voltmeter RVS-253, It will be seen from the drawing 
that, owing to the three-phase construction of the line, 


only odd harmonics are present in the network (n is the 
number of the harmonic). The relative level] of harmon- 
ics with respect to the line voltage is lower in the high- 
than in the low-voltage network. Repeated measure - 
ments of the harmonic voltages at different times of the 
day,and in various parts of the network,showed that the 
level of separate harmonics can vary by 20-50%, but 
that the nature of their distribution remains unchanged. 
The study of the oscillograms of separate harmonics 
showed that,with the rise in the order of the harmonic 
and an increase in the bandwidth of the measuring de- 
vice, the sinusoidal shape of the harmonic becomes dis- 
torted,owing to the pulsation of adjacent components 
(Fig. 1c). 

If the measuring frequency is raised above 5 kc,and 
the bandwidth increased, the received noise acquires a 
fluctuating nature. The amplitude probability density 
distribution of this noise approaches the normal law (in 
the bandwidth of the order of several kilocycles). This 
range extends from 5 to 30 kc. It should be noted, how- 
ever, that,from 15 kc up,in addition to the fluctuation 
noise there appear occasional noise pulses, whose over- 
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shoots exceed considerably the leve! of the fluctuation 
noise, The noise levei in this range can be measured 
directly by means of a normal tube voltmeter type LV-9 
connected to the output of a selecting device. It will 
then read the root -mean-square value of the noise de - 
creased by 10%. Obviously,in this measurement the 
level of occasional noise pulses is not taken into account. 
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Fig. 2 shows the curves of the rms values of the fluctua - 
tion noise in the range of 5-15 kc for the 380 vy network 
(dotted line) and for the 6 kv network (full line) at a 
bandwidth of 2 ke. 

From 30-35 ke up to several hundred kc, with a 
sufficiently wide bandwidth (about 10 kc), the noise 
acquires a pulse nature both in the high-and low-volt- 
age networks (Fig. 3a and b). Obviously.the boundary 
between the two ranges is nominal. 

The study of the oscillograms of pulse noise in this 
range led to the establishment of a synchronous relation- 
ship between the appearance of the overshoots and the 
voltage of the network. The noise pulses have a repeti- 
tion frequency of 100 cps. Owing to the three-phase con- 
struction of the line the pulses which correspond to the 
other two phases appear attenuated to 1/2 or 1/3 of the 
value of the main pulses, At much less frequent inter- 
vals, single pulses or groups of pulses appear, which are 
not synchronized with the frequency of the network. 
They are due to switching phenomena in the line. Tests 
have shown that such a noise structure is peculiar not 
only to oilfield networks, but also to networks of mines 
and industrial establishments, i.e.,to networks with a 
large number of electrical installations, which are simul - 
taneously connected to the mains. 

The oscillographic investigation of the noise has 
also shown the nature of the noise due to separate 
sources, and the simultaneous recording of calibrating 
signals provided a scale for their measurement. 

The main source of noise in the 380 v network is 
due to welding and therefore the noise due to the opera - 
tion of welding equipment was investigated with partic- 
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ular care. The study of the oscillograms revealed that 
the nature of the noise during welding remains the same 
as before the welding, but its amplitude increases by a 
factor of 20-30. 

Fig. 3a shows an oscillogram of noise in a 380 v 
network at 45 kc, The gain of the recciver was reduced 
to 1/10 when recording noise during welding. The top 
oscillogram shows the noise recorded directly at the out- 
put of the intermediate frequency amplifier with a band- 
width of A f = 200 cps, and the bottom one shows the 
recording of the noise envelope at the output of the amp- 
lifier with a bandwidth of Af = 2 kc. The amplitudes of 
separate noise overshoots at a bandwidth of about 10 kc 
in the 30-60 kc range amount to fractions of a volt dur- 
ing welding. The noise spectrum during welding extends 
to hundreds of kc, but the level of noises drops rapidly 
with frequency, owing to a sharp rise in attenuation. 

High level of pulse noise in the 6 kv and the 380 v 
networks is also caused by the switching-in of heavy 
500 kw motors, It will be seen from the drawing that. 
during the starting of the motor.a characteristic sequence 
of 100 cps pulses appears. The switching-off of the same 
motor leads to a short splash of noise of considerable 
amplitude. Value of the overshoots during the transient 
period of motor operation is 20-50 times larger than 
during its steady-state running, and amounts to several 
tenths of a volt at frequencies between 30 and 50 kc in 
a bandwidth of about 10 kc. 

The recording of noise during the switching -in of 
the motor was carried out simultaneously on two fre - 
quencies. The top oscillogram in Fig. 3c shows noise 
recorded at 60 ke at the output of a straight amplifier 
with a bandwidth of 10 kc. The middle oscillogram rep- 
resents the recording of noise at 120 ke directly at the 
output of the intermediate frequency amplifier with a 
bandwidth of 200 cps and,finally,the bottom one repre - 
sents the recording of the envelope of this noise at the 
output of the amplifier with a bandwidth of 2 kc. 

It is convenient to estimate the value of noise pul- 
ses by the number of pulses which exceed a given level 
per unit of time. By counting the number of pulses 
which exceed various levels,it is possible to plot the 
amplitude characteristic of pulse noise in the network 
at a given frequency. For this purpose,a stepped device 
with a controllable level of operation is connected to 
the output of the constant bandwidth selective receiver. 
The counting of pulses which exceed a given level can 
be carried out either by means of compact counters 
working on cold-cathode thyratrons,or by measuring the 
mean repetition frequency of the pulses. The amplitude 
characteristic can also be determined from the oscillo- 
grams of the pulse noise. 

The frequency characteristic of the noise level 
variations can be plotted on the basis of a family of 
amplitude characteristics whose parameter is the fre - 
quency. Fig. 4 shows curves of the variations of pulse 
noise level with respect to frequency for 380 v networks 
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Fig. 3. Oscillograms of pulse noise: a) in the 380 v network, b) in the 6 kv network, c) 
when connecting and disconnecting a motor in the 6 kvnetwork. 
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(dotted line) and for 6 kv networks (full line). The curve 
parameter consists of the number of pulses which exceed 
a given level in a unit of time. The bumps in the curve 
are due to resonance phenomena in the line. For other 
lines, resonances may occur at other frequencies; how - 
ever, the characteristic of all these curves consists in 

the reduction of the noise level with a rising frequency. 


2. Statistical Characteristic of Pulse 
Noise 





It is known that a random process such as the pulse 
noise is completely described by a multidimensional 
distribution law. .In practice, however, it is possible to 
use simpler statistical characteristics. In the first in- 
stance,it is possible to find the distribution of the over- 
shoot noise amplitudes, which is fairly fully represented 
by a single-dimensional function of the amplitude prob- 
ability density distribution. This characteristic also 
provides the most objective means for comparing the re- 
sults of various noise level measurements in the channel, 
and makes it possible to evaluate their interference for 
various transmission methods of the signal. 

The problem of the statistical analysis of the 
measurement results amounts to finding the type of dis- 
tribution law which corresponds to the data obtained. 
The most simple and best known is the norma! distribu- 
tion which is completely determined by two numerical 
parameters: the mean value and the root-mean-square 
deviation. Hence,in the case of asymmetrical distribu- 
tions,one aims at finding a function f(x) which will 
transform the given distribution into a normal one. 

This means that,irrespective of the type of distribution 
of x the distribution of the converting function f(x) is 
normal. Such a conversion is of practical interest,since 
the distribution of the converting function f (x) is fully 
determined by the root-mean-square deviation o and 
the mean value f(x). 

The amplitude characteristics of the pulse noise in 
the 380 v and 6 kv networks, and the data obtained from 
the statistical analysis of the oscillograms of pulse 
noise, give reasons to believe that the distribution of 
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pulse amplitudes is represented by logarithmic-normal 
distribution laws, This means that the logarithm of the 
random value x is distributed according to the normal 
law (the converting function is f (x) = log x). The densi- 
ty of distribution has the form 
{ log: —logn)? 

pone = i x0 SE]. 
where log @ is the mean value of the distribution of 
log X. Since log 4 is the median of the logarithmic-nor- 
mal distribution, a must correspond to the median in 
the distribution of x. It can be shown that 


o 


z a-10™ 
and the abscissa,which corresponds to the maximum dis- 
tribution density value (mode), o 


Tmax=a-10 ™, 


where M = log e= 0.4343, The values of the median, 
the mode,and mean distribution of x are related to 
the parameters of the transformed distribution by the 
formulas 


log T max==log a —2.3 3? (2) 
log To,5=log a, (3) 
log x =loga +1.150?. . (4) 


The difficulty in determining the theoretical dis- 
tribution parameters from the amplitude characteristics 
of pulse noise is due to the presence in the channel of 
a fluctuation noise level, which is small as compared 
with the pulse noise, and the existence of internal re- 
ceiving device noises. This makes it impossible to count 
at zero level the number of overshoots which are only 
due to pulse noise. Nevertheless, it is possible to assume, 
owing to the difference in the levels of fluctuation and 
pulse noise (the difference increases with rising frequen- 
cy), that a large part of the amplitude characteristic in 
the region of large amplitudes is entirely cue to pulse 
noise. The distribution in that part of the amplitude 
characteristic after logarithmic transformation can be 
regarded as a truncated normal distribution with an un- 
known cut-off level x» (it is unknown how many pulses 
there are below the xq level). The parameters of the 
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full distribution can be determined either from the tables 
or from three arbitrary points of the empirical probabil - 
ity density distribution curve of the number of over- 
shoots N*(u) (Fig. 5) taken from the formulas 
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where uj = log xj, N* (uj) is the value of the probability 
density at point uj, and 
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Fig. 6. Diagram of the distribution of pulse noise 
amplitudes, 1) 380 vnetwork, f = 35 kc (i = 0.75, 
0 = 0.35); 2) 380v network,f = 50 kc (u = 0.55, 

oO = 0.26); 3) 380v network, f = 70 kc (i = 1.31, 

o = 0.255); 4) 6 kv cable network, f = 115 kc 

(u = 0.03, o = 0.17). 
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It will be seen from formulas (5) - (7) that the de- 
termination of the parameters of the converted distribu- 
tion is obtained by means of the ratio of N* (uj) values; 
hence, it is not necessary to normalize the curve N* (u), 
and the scale along the y axis can be arbitrary. 

The norma! distribution law was checked by plot- 
ting the distribution function points on special probabil- 
ity paper (Fig. 6) on which a norma! distribution curve 
is represented by a straight line. The method of 
straightened curves[1] provides a fairly accurate indica - 
tion of even the smallest deviation from the normal dis- 
tribution law. Some 20 amplitude characteristics of 
pulse noise were investigated in 380 v and 6 kv power 
networks in the range above 30 kc, In all the cases,a 
good approximation to the logarithmic-normal law was 
obtained for the range of the large amplitude overshoots. 
Fig. 6 shows several straightened -out distribution dia - 
grams. The deviation from straight lines in the region 
of small amplitudes is due to the appearance of over- 
shoots caused by fluctuation noise. 

The above research makes it possible to arrive at 
the conclusion that the amplitude characteristics of 
pulse noise in 380 v and 6 kv power networks are correct- 
ly approximated by a logarithmic-normai distribution 
law, for whose complete definition it is sufficient to 
know only two numerical parameters, The knowledge 
of the distribution parameters simplifies considerably 
the measurement of noise and the comparison and de- 
sign of noiseproof communication and remote-control 
systems for power networks. 

The author wishes to express his gratitude to G. A. 
Shastova for her continual interest in the work,and to 
V. A. Zakharov for his assistance in testing. 
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The application of the dynamic characteristics of ferrite cores which have rectangular hysteresis loops, for 
designing electromagnetic circuits,is described. A formula for the transient process in an electromagnetic circuit 
of known parameters is given. The parameters of a complex circuit with ferrite cores,which provide the required 
working conditions,are determined. Experimental results obtained in checking the above technique are cited. 


1, Introduction 

Modern high-speed devices employed in automation 
and computers use circuits which contain ferrite cores with) 
rectangular hysteresis loops on awide scale. The cores 
of these devices work under conditions of dynamic re- 
magnetization which take microseconds or fractions of 
a microsecond. Confrary to the static conditions for 
which the behavior of a ferromagnetic core can be rep- 
resented by the characteristic B = B(H), there is no simi- 
lar characteristic for dynamic conditions, if the law of 





variation of the magnetic field H(t) or B(t) is not known. 


The determination of the dynamic characteristic 
which represents the behavior of a ferromagnetic core 
under the effect of an arbitrary pulse is the main prob- 
lem for any computation of the above electromagnetic 
circuits. 

Physical processes occurring in dynamic remagnet- 
ization can be accounted for by means of generalized 
electrodynamic factors such as eddy currents and mag- 
netic viscosity. In ferrite cores magnetic viscosity is 
the determining factor,since the low conductivity of the 
core materials practically eliminates any eddy currents. 

The laboratory of the Theoretical Foundations of 
Electrical Technology Department of the Moscow Elec- 
trotechnical Institute derived-on the basis of experimen- 
tal, and at a later date, theoretical investigations of dy- 
namic remagnetization of ferrite cores with rectangular 
hysteresis loops-an equation which represents the be - 
havior of ferrite when it is arbitrarily remagnetized in 
{1-3}: 


dB 7 
a =z (28. — 


2B — B, — B, \2 
[1— laa ) Ju — Hs), 
where B, By, By, and B, are the magnetic inductions at 
instant t, respectively, the initial, residual and saturated 


B, — By) 
(1) 





inductions; H is the strength of the external field at in- 
stant t; Hg, is the field strength of B,; = B on the static 
hysteresis loop which corresponds to the initial magnetic 
condition of the ferrite; 5 is the coefficient of magnetic 
viscosity which is determined experimentally. 

If we denote the position of the center of a particu- 
lar half-cycle as 


i 
Be = > (Br + Bu), (2) 


equation (1) can be represented in a differential form as 
follows: 


% = 3(B,— Be) 


7 (3) 


[1 ( 


or in an integral form 


BB, \2 
Fe) | (4 — Ha) =r (B)(H — Hy) 


t 
Q = Q(B) =\ (H — Hx) at = 





(4) 
| —— Be ~1Pu — 8 \ 
= = (tanh B,—B,~ ah pp.) - 


The left-hand side of the latter expression by anal- 
ogy with mechanical terminology can be called the in- 
tensity pulse of the effective field. Under the intensity 
of the effective field.we understand the excess of the ex- 
ternal field strengh H over the field strength Hst on the 
corresponding static hysteresis loop at B = B,,. 

Equations (3) and (4) represent the dynamic proper - 
ties of the ferrite cores and confirm sufficiently well the 
experimental characteristics of remagnetizations of the 
ferrite cores, starting with pulse fields which exceed 
the coercive force by a factor of 2[3]. The value of 
5 at room temperature is approximately equal for differ- 
ent ferrites under both complete and partial remagneti- 
zation conditions: 








6 = (6 — 7) 10° cm/coul 


Of the many problems encountered in the calcula- 
tions of electromagnetic circuits ,we shall only discuss 
two: the calculation of the transient process in the cir- 
cuit with known parameters and the determination of 
the parameters of a complex circuit with ferrite cores, 
which provides the required condition of operation. 

The solution of these problems is made by means of 
the above dynamic characteristics of a ferrite core. 


2. Calculation of a Transient Process in 
a Circuit Which Contains a Toroidal 
Ferrite Core 

Let us first determine the voltage pulse across the 
winding of a core in the most simple case, when the 
core has one winding with a number of turns w, in 
which current i, (t) is flowing. The geometrical dimen- 
sions of the core are also known (and hence the length of 
the mean line of force 7 and the cross section S) as well 
as its magnetic properties ( Ho, B,, Bs) and its initial 
magnetic condition By;. 

In the absence of skin effect and stray fields, the 
magnetic field strength and induction in a ferrite core 
can be represented as mean values: 











Li 

H = Hy =—*» (5) 
t 
| uyae 

B= Bn= 5 + Be. (6) 


Let us replace,in our calculations the actual static 
hysteresis loop of the core by a rectangular one with a 
coercive force 


_~ 


Hy =~ \H, dt, (1) 
0 


where rt is the time of nominal remagnetization. 
From expression (4) it follows directly that 


B—Be -1 3, — 8, 
B.— Be = th |8@«) + tanh rae | ° (8) 


By substituting (8) in equation (3) we obtain an ex- 


pression for the voltage pulse across the winding of the 
core: 


St (BL, — B.) (H(t) — FE 
#1 () = 058 8) Or me 


2 -1 5y— Be 
ch [se () + tanh ed 
The field strength H (t) in the core is completely 
determined by the current i; (t) under the above condi - 
tions, and therefore the voltage pulse can be calculated. 
Let us now examine the case when the ferrite core 
N has another winding with w, turns loaded by an effec- 
tive resistance R, (Fig. 1a). Now the strength of the 
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magnetic field H(t) in the core will be determined by 
an algebraic sum of the field strengths established by 
the known currents i; (t) in winding w, and the unknown 
current i, ft) in winding w. 

In order to calculate the voltage pulses u, (t) and 
u, (t) across both the windings of the core, let us first 
refer the parameters of the electromagnetic circuit to 
the material of the core. Such a change-over involves 
the following recalculation of values: the voltages are 
multiplied by a coefficient of the type V/ w,§, the cur- 
rents are multiplied by a coefficient w,/ land the re- 
sistances by 1/w%,S. The dynamic condition equa- 
tion (3) makes it possible to represent the ferrite core 
in the form of an equivalent circuit, in which a nonlin- 
ear resistance r(B) is connected in parallel with the 
source of the static field strength Hy. The complete 
equivalent circuit of the above core is shown in Pig. 1b. 

Let us replace the parallel connected resistances 
r (B) and r, by a single resistance R(B) and let us calcu- 
late the voltage pulse u,(t) for a certain equivalent 
core Ne, which has only one winding with a number of 
turns w, in which a known current i,(t) is flowing (Fig. 
1c). The dynamic characteristic of core Ne is easily 
found from 


we 
Qe(B) = Q(B) + 3-3 (B — Bu) = Q(B) + Qr29) 


If core N had several output windings with effec- 
tive resistances, the second term of the right-hand side 
of expression (10) would consist of the appropriate sum. 

Further determination of the voltage pulse u;(t) can 
be conducted more conveniently by a graphical metnod, 
since in order to obtain the analytical relation dB/dt = f(t 
it is necessary to solve the transcendental equation 
(10). We determine the value of H,(t) for known values 
of wy, 1 and iy(t) and plot the curve for 

t 


Qe (t) = \ (H, — Hy) dt = Q,(t). (11) 
0 

By a simple addition of the ordinates of the dynamic 
characteristic Q(B) of the core and of the straight line 
C,.we find the characteristic of Q.(B). From the rela- 
tionships Qe(t) and Qe (B) thus obtained,we find graphi- 
cally the relation B(t). The further calculation of the 
voltage pulse u, (t) is carried out from the formula: 

nO) = yt — H)wS, 
into which are inserted values of H, and B obtained for 
a definite value of t from the curves of the relation of 
the above values to time. 

The electromagnetic circuit with the voltage pulse 
source is transformed to the cases examined above by 
converting the source of voltage to a source of current. 
This statement is made clear by Pig. 2. 

The nonlinear nature of the effective resistance of 
the input or output circuit of the core provides further 
difficulties, which can be resolved by methods for the 








solution of nonlinear equations. Similar difficulties are 
encountered when the input (output) impedances of the 
core windings are inductive or capacitive. 

Figure 3 shows the voltage pulses u,(t) for various 
laws of the variation of the magnetizing current i, (t) 
and various initial conditions of By; of a ferrite core 
BT-2 with dimensions D x d x h = 4.85 x 2.84 x 1.44 
mm? and magnetic characteristics; B, = 2200 gauss, 

B, = 2500 gauss, 6 = 6.35: 10%cm /coul, Hy = 0.8 amp/ 
/cm, Curve 1 corresponds to H,,,= i,w;/7 = 2.8 amp/ 
/cm= const and By; = =Br; curve 2corresponds to H,,= 
= 2.8 amp/cm = const and B,, = 0; curve 3 was obtained 
for a linearly rising current i, = kyt(ky = 8.9 amp/psec 
and By = —B,; curve 4 was obtained for the effect 
(Byy = —B,) of a transcendental current pulse on the fer 
rite core, with the current rising up to the value of 

im = 0.9 amp = const according to the same linear law 
i; = kyt. 

The voltage pulses were calculated from formula 
(9) for Wy = 5, We= 3,and Rp = @. 

Figure 4 shows the oscillograms of the magnetizing 
current i(t) and voltage u,(t) for the corresponding 
windings of the same core,both for open-circuit and 
loaded conditions,with an effective resistance R,. In 
the initial condition (t = 0) the core is characterized by 
induction By; = —Br. Curves 1 and 2 were plotted for 
R, = ©; curve 3 for R, = 50 ohms; curve 4 for R, = 10 
ohms. Points show the calculation results. Figures with 
dashes refer to voltage oscillograms. 


3. Determining the Parameters of Elec- 
tro- Magnetic Circuits Which Contain 
Several Ferrite Cores 

In the controlling and logical elements of comput- 
ers, circuits with ferrite cores are used. 























Fig. 1. a) initial schematic of the 
electromagnetic circuit, b) its equiv- 
alent circuit, c) a circuit with an 
equivalent core. 
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The basic problem in designing these circuits con- 
sists in determining the winding data for the selected 
cores, the values of resistances in coupling circuits,and 
the required magnetizing pulse -which would provide 
the operation of the device for the given conditions. 

Let us show how it is possible to solve the above 
problem by means of the dynamic characteristic of the 
core. 

For a concrete example let us examine the widely 
used circuit of a shift register, constructed with trans- 
former elements (the ferrite core operates as a transform- 
er). Figure 5a shows a circuit of a two-phase shift register 
and Figure 5b and c that of a three-phase register, The param- 
eters of all the cores and coupling circuits are the 
same. The operating (clock) pulses iq, , iy, and its 
have similar values but arrive at the input windings of 
the cores at different times, establishing negative field 
strengths. In the initial position all the cores of the 
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Fig. 2, a) initial schematic of the 
electromagnetic circuit with a volt- 
age source, b) transformation of the 
voltage source into a current source 
i, = u/R, c) equivalent circuit. 
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registers are in the condition of a negative remanent in- 
duction —B,. The operation of the clock pulses does 
not affect the magnetic condition of the cores. Let us 
now assume that at the input of core 1 there appears a 
pulse which can remagnetize this core to a value of 

+B,. It is obvious that the remagnetization of the first 
core will not affect the remaining cores of the register, 
since a flow of current in the coupling loop between 
cores 1 and 2 is prevented by the rectifying element 
(diode). The next instant, under the effect of the clock 
pulse iy,,the first core will be remagnetized to its 
initial condition (—B,),and the current produced by this 
action in the coupling loop between cores 1 and 2 will 
remagnetize core 2 into condition +B,. Under the ef- 
fect of pulse it, the unit information which corresponds 
to induction +Br is shifted to the third core. In this pro- 
cess currents arise both in the coupling loop between the 
second and the third cores (forward movement of infor- 
mation) and between the second and first cores (reversed 
movement of information). In the two-phase register 


a 
bY, V 2 
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(Fig. 5a) the reversed current produces a certain change 
in the condition of the first core. 

In the three-phase register the magnetic condition 
of the first core is not altered because the clock pulse 
is applied simultaneously to cores 2 and 1 (Pig. 5b), thus 
keeping the first core in condition —B,, or because pulse 
iy, overlaps in time pulse iy, (Fig. 5c),thus keeping 
core I in condition —B, when the second core is being 
remagnetized. 

Since the condition of work of all the links (stages) 
is the same, the designing of the register can be based on 
a single link. One link consists of two cores connected 
by their coupling circuit, and its computation consists 
in analyzing the operating conditions of the circuit, i.e., 
of the forward and reversed movement of the informa - 
tion pulses. Providing all the link parameters are the 
same ,the consideration of only one condition of a cir- 
cuit consisting of three cores 1, 2 and 3 will lead to the 
same results. 

For the sake of simplicity let us neglect the differ- 
ent resistances of the diode for different currents and rep- 
resent the impedance of the coupling circuits by means 
of the effective resistances: 


r=R+#R,. (13) 


Let us first examine the operation of a two-phase 
shift-register link. Let it be necessary according to the 
operating condition of the register to transmit unit infor- 
mation from the second core to the third in time 7, and 
let us assume that at the same time a variation in the 
condition of the first core of AB = b2B, is permissible. 
The calculation is made for the condition of a simul- 
taneous remagnetization of all the three cores. For a 
complete remagnetization of the second and third cores 
by the value of the induction change * 2B,,it is neces- 
sary to apply to them a field-strength pulse equal to 















































12 














— WwW. 
wo SS ° 


Fig. 5. Schematic of magnetic shift registers consisting of transformer elements: 
a) two-phase register; b and c) three-phase registers. 
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2, = 0(2B,) + Agt = Sy + Ht, 


where S,, is the value of the effective field-strength 
pulse for complete remagnetization, given by equation 
(4) for Bry = —B, and B =+B,. The value of S,, is 
usually called the coefficient of remagnetization. 

The value of the variation of the inductance 4B in 
the first core corresponds to a pulse of 


Q (5) = q(b) Sy» + h(b) Ht. (15) 


Here q(b) is the dynamic characteristic (4) in rela- 
tive values and 


(14) 


—B,+b2B, 


h(b) = Hy, dB (16) 


{ 
523, A, \ 


r 

is the mean value of the mean strength of the static 

field for partial remagnetization (in relative values). 
The value of the clock field-strength pulse for 

time t can be represented by the expression: 


t . 
?,= ; \ i,dt = a hat 
0 
Let us denote the ratio of the number of turns in 
the coupling loop n= w,/wy,,and let us refer all the 
values which characterize the condition of the link to 
the second core. The field-strength pulse Qy which pro- 
vides in time f a variation of induction in core 3 of 
2B, is multiplied by w,/w, when it is referred to an 
equivalent circuit of the link (Fig. 6a); i.e.,it is multi - 
plied by n and the variation in induction by w,/wy, ie., 
by 1/n. The referred admittance of the coupling cir- 
cuit of cores 2 and 3 will then assume the form: 


Ss Ss 
G = gui —- = gus—-n’. (18) 


The field-strength pulse Q(b) in the preceding coup- 
ling loop is multiplied by 1/n when it is referred to the 
equivalent circuit, and the induction variation b 2Br 
occurring in the first core is multiplied by n. The ad- 
mittance of this coupling loop in the equivalent circuit 
is equal to 


(17) 


. S 
G=guiz- (19) 


From the equivalent circuit we derive the following 
obvious expressions on the basis of Kirchhoff's laws 


2B, 
nQ,, = G'n?® 9 (n — 1), (20) 
20 = G’2B,(1 — bn), (21) 
b 
2, = 2, + nQ,, + 20). (22) 


Combining expressions (20) and (22) we obtain an 
equation which determines the ratio n of the number of 
turns as a function of b; 


n— 1 (1-2) n_ 190 _ 0, 


b maa « ee) (23) 
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Expression (22) assumes after transformation the 

form: 
Q, n—b 
a= 1" oe 

Thus the controlling factor in determining the wind- 
ing data of cores and the magnetizing clock pulse for 
two-phase registers with one diode is the variation of 
the induction due to the reverse movementof the informa- 
tion pulse. 

For the above conditions we find the value: 

Qu! 
bY = Ben 1) ° 

The value of one of the quantities, for instance 
W2fepends on the value adopted for the other quantity, 
g- Let us assume that the impedance of the coupling 
loop is known. Under these conditions we find w, from 
(25). Now.it remains to determine the value of the 
building-out resistance R which it is necessary to connect 
in series with the selected diode. Experimental investi - 
gations of diodes show that their direct resistance under 
pulse conditions can be represented by an effective re- 
sistance R, Which in practice has single-value relations 
to the current pulse flowing through the diode over a 
given time interval [4\,i.e., 


(25) 


Re= f (\ idt) = 7 (q). (26) 


0 
Since, by taking expressions (23) and (29) into 
account, the value of the current pulse q can be assumed 
to be known 


q = gw,S2B,(n — 1), (27) 


it is possible to find the value of R, from the character- 
istic of R, (q), and then find the value of the building- 
out resistor R in the coupling circuit. 

If the difference in the equivalent resistances of the 
diodes is taken into account ,the value of n as determined 
fram equation (23) is decreased. Other factors, however, 
such as the simultaneous remagnetization of the cores 
increase n, This hoids sufficiently well for the :econd 
and third cores, but the rise of the magnetic induc- 
tion in the first core is completed before the other two. 

The opposite effects of the above two phenomena 
make it possible to neglect them in the first approxima- 
tion and use formulas (23) - (25). 

In the separate link of a three-phase register, 
whose equivalent circuit is shown in Fig. 6b, contrary to 
the case examined above, the variation of induction in 
the first core is zero and the value of pulse Q is unknown. 

We obtain the following design expression from the 
equivalent circuit: 


(28) 


Qw! 
(29) 








The ratio Q7/Qy = f(n) has a minimum at n = 2. 
Hence expressions (28) and (29) take the form: 


Q; — AQ ’ (30) 
Q,,! 
ews = SoB, (31) 


Thus the parameteis of a three-phase register cir- 
cuit are calculated from the minimum of the field 
strength pulse for time T. 

Expression (30) determines the full value of the 
magnetizing pulse only in the first circuit (Fig. 5b). In 
the second circuit (Fig. 5c) this expression determines 
only that part of the pulse which is produced by 
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Fig. 6. Equivalent circuit for calcu- 
lating parameters of a register link: 

a) for a two-phase register, b) for a 
three -phase register. 


the generator during time tr. The full value of the sym- 
metrical magnetizing pulse, whose duration is twice as 
great as T, is equal to8 Quy. 

Let us note that in the coupling loops during time 
T for n = 2 there flow equal current pulses, and hence 
the equivalent resistances of the diodes can also be con- 
sidered equal. 

Methods of integral calculations are also applicable 
to other circuits with ferrite cores, which perform differ- 
ent functions in digital computers [5, 6]. 


3, Conclusions 

It follows from the above that by means of equa- 
tions (3) and (4) it is possible to evaluate circuits with 
ferrite cores for different pulse effect. 

The oscillograms of the pulse voltages (Fig. 4) ap- 
pended to the article and the experimental checking of 
model ferrite -diode circuits designed on the basis of the 
theoretical considerations showed that the computations 
were sufficiently accurate. As compared with the experi - 
mental data,the calculation error does not exceed 10 - 


20%, 
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The relations for the optimum number of pole pairs of electromagnetic slipping clutches with massive 
steel armatures are found. The change in permeability is taken into account. A method for choosing the main 
dimensions of the clutch is described. 


Electromagnetic slipping clutches (e.s.c.) and brakes The results of these formulas do not depend on the choice 
with massive steel armatures find ever wider use in drives of 4, and are uniquely determined by the given technical 


for automatically controlling speed [1] and torque [2], conditions (the torque M, and the relative speed ng). 

to provide smooth engagement and braking, and for ob- This article is a further development of works published 
taining stable slipping speeds in drives with asynchron- on the problems of designing e.s.c. and brakes with mas- 
ous motors [3], as loading devices for bench tests, etc. sive steel armatures [4-7]. 


Among their advantages are simplicity of construction, 
low consumption of nonferrous metal, low cost, ease of 
control, and the ability to dissipate a large amount of 
heat in the armature. 

The methods which have been developed for deter- 
mining the main dimensions D, 2, and the optimum 
number of pole pairs, po, of an e.s.c. with a massive 
steel armature [4, 5],are based on relations derived on 
the assumption that the permeability of the armature 


1, Optimum Number of Pole Pairs 

The expressions for the maximum torque M, and 
critical speed n,;, of an e.s.c. with a massive steel arm - 
ature and poles of alternating polarity, and also the 
equation for the mechanical characteristic, considering 
the change in», are derived in [7): 





t 
(i:17y)* pap DI,"sin (= 5) 














material, 1, does not change. In these methods, there x 55,2 Ragksd (1+ %)(1 Fin 9) ' @) 
is a certain arbitrariness in choosing » , and this leads to aad 7 
different results in calculations (see example below). , ~~ ggg th -_ np? 
; +x)? kK 

The design formulas for determining py, D and | 5 ny = [oe] fv - = ’ 
presented in this paper take into account the change in pik! ( Dkig) * D* sin® (= 3 
permeability as the magnetic field strength changes. . 
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Here M, is in kg-m; ny. —rpm; ijw, are the am- magnetization curve [8]; k, and ky are coefficients 
pere-turns of the excitation winding, in amps; p is the which depend on the exponent f [6], and can be set 
number of pole pairs; D, and1g are the diameter and equal to k, = 0,845, and’k, = 0.87 for ferromagnetic 
length of the machine along the air gap, m3; 6 — the materials, taking hysteresis into account; « is a coef- 


air gap , m;t/r—ratio of the designedpole width to pole ficient equal to the ratio of the drop in mmf in the steel 
spacing; kaq is the reaction coefficient of the armature inductor, F,, to the mmf in the air gap; F 5; 
along the longitudinal axis; ¢ — the desigi value of 


armature material resistivity, considering the edge ef- sing = ] fati- 
fect ohm-m; K is a coefficient in the equation of a dn + 1 


1 
parabola B,, = KH," » which approximates the basic In the given relations, one should assume that 
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t/*=1, kea = ; , iw, = i F,, 


for e.s.c. with nondistinct poles, where F, is the mmf 

of the excitation winding, per pair of poles. The rela- 
tions given in [5] follow from (1), (2), and (3) in the par- 
ticular case when n = 1 [ p = const }. 

The numerical coefficient in the formula for My, 
as further analysis will show, should be approximately 
53.6; this takes into account the loss-in the armature 
from the higher flux harmonics in the air gap. 

As the relative speed is increased from zero to Ng 
the reaction mmf of the armature, F2, grows, and Fs de- 
creases, It can be shown, that in the extreme 


Fs(1+%) = Fokoa = cn 


Vilitees > ag 0,56 i,w,. (4) 


The change in the mmf's Fs and F, as a function of 
speed obeys the following rules: 
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Here, Bs is the flux density in the air gap of the 
e.8.c., which results after the field in the gap is brought 
to a pole generated by the excitation winding, webers/ 
m?; Bs max is the maximum flux density in the gap for 
a given iyw, corresponding to n, = 0, 

Curves of ) and 6 as a function of M/My, for the 
stable region of the mechanical characteristic, where 


O<ns<N, are sliown in Fig. 1. 


The flux density under a pole, B_s, is approximate- 
ly twice as much as the variable component of the flux 
density in the gap, Bg, for induction clutches , since 


B., = 2Bsky , where kp = 1.05 to 1.1 [4], Moreover, 


<a 1,andk,4 0,785. for the clutches mentioned, 


Taking these relations into account, and doing a 
calculation analogous to that done in [7], we can gui 
the expressions: 


M (i)? popDI 








k= ——, (7) 
55,2-0, 785-zkgk, (1 + *) (1 + 8ing) 
2-1 anti : 
ne < [ ats see sg) | 11589(1 + »)* p ™ KIKCp (2kq) 
iw, e+! 
uzk! (0,785D) ™ D 
(8) 


Hence, the maximum torque of an induction clutch, 
for an identical number of excitation ampere turns, is 
approximately half as much as for clutches with poles 
of alternate polarity, and the critical speed is four times 
as great. It should be mentioned that the above-men- 
tioned facts are justified for nonsaturated machines 
with the same «x. In iron-clad inductors, x can be 
greater than in induction clutches,as a result of the 
large leakage flux, and this, of course, somewhat alters 
the quantitative relations obtained, as they pertain to 
induction clutches. 

The effect of p on the mechanical characteristics 
of the clutch is seen in Fig. 2, where the mechanical 
characteristics for a nonsaturated inductor (K = const ) 
are plotted in relative units (the values of My, and 
Nsky When p = p, are taken as references) for different 
values of p but with the other parameters remaining the 
same, and the excitation remaining unchanged, Figure 2 
shows that.as p decreases in the region of low slip, the 
torque increases, but becomes less when ng is large. 

Thus, if, for example, the given operating speed 
lies to the left of point A (Fig. 2.), the value p = p,/1.5 
is preferred over p = p,. If the given n, lies to the left 
of point C, then p = 0.5p, is preferable . In the zone 
from C to A, the optimum aumber of pole pairs is 
p = p;/1.5, for the given three values, Consequently, a 
certain value of p, for which the torque will be greatest, 
will correspond to each speed n,. Physically, the in- 
crease in the torque when the slip is small and p de- 
creases is explained by the increase in the penetration 
depth of the field, and consequently, by a lowering of 
the electrical resistance of the armature. The decrease 
in the torque, when the slip is large and p lessens, 
occurs because of an increase in the demagnetizing ef- 
fect of the armature reaction. 

In order to determine the optimum number of pole 
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pairs, relations are used which associate the magnetic 
quantities (the first harmonic of the field strength at the 
surface, Hr xo, and the flux density in the gap) with the 
torque M and the speed ng: 
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t 
Vi Anse = . a © 








(here p{m™) is the permeability determined from the 
basic magnetization curve when H = Hm xo), 


13,4Mkz . 


D4 peqBy sin (=) 
and also from the balance -of-mmf'‘s equation which 
follows from the vector diagram of the e.s.c. [7]: 


(10) 





k= 


(i,~)? = (11) 
= F§ (1 + x)® + kaF3 + 2Fs (4 + *) Fokaa sin 9. 
In expression (11), we get 
F;= B22 . (12) 
Ho 


for clutches with poles of alternate polarity; for induc- 
tion clutches: 


23 4Bk, 





F.= —= (13 
. Bu Bo Bo ) 
D 
F, 1 > quape (14) 
We. will also consider the relation 
A 
ug = KES, as) 


Eliminating Bg and Hmyxo from equations (9), (10), 
(12), and (14), and putting the values obtained for F¢ 
and F, into (11), and then examining the function M = 
= f(p),which results, at the extremes, we can get a value 
of py at which the torque is maximum. 

For clutches with poles of alternate polarity, 
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For induction clutches, n—1 
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The dependence of A on the exponent, pn, of the 
parabola which approximates the actual function, is 
shown in Fig. 3. Since the steel used for making the 
armature has n = 7 — 12, one can set A equal to 5, on 
the average. 

Comparing expression (16) for py) with (2) for the 
critical speed ng}, we find that the coefficient A is the 
ratio of the speed at which (for the given p = p,) the 
torque is greatest (in comparison with any other value 
of p) to the critical speed when p= py. Since A = 5, 
the maximum torque (which is the best operating point 
on the mechanical characteristic) occurs when ng © 
0.2ngk, and this gives M/My, * 0.7. The computation 
for 4 = const gives the best point when ng = 0.0607ngk 
(A = 16,5), and this corresponds to M/M, = 0.59. 

When the clutch is designed, the excitation ampere- 
turns iywy are initially unknown, and hence it is best to 
select them on the basis of the flux density in the air 
gap at no-load 

28 (1 + x) 


iw, = Bsmax — —* (17) 


In the determination of the main dimensions D and 
2 5, some nominal value of Bé is ordinarily given as the 
flux density in the air gap. The relation between Bs 
and Bs ax is given by the curve A = /(M/M,) in 
Fig. 1. The armature reaction mmf is low up to M/M,= 
=0.7, and the decrease in the flux density in the air gap 
is not more than 12% of that at no-load. Therefore 
Bs max = (1-0—1.14)Bs depending on the given My/M,. 

One further thing should be mentioned. If the over- 
load ability, My;/My = 0.7, meets technical require - 
ments, then ng should be substituted directly in the ex- 
pression for po. If My;/My must be less than 0.7, then 
the quantity € ngty should be used in (16), where 
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*"In the given expression, the numerical coefficient, 
13,4, takes into account the loss in the armature due to 
the higher harmonics of the flux density in the gap, 
whereas in [7], this coefficient, equal to 13.8, is derived 
without considering the latter. 
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e= = 
Ney / "x fora given M,/M, 
e 0,2 ’ 
~~ Nyy / Mg, fora given M,/M,, 


The ratio Ms, /n,, fora given M,/M, 
is determined from the mechanical characteristic (3). 
The curves given in [7] may be used for this purpose; they 
are plotted for various values of n, as a function of the 
given (3). 

Considering all this, we finally obtain: 
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for a clutch with poles of alternate polarity, and 
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for an induction clutch. 


2. The Procedure for Determining the 
Main Dimensions 

On the basis of the information given above, the 
following procedure is recommended for determining 
the main dimensions of an e.s.c, with a massive arma - 
ture, for given values of My; and nsH. 

1. The quantities D and 6 are chosen according to 
the practical demands of construction, and pp» is deter- 
mined, One can assume that Bs = 0.6 — 0.7 webers per 
m*, and t/r = 0.6 —0,75 for an iron-clad inductor. 

The drop in the mmf in the steel, F,, is evaluated 
while the steel is in an oriented state, Computations 
on actually constructed clutches show that on the aver- 
age, kK = 0.5 —1.0, The basic recommendations are 
given in [4] for induction clutches, 








2, The curve VpH =/(H). is plotted from the 


basic magnetization curve of the material chosen, 


3. The quantity Vm Himxq is computed from 


ele 


formula (9), and H,,,y9 is found from the curve 


Vp = /(H) 

4, The quantity 15 is determined from (10). 

5. The excitation ampere turns iyw, are found from 
(11), or can be found more accurately by the method 
given in [6}. 

Later, the geometric dimensions of the inductor 


are made more precise, and a verifying calculation is 
made, according to the method described in [6]. The 
dimensions D and 15, which are determined from elec - 
tromagnetic conditions, must be large enough to dissi- 
pate the heat caused by the losses in the anchor and ex - 
citation winding of the clutch, Since the heat transfer 
coefficient and the permissible overheating of different 
parts of the machine depend on a number of physical 
and constructional considerations, the problem of cal- 
culating D and Is from thermal conditions is very diffi- 
cult, Therefore, we recommend that experimental data 
on the determination of the permissible heat load, Aps, 
per unit area of the armature surface along the air gap, 
be used, 

In Fig. 4, curves are shown of J pg as a function of 
the linear velocity of the armature; they were plotted 
from the data of a number of studies on clutches with an 
external armature and air cooling, Curve 1 (for an ar- 
mature with a smooth external surface) and curve 2 
(for an armature with cooling fins) are plotted from the 
data of the MS [11]}, and curve 3 was plotted from 
data of the Stromag firm [10] for clutches with longi- 
tudinal channels along the external surface of the arma- 
ture, These clutches take an intermediate position be- 
tween clutches with a smooth armature and clutches 
with an armature having cooling fins, insofar as their 
ability to dissipate heat is concerned, Curve 3 also 
takes an intermediate position between the Enims 
curves, Experimental dependencies of 


mami (2) 


are presented in [9] for different power ratings of the 
actuating mechanisms; these dependencies are plotted 
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for the case of loss dissipation by the armature when the 
speed regulation is close, The relation 
kw + sec 


AP, max r 
——  =( 0.65 — 1.30 y 
met | i 


can be obtained from them. 
It follows from curve 3 of Fig. 4 that 


Ap, 
v, 





= (0.8 - 1.43) see, 


Hence, the experimental data presented mutually 
agree, and the curves in Fig, 4 can be recommended 
for direct calculations, 


3, Conclusions 

1, The formulas presented for determining the op- 
timum number of pole pairs and the main dimensions 
of electromagnetic slipping clutches with massive stee] 
armatures were derived from relations in which the 
change in permeability, 4, as the magnetic field 
strength in the anchor changes was taken into account, 
The correctness of these relations is confirmed by ex- 
perimental studies (6, 7]. 

2, When pp, D,and1¢ are found from relations 
based on the assumption that yp is constant [4, 5], the 
results depend on the average value of y selected, The 
method presented in this article gives a unique solution 
for given values of (My, men, M,/M],) . and the data 
obtained occupy an intermediate position in relation to 
the results obtained from the method explained in [5] 
(see example below). 

3, It has been theoretically established that induc - 
tion clutches have greater size than clutches of alter- 
nate polarity, when the same torque is transmitted by 
each, They are best used in small diameters (D = 10 
to 15 cm) and with an inside inductor, when the leak- 
age flux of an iron-clad inductor is significantly large 
in comparison to the effective flux, and, moreover, dif- 
ficulties arise in locating the excitation winding. Induc- 
tion clutches can be profitably used also when the di- 
mensions are determined by thermal conditions (for ex- 
ample, when the speed regulation is close), since their 
construction is simple. 

4, Average values of the permissible thermal load 
which can be used for practical calculations are present- 
ed on the basis of the analysis and generalization of 
data appearing in the literature, 





Example, Let us determine the main dimensions of an 
electromagnetic brake with an iron-clad inductor which 
develops a torque Mj; = 28 kg-m at ng = 120 rpm and 
M,;/M, = 0.7, The armature material is ST-20L, p = 
=Q16+10°¢ohm-m at 20°C. An approximation of the 
basic magnetization curve gives n = 10 and K = 0,675, 
The coefficient ¢ is determined from [12], We will 


take €= 1,7 as a first approximation, and an armature 

temperature of 120°C, Then €p oes = &p(1 +a At) = 0.408 x 

x 107° ohm-m, The curve v pH = f (H), plotted on 

the basis of the magnetization curve, is shown in Fig. 5, 
We will assume that 


D=0,4m,6=1,0mm, t/t=0,7; / 5 
= 0,65 webers/m?; x = 0,5. 
We find that kg = 


Nia 


) = 0,89. 


t 
= 0,85; 4 = 0,88; sin @ = 0,595; e = 4; sin (— 


Moreover, we will consider that 


? 
A = 5; k, = 0,845; ke = 0,87; pio = 1,255-10-* h/m. 


We find the optimum number of pole pairs from 
(18): p = 11,2, 
We will use p=12, We find from (9) that 


V oho? Ameo = 175 |/ 2 ajm 


We get Hm xe = 17500 a/m from the curve Vp H = 
= §{(H). From formula (10), we determine that 15 = 
= 0,2 m, Ultimately, the geometrical dimensions of 
the inductor are determined and a confirming calcula- 
tion is made, 

For comparison, we will compute p» according to 
T. A, Glazenko's method [5], There is a parameter 


Dp Tito, _ 165 
K= V us oo, = =. 
4de oy ’ y p 


Consequently, K depends on the chosen value of 
the relative permeability. If one assumes that p = 1000, 
as in the examples of [5], then K = 5,22 and pp = 6. 

This average value of » has no real basis, For ex- 
ample, in our case p. = 95.5 at the surface, If u = 150 
is taken as an average value, then K = 13,5, which gives 
Po = 12, 

Let us determine the main dimensions D and 1g by 
V. S, Sharov’s method for p = 12 from the relation pre- 
sented in [4] on page 18; 





Dl, = M pV up -10¥ ; cm’. 
7.1-Bl OV fy ky 
Here, [ p] is in h/cm, It follows that we take . 
p = 0,24°10™ ohm-cm, for our case: 





webers 
m2 





4 t x 
By = sin (+ 5) B, = 0.737 





= T3TagaUss/, = F7* -- 24 cps, 
k, = 1.2 ([4] ,page 77), 
D*l, = 7,3-408 Vp. 








In determining », V. S. Sharov recommends either 
that one assumes the value B A = (14 — 16)-10° gauss 
for the average magnetic flux density in the active layer, 
and finds p at this B , from the magnetization curve 
({[4}, pages 22, 76, 88), or that one is given By» as the 
flux density at the surface of the armature ((4], pages 
43, 61). 

We will use the first method, 

For a B Aequal to 14000 gauss, one can find that 
p = 0,76°10~° h /cm from the magnetization curve for 
ST-20L,andD? ,= 20,1-10° cm*, When D= 40 cm, 
lg = 31 cm, 

If B 4 = 16000 gauss, then 


yt = 0,296- 10-5; D%, = 12,5-108cm’, J, = 19,5 cm. 


Thus the results of the methods mentioned depend 
substantially on the value of y which is chosen, 
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Time relays which use the inertial properties of magnetic amplifiers are examined. It is known that under 
certain conditions, time relays with ashort-circuited network possess asatisfactorily high stability of operating 
time, An appoximate design method is presented for this type of relay. 


Contactless magnetic time relays are one of the 
varieties of magnetic amplifiers which operate linearly 
and generate a given time delay from the moment a sig- 
nal is received until the relay operates. 

The inertial properties of magnetic amplifiers are 
utilized in contactless magnetic time relays. Ordinari- 
ly, in using magnetic amplifiers, one attempts to lower 
their inertia as much as possible, since it is known that 
this is one of their greatest deficiencies, but special 
measures are taken to increase it in contactless magnet- 
ic time relays. 

It is known that the time characteristics of contact- 
less magnetic relays depend substantially on the feed- 
back coefficient, the amplitudes of the control and bias 
currents, and the circuit connection of theac winding. 
By using these factors, the delay time of ordinary con- 
tactless magnetic relays can be sharply increased (see, 
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for example [1-4] ). However, the length and especially 
the stability of the operating time of such devices is not 
always satisfactory, and this forces one to take special 
measures to increase the inertia of contactless magnetic 
relays. The following things are most frequently done 
to accomplish this; 

1) increase the inertia of the controlling network; 
this is done by connecting additional inductances, capa- 
citances, thermistors, etc. (Fig. la, 1b, 1c); 

2) variable negative feedback (Fig. 1d, le, 1f) 

3) short-circuited winding (Fig. 1g). 

Without going into detail on relays of the first type, 
we will note that the connection of additional elements 





* The work was reported at the All-union Conference 
on Magnetic Elements of Automation, Telemechanics, 
and Computer Technique, November 29, 1957. 
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in the controlling circuit means an increase in the con- 
trol power, and also in the size and weight. 

Here it is appropriate to mention a circuit suggested 
by V. N. Tkachev and S. V. Solntsev (Fig. 2a), The re- 
lay is activated by breaking contacts K, and K,. When 
this is done, the operating point moves from A to B (Fig. 
2b), and then, as the capacitor charges, the current in 
the control winding drops to zero, and the relay operates. 

The investigations of those authors have shown that 
the instability of the relay operating time does not ex- 
ceed 0.7% for a 10% change in the line voltage, and 
6.8% for a temperature change from +70 to -50°C, 

It should be mentioned, however, that the capacitance, 
C, is 15uf for an operating time of 3.24 seconds. 

Time relays using the transient in a variable negative 
feedback network have been examined in a number of 
papers, including those of I, B, Negnevitskii and I, S. 
Mikadze [6, 8]. 

We will note that besides such advantages of these 
relays as relatively large time delay, and the possibility 
of controlling the operating time over wide limits with- 
out significantly increasing the control power, there are 
also important disadvantages. In the first rank of these, 
one should mention the necessity of using high value 
capacitances (hundreds of microfarads for operating 
times of 5 to 10 seconds) or powerful differentiating 
transformers, which means a significant increase in size 
and weight. One should also mention the substantial de- 
pendence of the operating time on the line voltage fluc- 
tuations and control current fluctuations [6]. 

I. S. Mikadze suggested that the variable feedback 
signal be amplified by an additional magnetic amplifier 
in order to eliminate these disadvantages; this natural- 
ly complicates the circuit a great deal. 


1, A Contactless Magnetic Time Relay 





with a Short-Circuited Winding 











Contactless magnetic time relays with short-circuit- 
ed windings are very simple and reliable, and at the same 
time have small dimensions. However, the maximum 
time delay, and especially its stability, as a rule, is not 
high for relays of this type, and this limits their use. 

Practice has shown that time relays of this type can 
be successfully used when stability requirements are com- 
paratively stringent, provided the delay time is relative- 
ly small and special measures are taken for stabilization. 
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Fig. 2. 


Thus, for example, the relay developed by N. A. Polya- 
kova [7] ensures a time stability of + 5% for a line volt- 
age change of 4 10%, when the operating time is con- 
trollable from 5 to 10 seconds, 

The greatest error of the relay developed in the 
ITM Electric Model Shop and the VT of the Academy 
of Sciences, USSR [3], is 7% for an operating time of 3 
to 20 seconds, and a supply voltage change between 190 
and 240 volts (nominally, 220 volts). 

It should be noted, however, that the signal source 
voltage for this relay was stabilized by a special stabi- 
lizer, The greatest operating time instability of the re- 
lay was caused by temperature change, Unfortunately, 
no data on operating-time temperature stability were 
given in [3,7]. 

Let us examine the time relay whose circuit is 
shown in Fig. 3. The saturable reactor for this relay was 
made from alloy 79NM toroid cores, 42/30 mm in dia- 
meter, and 3 mm thick. 

Two windings, which were meant to increase the 
relay inertia, were wound on the saturable reactor, as 
well as the control, feedback, and bias windings. One 
of these windings (W),) was wound with 1 mm diameter 
wire and was short-circuited, and the other (w}3;) was 
wound with relatively thin wire and was shorted through 
a wire-wound rheostat, Ry, which permits one to control 
the operating time of the relay from 0.2 to 0.5 seconds, 

The rectifier bridge was assembled from type DGTs- 
24 germanium diodes, A small-sized type RSM-2 relay 
was used as an actuating relay. Resistance R, in the con- 
trol circuit, and Rg in the bias circuit serve for regula- 
tion. Resistance R, was made from copper wire, and is 
intended to compensate for temperature. The tempera- 
ture stabilization of the bias circuit is ensured by the 
fact that the bias winding is wound with Constantanwire. 

The results of the experimental investigation of the 
operating time stability of dozens of relays is given in 
the Table. 

The static characteristics of the time relay being 
studied are shown in Figs. 4, 5, and 6; they were taken 
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TABLE 





Greatest relay operating time change, in % 




















Test conditions of its nominal value 
t 0.5 | t =0.3 
Change of supply voltage 
(40 volts, 1000 cps) by +2 volts 1,9 2,7 
Temperature change from +20 
to+ 70°C 6 6,8 
Temperature change from + 20 
to - 50°C 6 3,5 
lyme 





6-6. 3% %. 1 2i,s0 





Fig. 5. 1) f= 1000 cps, U= 40v; 
Fig. 4. 2) f= 1010 cps, U = 40 v; 3) f= 
=990 cps, U = 39 v. 
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Fig. 6. Fig. 7. 
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at various voltages and frequencies of the source, and 
also at different temperatures. These characteristics 
were taken with the bias circuit open. 

The dependence of the operating time of the relay 
on the control current is shown in Fig. 7, as is also the 
reserve factor k,y = I¢/lop; the bias current was held con- 
stant. One's attention is drawn to the fact that as I, and 
consequently ky, grows, the dependence of the operating 
time on these quantities drops sharply. 

The dependence of the operating time on the bias 
current is given in Fig. 8, and the dependence of oper- 
ating time on control current is given in Fig. 9 for dif- 
ferent resistances connected to the wks; winding. 

The weight of the entire apparatus together with 
the actuating relay is 420 grams, The dimensions are 
55 x 80 X 72 mm. 


2. Methods for Raising the Stability of 





the Relay Operating Time 








It is evident that it is essential to ensure the stabili- 
ty of the relay static characteristics, primarily the pull- 
in current, in order to increase the operating time 
stability. Since the pull-in current is somewhat more 
stable than the drop-out current, it is essential that a 
contactless magnetic relay with normally open contacts 
be used. 

The basic way of stabilizing the static characteris- 
tics is to select the operating mode of the magnetic 
amplifier so that it approaches the ideal amplifier [9, 
11}. As is known, the latter is achieved when the follow- 
ing two conditions are met [9, 10, 11]. 

1) The total resistance of the load circuit, [R,, 
must be substantially less than the inductive reactance 
of the ac winding, X-h, in the no-load state. This con- 
dition will hold if Xch = (160— 200) ER,. 
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Fig. 8. 


2) The semiconductor rectifiers must be designed 
so that the reverse currents will be negligibly small, 
that is, Kp © 1/kg. Here Kp is the rectification coeffi- 
cient, and k¢ is the form factor. The latter conditions 
ordinarily occur when the reverse voltage on the recti- 
fiers is significantly less than the maximum permissible 
rating. 

If the above conditions are met, the necessary sta- 
bility of the relay operating current without bias is en- 
sured (see Figs. 4-6). 

The stability of the operating current often drops 
when bias is introduced, and one is obliged to take 
special measures to improve it, for example, stabiliza- 
tion of the bias voltage, temperature compensation of 
the bias circuit, etc, [1, 2]. 

Stabilization of the operating current against supply 
voltage changes can be accomplished with automatic 
bias, that is, bias supplied by a rectified current from 
the same source which feeds the ac relay circuit. 

As is known [1-3], the operating time stability of a 
contactless magnetic relay increases substantially as the 
operating current reserve factor increases (Fig. 7). How- 
ever, the increase in the reserve factor is limited by the 
power of the signal source, and by the relay properties. 

The choice of the feedback coefficient, Kq,, is 
equally important for the operating time stability of a 
contactless magnetic relay. 

It is known that the relay operating time sharply in- 
creases when Kgp is nearly unity, but at the same time, 
the instability both of the contactless relay static 
characteristic and of its operating (and drop-out) time 
also increases, The choice of the parameters of the 
short-circuited winding is very important. Indeed, the 
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temperature instability of the relay operating time is 
determined to a significant degree by the temperature 
instability of the resistance of this winding. It is not 
possible to make this winding from alloys with a low 
temperature coefficient of resistance because of the high 
resistivity. Therefore, the winding resistance must be 
made as low as possible. This can be done by increasing 
the wire diameter, It is evident that better utilization of 
the winding window is ensured as the wire diameter is in- 
creased, because the copper fill factor, k.,,, becomes 
large. 

3. Design of a Time Relay 





An approximate method for designing a contactless 
magnetic time relay with a short-circuited winding is 
presented below. 

The design procedure takes the following basic 
steps: 

a) select the actuating relay (repeater-relay), 

b) plot the static characteristic of the contactless 
magnetic relay 

c) design the contactless relay from its static 
characteristic, taking into account the operating time, 

d) design the mechanical construction of the 
saturable reactor. 

Ordinarily, the given quantities in the design of a 
contactless magnetic relay are: relay operating time, 
top, the number of contacts on the actuating relay and 
their basic characteristics, power source frequency, f, 
pull-in current, Iop, drop-out current, 1,,, and the am- 
bient temperature change 0-9. 

In choosing an actuating relay, one should be guided 
toward a relay with the least winding resistance compat- 
ible with minimum operating power. 

Starting with the operating (pull-in) current, I r,0p 
and drop-out current, I; oq of the actuating relay select- 
ed, and the given pull-in and drop-out currents of the 
contactless relay, its static characteristic is plotted. In 
doing this, the necessary reserves which ensure reliable 
operation and drop-out of the relay are provided (Fig. 
10). Ordinarily, IL,op= 1.25 Ir,op and I, do = 0.75 x 
X Te do* 

The bias is chosen on the basis of the preceding de- 
sign steps, and the schematic diagram of the contactless 
magnetic relay is drawn. 

The relay design procedure then follows the next 
series of steps, 

1, Rectifier Design 

a, The type of diodes is chosen from the quantities 


TL op & lot and UL op = I opi: then the number of 





diodes connected in parallel in one arm of the rectifying 
bridge is determined 


ree ~ (1) 
2lnc 


where I4, is the average amplitude of the permissible 
dc of one diode at maximum ambient temperature. 
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The number of diodes connected in series in one 
arm of the rectifying bridge is found; 
4,57] R 


Iv 


where Ur, is the amplitude of the reverse voltage on one 
rectifying element which is permissible at the maximum 
ambient temperature. 

b. The resistance of a bridge diode is found: 


Raf Un (3) 


Ra= eee ot eles 
d= 2m ( ~ + ies 
where Rg¢ is the differential diode resistance in the for- 
ward direction, and Uzp is the threshold diode voltage in 
the forward direction [10, 11]. 
2. Calculating the Parameters of the Saturable 
Reactor 
a. The inductive reactance of the saturable reactor 
without excitation is found: 


Xen = (160 - 200) ((1.2 - 1,4)R, +R). 





b, The necessary power supply voltage is computed; 


U.=1.11 Ty op I(1-2 — 1,4) Ry + Ral. (5) 


c. The volume of steel in the flux path is calcu- 


lated: t * . 
‘ Us hy 10% 


Vig mc O—Lis an. (6) 
st ; ’ 
{ Br Age ‘ch 


where jig is the differential permeability for the alter- 
nating component of the flux [10], B,, is the design 
value of {.ux density, gauss; and k,, is the packing fac- 
tor of the steel laminations making up the flux path 
(ordinarily, ks¢ = 0.85 — 0,9). 

At 
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+t The formula is derived in the Appendix, 
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One of the existing types of flux paths, for which 
Vst is greater than or equal to the value calculated from 
(6) is chosen. 

d, Thenthe number of turns inthe ac windings is com- 
puted:} 





wre 1,125 By Xop gt 
* Ha U (7) 
For a series connection of the ac windings, W~ = 2w~, 
and for a parallel connection W~ = ww. Here, wi is 
the number of turns on one core. 

e. The number of turns in the feedback winding is 
found; 


wp, = (0.505 ---0.525) W (8) 


for high-permeability permalloy alloys; and 


wp, = (0.525 — 0.65) W.. (9) 


for transformer steel flux paths, and also for high permea- 
bility alloys having an air gap. 
f. The no-load current is computed: 


U 


fa “Tk (10) 
g. The number of turns in the control winding is 


UW. = (WR — -— Lop — “nl = 
: - . top — fae | 





I —] 
- W.(Kp— 1) 22, (11) 


where A is the width of the relay characteristic. 

h. The no-load current corresponding to the given 
values Of Iop and Igo is calculated from the ly, go of the 
actuating relay which was chosen: 





| lon — lw 
& Kio c (12) 


Tay =I L,do ~ 
+r oo Ub 
If In) is greater than Ip), , a new value of Xop is 
computed; 


U 


Xa ile’ (13) 


Then a new value of V,, is computed, a new flux 
path is selected, and the number of turns in the ac, 
feedback, and control windings are recalculated, 


j. The ampere-turns of the bias winding are found 
[2]: 


Tyg Ms = Togpe— In Wf (14) 


Having been given I}, or Wps , we find the second 
factor. 

k. The diameter of the winding wires is calculated: 

ac winding 





I 
V 4 C, = 


where C = 1 for series connection, and C = 4 for parallel, 
control winding 


I, 
d, = 143// s. (16) 


The diameters of the wires in the feedback, dg, , 
and bias, d),,, windings are computed analogously. 

1, The area of the winding window, which is occu- 
pied by the ac , feedback, control, and bias windings, is 
found: 


Q, =< wate t metot edet webs 
oa 


, a7 
hea (17) 


where k,.,, is the copper-fill factor (ordinarily, koy = 
= 0.2 — 0.3). 

3. Designing the Short-Circuited Winding 

Starting from the given operating time, the area of 
the winding window occupied by the short-circuited 
coil is determined; * * 








t_(Kga—1) ple dy 10° 
Ce =e ee 


mee OR Aen (18) 
ip lao 


where p is the resistivity of copper (0.0175 ohm-mm" per 
m); / y is the average length of the wire in the winding, 
cm; } is the differential permeability for the magnetic 
field of the control signal (a method for determining pu 
is given in [1]), 


2wfh 








iy ae: (19) 
_%%c 20 
Qoe “Vie = 





t The formula is derived in the Appendix. 


** The formula is derived in the Appendix, 








Then the total window area occupied by all the 
windings, including the short-circuited one, is found. 
When this is done, an additional window area Q), is pro- 
vided to pass the reel of the winding machine: 


Qozr = Qy + Qos 7 Q,- (21) 


If the value of Qoy which is obtained exceeds the 
winding window area of the flux path chosen earlier, a 
flux path with a larger window is chosen, and the satur- 
able reactor must be redesigned. 

If the area Qos is less than or equal to the winding 
window area of the chosen flux path, the number of turns 
in the short-circuited winding is found from the formula, 


wo, — eos (22) 
$s d, 

The wire diameter, d,, is made as large as possible. 

4. Constructional Design of the Saturable Reactor 

The constructional design problem is one of deter- 
mining the size of the winding window necessary to 
contain all the windings and insulating layers, and com- 
paring it with the window of the type of flux path cho- 
sen, For toroidal reactors, the external diameter of the 
saturable reactor and its height are determined in addi- 
tion to this, A constructional design of a saturable 
reactor is given in [11]. 


APPENDIX 
1, Derivation of formula (6): 








0,4 x we Qst kee Ba 
2p 10° SCS 





X= ob = 2nf (1D 


Up =4.hh FW, Bin Qe ky 10% = (ID) 


Squaring both sides of equation (II), and solving for 
W2,, we get 


ose 


~ ; (1m) 
CPB Oe FP 


Substituting (III) into (I), we find 


ee 
ch Gy 10° G78, OF OY) 


Taking into account that Vt = Qs! ,,, we find 





Un, . 107 
f Binks X oh 





Vet =2 


after simple transformations. 
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2. Derivation of formula (7). Solving equation 
(II) for Qst, and substituting it into (I), we get 


2 
ch~ Tit, 6, 





whence 


o Hal, 





8. Derivation of formula (18). 
Ignoring the delay in the 4c circuit, one can write 


di U, 
t Cc +i — Cc 

ee ee $i gm , 

i—K,, @ R. (vd 
for an idealized amplifier [1], where r is the time con- 
stant of the magnetic amplifier control circuit with the 
feedback loop open. 

Taking the short-circuited winding into account 
(1) 





<5). wil 


t=te(14 
e\ | wes 


where R, is the resistance of the short-circuited winding, 
and R, is the resistance of the control winding. 














As is known, 
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Taking equations (IX-XI) into account, we will re- 
write formula (VII) in the form 


— wee Key 1 4 Qos 
10’ p EF ce (1+ G—) 


Ost Kee Koy Xx 
T= i055 hr Icy (Poe + es be (xm) 


Integrating equation (VI), and considering that 
i, = lop when * = t, and i, = Ido when t= &, we get 
(1) 


le J 


Te—! do 


t == . In 


7 Kp—t 





(XM) 


Substituting the value of tT from equation (XII), we 
find that 


_ Qt Kst%cu Poe t Qos 4, te = op : 


op “WW pllo, Ag—! 1, = Tgp 





whence 


(Bp — 1) e & Ing 10° 
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This article presents a method for determining the parameters of linear combined control systems and 
servosystems with respect to the transfer loci of the system to be controlled, the analytical expression of which 
can be unknown. 

















































































































The application of the disturbance-compensation In connection with this, the controlled variable in 
principle makes it possible to increase to a considerable the systems given inFig.la andb can be represented 
extent the accuracy of the operation of automatic con- by the following general expression: 
trol systems, and, in particular cases, even to secure the 
complete independence (invariance) of the controlled w= @ (p) [zy + WG, (Pp), test Mik (Pe td 
variable with respect to measured disturbances [1]. In 
contrast to ordinary systems (Fig. 1a) besides deviations (3) 
of the controlled variable, one or several additional where 
disturbances, suitably transformed in the introducing de- Wee (P) Me (Pp) 
vices, are fed to the regular input in combined control ow = 1 + Wee (p) We (p) 
systems (Fig. 1b). As in the case of ordinary systems, 
the aim in the design of combined systems is the deter- 
mination of the regulator and the disturbance introducing ~ 7 | Mec{P) 
device parameters for which the system sacrifices the a *cre cof % L 
performance specifications. In this, the practical reali- - 0) eA 
zation of the obtained transfer functions must be as sim- . KP) 
ple as possible. Xq 

The control systems in Fig. 1a andb can be for- b W.(p) 
mally transformed into the system shown in Fig. Ic, . . Woda) , 
where the disturbances are reduced to a single control < % (p) ¥ 4 Se 
input. It can be readily shown that the controlled varia- I 
ble x in this system will change in exactly the same way Sc ies 
as in the system in Fig. 1a,if the transfer function of the i Me (p) 
filter through which the disturbance is transmitted before 
it is introduced in the system as an additional control in X. Le X, £ 
put in the following form: = Mr(p) [>| Kale} - 

z 
~ Wea (p) Fig. 1. Structural diagrams of control systems. 
¢ (P)= ho) (1) X_) Command; xg)disturbance; xq) control 


input; x) controlled variable; ¢ ) deviation; 

Wsgc (p) and Wgq (p) are the transfer functions 

For the combined system in Fig. 1b, the filter-trans- of the system to be controlled with respect to 
fer function is determined by the equation, the control and disturbance inputs, respectively; 
Wr (p) is the regulator transfer function; W, (p) 
is the transfer function of the disturbance intro- 


Wea (P) 
wy p= 2. : ducing device; and Wf (p) is the filter-transfer 
‘ Wye (p) Wy (Pp) calle - function. 








(3) 


y3 
») 





is the servosystem transfer function with respect to the 
control input, and Wg, (p) is the k-th disturbance filter- 
transfer function. 

The more accurately the controlled variable responds 
to the control input and the lesser its response to disbur- 
bances, the better the performance of the control system. 
From this point of view, the part of the control system in 
Fig. lc which is characterized by the transfer function 
&p) must transmit the control input Xy with as little dis- 
tortion as possible. 

The output quantities of disturbance filters can be 
considered as interferences imposed on the useful control 
signal by the disturbances. Therefore, in adjusting the 
system, one should find such a filtertransfer function for 
which the best filtration of disturbances in the range of 
frequencies which are readily transmitted by the servo- 
system is secured. This condition is given by the expres- 
sion, 


| W¢ (iw) | = min, (4) 


i.e., the system response to disturbances decreases as the 
deviation from zero of the filter frequency response de - 
creases in the frequency range from w = 0 to the cutoff 
frequency W= Weo of the servosystem, which can be 
found from the frequency response curve | (iw) | (Fig. 
2.) 

The endeavor to satisfy optimum performance con- 
ditions (4) must not result in a reduction of the servo- 
system stability margin, which can be estimated with 
respect to the fluctuation characteristic, 


| O (iw,,) | 


j@uoy, = “ma. (5) 


Here, | ® (i0)| and | (iw,) | are the frequency 
response values for the servosystem at zero and resonance 
wy frequencies (Fig. 2), and Mma is the maximum allow- 
able value of the fluctuation characteristic (usually 
Mma = 1.4-2). 

The introduction of additional disturbances in the 
system shown in Fig. 1b does not affect the system stabil - 
ity, since the number of closed loops in the system does 
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Fig. 2. Frequency response curve of the 
servosystem. 


not increase, and the disturbance introducing devices 
do not enter the existing loops. This fact makes it 
possible to vary the filter-transfer function by changing 
the transfer function of the disturbance -introducing de- 
vice without causing the system to lose its stability. 
Consequently, the design of a combined system can be 
divided into two stages. 

1. The regulator is adjusted under the assumption 
that the system operates only in response to the controlled 
variable deviation, ‘rhis is necessary because it is 
impossible to introduce additional effects from all dis- 
turbances in the system. 

2. The transfer function of the introducing devices 
is determined with respect to condition (4). 

The first-stage calculation is performed in the usual 
manner [2 and 3]. By taking into account Eq. (2), the 
optimum performance condition (4) is rewritten in the 
following manner: 


Wy (ie) Ww. (ie) Wi (i@)| = min. (6) 





In the particular case, where it is possible to secure 
the following form of the introducing-device transfer 
function 


: Wi (Pp) D 
i O- WeDo’ ' 





the controlled variable is not at all affected by the dis- 
turbance. It should be noted that in order to satisfy the 
absolute invariance conditions, the transfer function of 
the introducing device must coincide with the transfer 
function (1) of the filter in a system which responds only 
to controlled variable deviations. 

In actual automatic control systems, it is usually 
difficult to satisfy the absolute invariance conditions, 
since the technical realization of the introducing-~device 
transfer function, which is necessary for satisfying the 
above conditions, is either difficult or physically impos- 
sible. 

In this case, one should select a simple introducing 
device, the transfer function of which, in a certain 
sense, would deviate as little as possible from the transfer 
function of an ideal introducing device. In order to 
solve this problem, in particular, we can use the simplest 
approximation method, the interpolation method, accord- 
ing to which it is necessary that the functions to be ap- 
proximated coincide in several points assigned before - 
hand. In practice, in order that condition (4) is suffi- 
ciently well satisfied, it is usually sufficient thar the 
transfer loci Wy (iw) and Wij (iw): 


Neg (0) 
Wye (0) W, (i0)" 


ia. 
Wee (io, Wy (iq) = (8) 


Wy (0) = Wig ©): 


Wy (i@y) = Wig (ia) 
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Fig. 3. Configuration of the transfer locus vectors for an ideal introducing 
device and a circuit of passive electric fourpoles for securing approximate 
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coincide for the zero and resonance frequencies of the 
system. 

The above conditions can be realized by means of 
relatively simple devices. Figure 3 shows several most 
typical configurations of the W,; (iw) locus vectors for 
w= 0 and w= w,. The same Figure also shows examples 
of passive electric fourpoles for which conditions (8) can 
be realized for each of these cases. 

The tuning conditions for a combined system with 
additional control input can be found in a similar man- 
ner (Fig. 4). The mismatching in this system is deter - 
mined by the equation 


i , 
e=xy—2 | (py (et |, (9) 


and, therefore, the optimum performance condition is 
written in the following manner: 


1 
Wee (p) Wr (p) P, 





Wiz (p) = min. (10) 


The considered method for determining combined 
system parameters can also be extended to the design of 
autonomous control systems with several interdependent 
variables. 


Example of Calculation 

An automatic control system has two additional in- 
puts; a disturbance and a control input (Fig. 5). The 
transfer functions of the system,to be controlled with 
respect to the control and the disturbance inputs,are de - 
termined by the equation: 





, _ 1‘ 1 
The regulator transfer function is of the form: 


i 
Trep 





\ 
W, (p) = hy we ). 
where k, and Ty. are the transfer coefficient and the re- 
set time, respectively. 
The calculation of the regulator tuning, performed 


by means of the method described in [3] for Mma = 1.62, 
yielded the following results: 


= 11, = 4.5min, @, = 2.8 min -!. 
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1. Determination of the transfer function of the 
disturbance -introducing device. By substituting in Eq. 














(7) the expressions for the transfer functions of the system 





to be controlled and the regulator, we obtain: 


Waig (Pp) = 0.004 rer 
i.e., the actual differentiating element makes it possible 
to satisfy the absolute invariance conditions for the con- 
trolled variable with respect to disturbance. 

2. Determination of the transfer function of the con- 
trol input introducing device. In correspondence with 
Eq. (7), the transfer function of the introducing device, 








p?+ 3p*+ p 


=0, 136 T5p+i 


Wo (p) 
for which the absolute invariance conditions are satisfied, 
cannot be physically realized. 

The vector 


Wie) — 0.935et42.6" 

is in the second quadrant of the complex plane, and, 
therefore, for approximate invariance, an introducing de- 
vice of the type shown in Fig. 3c must be chosen. If 
C,/C, = 2 in this system, the maximum ratio of time 
constants T,/T, will be equal to 0.41. In this case, the 
numerical values of the transfer coefficient kj and of the 
time constants T, and T, are determined from the equa - 
tions 


Ti, 
i —-Tle 





tg Pi (o,) = 


7? oF 


A — ’ 








which yields 


ks == 12.9; T, = 0.252; 7's -= 0.104. 
Figure 6 shows the variation of the controlled vari - 
able in a combined system for the adjustment found and 


a unit-step input. For the sake of comparison, this Figure 
also shows the variation of the controlled variable under 
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the same conditions, where, however, no additional con- 
trol signal input is provided (dotted line). 
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The equation for the hydraulic servomotor in automatic control systems with a jet amplifier is derived. 
Experimental dataderived in verifying the obtained result are provided. 


A hydraulic servomotor, which is controlled by a 
jet amplifier, is used as the final mechanism in many 
automatic control systems. Various designs of amplifi- 
ers with a flapper or a jet pipe as mobile elements are 
used, 

In certain systems, a reset device is provided in one 
of the piping branches connecting the servomotor to the 
amplifier, As the working liquid moves through the 
throttle of the reset device, a pressure drop arises across 
its membrane causing the membrane to displace. 

Thus, a considerable number of control systems are de- 
signed according to the scheme shown in Fig, 1. 

We shall determine the equations of motion of the 
servomotor in such automatic control systems, 

Every jet amplifier design is characterized by a 
well-defined static dependence of the mobile element 
position on the pressure drop in the receiving nozzles, 
This dependence is such that a deflection AS of the flap- 
per or the jet tube from the equilibrium position that 
does not exceed a certain given value AS» causes an al- 
most proportional increase in the pressure drop OP be- 
tween the receiving nozzles of the jet amplifier. 

Consider the system shown in Fig. 2, We shall neg- 
lect the effect of external loads on the servomotor mo- 
tion, the friction in the servomotor, its unbalanced state, 
the inertia of the servomotor mobile parts, and water 
leaks in the servomotor, as is often done in deriving the 
equations for a servomotor which is controlled by a cut- 














Fig. 1. (1) System to be con- 
trolled; (2) measuring com- 
ponent; (3) amplifier; (4) re- 
set device; (5) servomotor, 


off valve. Assuming that the water flow in the piping 
connecting the amplifier to the servomotor is turbulent, 
we obtain the following expression for the pressure drops, 
which are measured by pressure gauges, 


AL\ yu? 


where ¢, is the resistance coefficient of the servomotor 
throttle, A is the linear resistance coefficient, J is the 
piping length, d is the inside diameter of the piping, y 
is the working liquid specific weight, g is the accelera- 
tion due to gravity, v is the water flow velocity in the 
piping, and AP, is the pressure drop across the reset-de- 
vice membrane, 

From the equation of continuity, it follows that the 
water flow through any piping section is equal to the 
volume traversed by the piston in unit time. 














Fig. 2. (1) Servomotor; (2) servomotor 
throttle; (3) pressure gauge; (4) amplifi- 
er; (5) reset-device throttle; (6) reset de- 
vice. 








Consequently, 





= (2) 


where m¢ is the servomotor piston velocity, F, is the ac- 
tive surface area of the servomotor piston, and f is the 
inside cross-sectional area of the piping. 

From (1) and (2), it follows: 


a= He D(a 0 


For the servomotor motion, oo < AP always holds. 
In this case, it can be considered that the 4P— AP 


difference is equal to the pressure loss in the amplifier 
receiving nozzles, which is a quadratic function of the 
servomotor velocity, i.e., 
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AP -- AP,=t—3,— » (4) 


where 6, is the resistance coefficient of the amplifier 
receiving nozzles, 

This condition will always be satisfied if we assume 
that the , coefficient is not constant, but depends on a 
given manner on the parameters that define the state of 
the system in Fig. 2, Such parameters are obviously the 
position AS of the jet pipe, the static pressure drop AP, 
and the liquid flow resistances which are determined by 
the servosystem design. If we neglect the inertia of mo- 
bile parts in the reset device, the friction in it, and the 
initial tension in its spring, the effect of which is felt 
only for very small deviations of the mobile element 
from its neutral position, we can establish a well-defined 
relation between the pressure drop across the reset-de- 
vice membrane and its displacement: 


(5) 


where c; is the over-all rigidity of all elastic elements 
acting on the reset-device stem, which is reduced to the 
stem, Am, is the reset-device-stem displacement, and F; 
is the active reset-device membrane area. 

The following relation holds within the linear inter- 
val of the jet amplifier characteristic: 


AS 
AP - AP, AS, ’ (6) 


where AP, is the static pressure drop between the receiv- 
ing nozzles, which corresponds to the mobile element 
deflection AS». 


By taking into account (3), (5), and (6), equation (4) 


can be rewritten thus: 
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By introducing relative coordinates, 
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where Am, is the servomotor displacement, m, is the 
maximum servomotor stroke, and m,; is the maximum 
stroke of the reset device stem, we obtain: 
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Finally, by using (8), the servomotor equation in the 
considered case can be written in the following form: 
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When o changes from 0 to 1, T, does not remain con- 
stant. This is explained by the fact that the expression 
determining T, contains the ¢, coefficient, which, as was 
established experimentally, basically depends on 0. As 
the jet amplifier mobile element moves from the neutral 
to the extreme position, the change in Ty is determined 
by the ratio of €. to the over-all resistance coefficient 
over the path traveled by the working liquid. If the ¢ 
coefficient constitutes only an insignificant part of the 
over-all resistance coefficient ¢, it can be considered 
that T, = const throughout the entire range of the relative 
coordinate © variation. 

In this case, the dependence of the servomotor veloc- 
ity pu” on O— pE will be close to a parabolic dependence. 
If the value of €,, is larger than all the other resistances, 
the dependence of T, on o cannot be neglected. It fol- 
lows from experimental data that, as the mobile element 
deflection from its neutral position increases, the value 
of ¢ decreases. This means that, for o =0, the value of 
T, Will be at the maximum. Due to this, the deviation 
from linearity of the u = f(o— p€) characteristic will be 
less pronounced, 








This deduction can be confirmed experimentally. 
Figures 3 and 4 show the results obtained in experimental 
verification of a system which differs from that in Fig. 2 
by the fact that it does not contain the reset device 6 
with throttle 5 (i.e., the equation p& = 0 holds). 

During the tests, the jet-tube deflections did not ex- 
ceed the ASp value, In the case shown in Fig. 3, there 
was practically no resistance in the servomotor throttle. 
Due to this, the time for the servomotor displacement 
through the distance m, was 19 sec for the maximum jet- 
pipe deflection, It follows from the above graph that 
the servomotor motion can be sufficiently accurately de- 
scribed by the equation 


’ (10) 


where T; = const, 
In the second case (Fig. 4), there was a considerable 
resistance to the liquid flow from the amplifier to the 
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servomotor, due to which the minimum time for the ser- 
vomotor displacement through distance m, was 160 sec. 
In the given case, the servomotor motion can be suffi- 
ciently accurately described by the equation 


p’=signo vie J ; (11) 





where T, = const. 

Thus, in investigating the dynamics of an automatic 
control system with a hydraulic servomotor which is driv- 
en by a jet amplifier, the servomotor equation should be 
written in the form given by (9). 

The T, coefficient in this equation generally depends 
on the relative coordinate o, the pressure drop AP, be- 
tween the nozzles, which corresponds to the mobile ele- 
ment displacement AS», and the flow resistance of the 
working liquid, This dependence must be determined 
experimentally for every individual jet amplifier type. 
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